Thermosciences  Division 


MECHANISMS  OF  IONIZATIONAL  NONEQUILIBRIUM 

IN  AIR  PLASMAS 


Final  Technical  Report  for  the  Period: 

May  1,  1997  -  September  30,  2002 
Submitted  to 
Dr.  Robert  J.  Barker 
Air  Force  Office  of  Scientific  Research 
Multidisciplinary  Research  Program  of  the  University  Research  Initiative 

Grant  No.  F49620-97- 1-0316 

Submitted  by 

Professor  Charles  H.  Kruger 
Principal  Investigator 

October  2002 


Mechanical  Engineering  Department 
Stanford  University 
Stanford,  California  94305 


REPORT  DOCUMENTATION  PAGE 


AFRL-  SR- AR-TR-02- 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  e: 
the  collection  of  information.  Send  comments  regardinrrztNs  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions 
Operations  and  Reports,  T  21 5  Jefferson  ftvis  flight/,  Site  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  R 


1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 


01  May  97  to  30  Sep  02  Final 

4.  TITLE  AND  SUBTITLE  5.  FUNDING  NUMBERS 

MURI  97  Mechanisms  of  Ionizational  Nonequilibrium  in  Air  Plasmas  61 103D 

3484/RS 


i 


6.  AUTHOR(S) 

Dr  Kruger 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Stanford  University 
651  Serra  Street 
Room  260 

Stanford  CA  94305-4125 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Department  of  the  Air  Force 
Air  Force  Office  of  Scentific  Research 
801  N.  Randolph  St  Rm  732 
Arlington.  VA  22203-1977 

11.  SUPPLEMENTARY  NOTES 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

F49620-97-1-0316 


12a.  DISTRIBUTION  AVAILABILITY  STATEMENT  12b.  DISTRIBUTION  CODE 

Distribution  Statement  A.  Approved  for  public  release;  distribution  is  unlimited. 


13.  ABSTRACT  (Maximum  200  words) 

The  goal  of  this  MURI  program  was  to  discover  physical  mechanisms  for  significantly  reducing  the  power  budget  required  to 
create  and  sustain  large  volume  atmospheric  air  plasmas  at  temperature  below  2000K  with  electron  number  densities  of  the 
order  of  10  13  cm  -3.  The  Stanford  consortim  Experimentally  demonstrated  the  generation  of  over  10  12  electrons/cm  3  in 
atmospheric  pressure  air  at  2000K,  with  an  average  power  of  12  W/cm  3,  i.e.  250  times  smaller  than  a  DC  discharge  at  10 
12/cm3. 


14.  SUBJECT  TERMS 


15.  NUMBER  OF  PAGES 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 


19.  SECURITY  CLASSIFICATION  20.  LIMITATION  OF  ABSTRACT 

OF  ABSTRACT 


Unclassified 


Unclassified 


Unclassified 


UL 

Standard  Form  298  (Rev.  2*89)  (EG) 

Prescribed  by  ANSI  Std.  239.18 

Oesigned  using  Perform  Pro,  WHS/DIOR,  Oct  94 


/oG?  /  O  cf 


MECHANISMS  OF  IONIZATIONAL  NONEQUILIBRIUM 

IN  AIR  PLASMAS 


Final  Technical  Report  for  the  Period: 

May  1 ,  1 997  -  September  30,  2002 
Submitted  to 
Dr.  Robert  J.  Barker 
Air  Force  Office  of  Scientific  Research 
Multidisciplinary  Research  Program  of  the  University  Research  Initiative 

Grant  No.  F49620-97-1-0316 

Submitted  by 

Professor  Charles  H.  Kruger 
Principal  Investigator 


October  2002 


Contents 


2 


Section  Page 

Contents . 2 

1.  Abstract . 5 

2.  Introduction . 7 

3.  Thermal  Baseline  Experiments . 9 

3.1.  Plasma  Torch  Facility . 9 

3.2.  Plasma  Characterization . . . 11 

3.3.  Power  Budget . 17 

3 .4.  Conclusions  on  Baseline  Experiments . 1 8 

4.  DC  glow  discharges  in  atmospheric  pressure  air . 19 

4.1.  Introduction . 19 

4.2.  Two-temperature  kinetic  simulations . 20 

4.2 . 1.  Two-temperature  kinetic  model . 20 

4.2.2.  Results . 23 

4.2.3.  Analysis  of  the  ionization  mechanisms . 26 

4.2.3. a.  Ionization  mechanism  in  Region  A . 27 

4.2.3.  b.  Ionization  mechanism  in  Region  B . 27 

4.2.3.  c.  Analytical  solution . 28 

4.3.  Predicted  Electric  Discharge  Characteristics . 31 

4.3.1 .  Ohm's  law . 31 

4.3.2.  Electron  energy  equation . 32 

4. 3.2.  a.  Introduction . 32 

4.3.2. b.  Rate  coefficients  controlling  the  vibrational  distribution  ofN2  levels . 34 

4.3.2. C.  Vibrational  distribution  of  N2  ground  state . 40 

4.3. 2. d.  Inelastic  electron  energy  losses  in  air  plasmas . 42 

4.3.2. e.  Predicted  DC  Discharge  Characteristics  in  Atmospheric  Pressure  Air . 44 

4.4.  Experimental  DC  glow  discharges  in  atmospheric  pressure  air  plasmas . 45 

4.4.1.  Introduction . 45 

4.4.2.  DC  discharge  experimental  setup . 46 

4. 4. 3.  Spectroscopic  Measurements . 48 

4. 4. 3.  a.  Measurements  without  DC  discharge  applied . 48 

4.4. 3.  b.  Measurements  with  DC  discharge  applied . 51 

4.4. 4.  Current  Density  Measurements . 56 

4. 4. 5.  Electric  Field  Measurements . 57 

4.5.  Electrical  Characteristics  and  Power  Requirements  of  DC  Discharges  in  Air . 58 


3 

4.6.  DC  Discharges  in  Nitrogen . 62 

4.7.  Conclusions . 63 

5.  Repetitively  Pulsed  Discharges  in  Air . 65 

5.1.  Introduction . 65 

5.2.  Experiments  with  a  Single  Pulse . 67 

5.3.  Experiments  with  100  kHz  Repetitive  Discharge . 68 

6.  Computational  Analysis  of  Diffuse  Discharges  in  Atmospheric  Pressure  Air . 74 

6.1.  Introduction . 74 

6.2.  Flow  Field  and  Discharge  Model . 74 

6.3.  Numerical  Method . 77 

6.4.  Results . 80 

6.5  Approximate  Discharge  Model . 82 

6.6.  Conclusions . 88 

7.  Measurements  of  Ion  Concentration  by  Cavity  Ring-Down  Spectroscopy . 90 

7.1.  Introduction . 90 

7.2.  Cavity  ring-down  spectroscopy . 91 

7.3.  N2+  Measurements . 92 

7.3 A.  Experimental  Setup . 92 

7.3 . 1  .a.  Atmospheric  Pressure  Discharge . 92 

7.3. 1. b.  CRDS  Measurements . 93 

7.3. 1. c.  Electrical  Measurements . 95 

7.3.2.  N2+  Measurements  -  Results  and  Discussion . 96 

7. 3. 2.  a.  N2+  Ring-Down  Spectra . 96 

7.3. 2. b.  Spatial  Profiles  of  Ion  Concentration  and  Electron  Number  Density . 98 

7.3.2.  c.  Temporal  Profiles  of  Ion  Concentration  and  Electron  Number  Density . 100 

7.3. 2. d.  Discharge  Nonequilibrium . 103 

7.4.  NO+  Measurements . 104 

7.4.1.  Experimental  Setup . 104 

7.4.1.  a.  RF  Air  Plasma . 104 

7.4. 1. b.  CRDS  Measurements . 104 

7. 4. 2.  Results  and  Discussion . 1 06 

7.5.  Conclusions . 107 

8.  Summary  and  Conclusions . 109 

9.  Personnel . Ill 

10.  Refereed  Publications . 112 

11.  Interactions/Transitions . 113 


4 


11.1.  Participations  at  meetings,  conferences,  seminars . 113 

1 1.2.  Consultative  and  Advisory  Functions  to  Other  Laboratories  and  Agencies . 116 

12.  New  Discoveries,  Inventions,  or  Patent  Disclosures . 116 

13.  Honors/Awards . 116 

14.  References . 119 


5 


1.  Abstract 

The  present  report  describes  the  accomplishments  of  the  Stanford  Consortium  during  the  Air 
Plasma  Rampart  MURI  program  (May  1997-September  2002).  The  Stanford  consortium  includes 
the  groups  of  Profs.  Kruger  and  Zare  at  Stanford  University,  Prof.  Candler  at  the  University  of 
Minnesota,  and  Dr.  Kelley,  at  Boeing  until  April  2001  and  then  with  the  Stanford  group  until  Sep¬ 
tember  2002.  The  goal  of  this  MURI  program  was  to  discover  physical  mechanisms  for  signifi¬ 
cantly  reducing  the  power  budget  required  to  create  and  sustain  large  volume  atmospheric  air  plas¬ 
mas  at  temperature  below  2000  K  with  electron  number  densities  of  the  order  of  1013  cm'3.  The 
Stanford  consortium  has  achieved  the  following  results: 

•  Experimentally  determined  the  power  requirements  of  thermal  air  plasmas  with  electron  densi¬ 
ties  higher  than  1013  cm'3. 

•  Determined  the  chemical  reaction  pathways  that  control  ionization  and  recombination  in  two- 
temperature  atmospheric  air  and  nitrogen  plasmas  produced  by  DC  and  pulsed  discharges. 

•  Experimentally  demonstrated  the  production  of  stable  diffuse  glow  discharges  in  low  tempera¬ 
ture  (<2000  K)  atmospheric  pressure  air  and  nitrogen  plasmas  with  electron  number  densities  in 
excess  of  1012  cm"3,  more  than  six  orders  of  magnitude  higher  than  in  thermally  heated  air  at 
2000  K.  The  measured  discharge  characteristics  compare  well  with  the  predictions  of  a  two- 
temperature  kinetic  model.  Experimental  and  modeling  results  show  that  the  steady-state  elec¬ 
tron  number  density  exhibits  an  S-shaped  dependence  on  the  electron  temperature,  a  behavior 
resulting  from  competition  between  ionization  and  charge  transfer  reactions.  Non-Maxwellian 
effects  are  shown  to  be  unimportant  for  the  prediction  of  steady-state  electron  number  densities 
at  least  for  nitrogen  plasmas. 

•  Experimentally  demonstrated  the  generation  of  over  1012  electrons/cm3  in  atmospheric  pressure 
air  at  2000  K,  with  an  average  power  of  12  W/cm3,  i.e.  250  times  smaller  than  a  DC  discharge 
at  1012/cm3.  This  result  was  achieved  with  short  repetitive  high-voltage  pulses  (100  kHz,  12 
kV,  10  ns)  following  a  technique  originally  proposed  in  our  1999  Annual  Technical  report. 
Plasma  volumes  of  about  1  cm3  were  achieved  using  a  multi-electrode  configuration. 

•  Developed  spatially  and  temporally  resolved  diagnostics  of  N2+  ion  concentrations  in  DC  and 
pulsed  atmospheric  pressure  nitrogen  discharges  by  means  of  Cavity  Ring-Down  Spectroscopy 
(CRDS).  Spatial  resolutions  of  the  order  of  0.5  millimeter  and  temporal  resolutions  of  1  micro¬ 
second  were  achieved.  We  have  also  conducted  preliminary  CRDS  experiments  to  measure  the 
NO+  ion  concentration  in  air  plasmas  via  its  infrared  transitions  at  4.1  micron.  Interpretation  of 
the  measurements  is  found  to  be  limited  at  present  by  a  poor  knowledge  of  the  NO+  ion  spec¬ 
troscopy. 
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•  Developed  a  self-consistent  computational  fluid  dynamics  model  with  nonequilibrium  thermo¬ 
chemistry  of  the  two-temperature  discharge  experiments.  The  simulations  examined  two- 
dimensional  effects  associated  with  the  practical  realization  of  low  temperature  air  plasmas  with 
high  electron  concentrations. 
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2.  Introduction 

Nonequilibrium  discharges  in  air  and  nitrogen  plasmas  at  atmospheric  pressure  present  con¬ 
siderable  interest  for  a  wide  range  of  applications.  Desirable  conditions  are  electron  densities  of  the 
order  of  1013  cm'3  at  gas  temperatures  less  than  2000  K.  The  equilibrium  concentration  of  electrons 
in  low  temperature  atmospheric  pressure  air  is  extremely  small.  At  2000  K,  for  instance,  it  is  ap¬ 
proximately  equal  to  3.3xl06  cm'3.  To  produce  and  sustain  electron  number  densities  of  1013  elec¬ 
tron/cm3  in  air  while  maintaining  the  gas  temperature  below  2000  K,  energy  must  be  added  to  the 
plasma  in  a  targeted  fashion.  One  method  is  to  apply  the  energy  addition  to  the  free  electrons  by 
means  of  an  electrical  discharge.  Under  certain  conditions,  the  energetic  electrons  produced  by  the 
discharge  can  catalyze  plasma  ionization  reactions  and  inhibit  electron  recombination  channels 
without  significantly  heating  the  gas.  Power  budget  reduction  strategies  can  then  be  designed  by 
enhancing  ionization  channels  or  by  exploiting  electron  recombination  processes.  It  is  therefore 
important  to  understand  the  chemical  mechanisms  of  ionization  and  recombination  in  air  plasmas 
with  energetic  electrons.  Such  plasmas  are  sometimes  referred  to  as  two-temperature  plasmas,  with 
the  kinetic  temperature  of  the  free-electrons  (Te)  being  higher  than  the  kinetic  temperature  of  heavy 
species  (Tg). 

Section  3  of  this  report  describes  a  baseline  experiment  to  determine  the  power  required  to 
produce  1013  electrons/cm3  with  a  thermal  plasma.  In  this  case,  the  plasma  is  close  to  equilibrium 
and  the  gas  temperature,  about  4500  K,  is  higher  than  the  targeted  range.  Nevertheless,  measure¬ 
ments  of  the  power  requirements  provide  a  useful  point  of  comparison  for  power  reduction  strate¬ 
gies. 


Section  4  describes  experiments  and  theoretical  analyses  of  nonequilibrium  DC  discharges 
in  atmospheric  pressure  air  plasmas  at  temperatures  around  2000  K.  We  present  first  our  efforts  to 
develop  a  two-temperature  chemical  kinetic  model  of  air  plasma  discharges.  We  then  discuss  our 
DC  discharge  experiments  in  atmospheric  pressure  air  and  nitrogen.  We  conclude  this  section  with 
measurements  and  predictions  of  electrical  discharge  characteristics  (electric  field  vs.  current  den¬ 
sity)  and  power  requirements  (product  of  the  electric  field  and  current  density).  The  influence  of 
departures  from  a  Maxwellian  distribution  of  the  free-electrons  are  examined  by  means  of  a  simpli¬ 
fied  collisional-radiative  model  coupled  with  a  Boltzmann  solver  for  a  nitrogen  plasma. 

Section  5  describes  results  from  pulsed  discharge  experiments  in  which  electron  number 
densities  of  more  than  1012  cm"3  in  air  are  produced  with  approximately  12  W/cm3,  more  than  two 
orders  of  magnitude  lower  than  the  power  required  for  a  DC  discharge.  The  basis  of  the  technique 
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is  to  apply  short  (ten  nanosecond),  high  voltage  (~10  kV)  electric  pulses  with  a  repetition  frequency 
tailored  to  match  the  recombination  time  of  electrons. 

Section  6  presents  a  computational  methodology  to  analyze  the  DC  and  pulsed  discharge 
experiments.  A  complete  axisymmetric  simulation  including  an  11 -species  finite-rate  air  kinetics 
model  and  a  three-temperature  energy  formulation  was  used.  The  three-temperature  model  is  used 
to  allow  separate  heavy-particle,  vibrational  and  electron  translation  temperatures  of  the  plasma.  A 
self-consistent  effective  binary  diffusion  model  is  used  for  the  molecular  transport  properties.  The 
discharge  is  modeled  as  a  conducting  channel,  allowing  a  direct  enforcement  of  total  current  con¬ 
servation. 

Section  7  describes  experiments  conducted  to  measure  spatial  and  temporal  profiles  of  N2+ 
ion  concentrations  in  DC  and  pulsed  atmospheric  pressure  nitrogen  glow  discharges  by  Cavity 
Ring-Down  Spectroscopy.  Using  a  collisional-radiative  model  we  infer  electron  number  densities 
from  the  measured  ion  profiles.  The  values  of  electron  number  density  found  in  this  way  are  con¬ 
sistent  with  those  found  from  spatially  integrated  electrical  measurements.  The  spectroscopic  tech¬ 
nique  is  clearly  favorable,  because  it  offers  spatial  resolution  and  does  not  require  knowledge  of 
other  discharge  parameters.  The  measurements  also  demonstrate  the  feasibility  of  repetitively  puls¬ 
ing  the  discharge  as  a  means  to  elevate  the  time-averaged  ionization  fraction.  Measurements  of  the 
NO+  ion  in  air  plasmas  have  also  been  demonstrated.  The  accessible  spectral  features  of  NO+  are 
vibrational  transitions,  are  considerably  weaker  than  the  UV  electronic  transitions  used  to  probe 
N2+.  Nevertheless,  CRDS  data  from  air  plasmas  was  obtained,  and  spectral  features  attributed  to 
NO+  were  observed.  More  precise  knowledge  of  the  spectroscopic  constants  of  NO+  is  needed  to 
match  the  experimental  data  with  modeled  spectra. 


Section  8  provides  conclusions  and  summary  of  this  report. 
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3.  Thermal  Baseline  Experiments 

Thermal  plasma  processes  represent  a  convenient  method  for  producing  atmospheric  pres¬ 
sure  air  plasma  volumes  of  dimensions  greater  than  5  cm  in  all  directions  with  the  desired  electron 
number  densities.  Although  the  plasmas  generated  with  thermal  processes  are  in  a  state  of  Local 
Thermodynamic  Equilibrium  (LTE),  which  implies  that  the  gas  temperature  (>-4500  K)  is  higher 
than  the  targeted  range  for  applications  of  interest,  it  is  interesting  to  establish  the  power  budget 
required  for  generating  such  plasmas.  This  section  of  the  report  describes  the  thermal  baseline  ex¬ 
periments  conducted  at  the  High  Temperature  Gasdynamics  Laboratory  of  Stanford  University.  In 
this  work,  the  generation  of  steady-state,  atmospheric  pressure  air  plasma  volumes  greater  than  5 
cm  in  all  directions,  with  electron  number  densities  greater  than  1013  cm'3,  was  demonstrated  and 
the  power  budget  requirements  were  measured. 


3.1.  Plasma  Torch  Facility 

Atmospheric  pressure,  thermal  air  plasmas  were  generated  by  a  50  kW  inductively-coupled 
plasma  torch  (TAFA  Model  66)  powered  by  a  radio-frequency  LEPEL  Model  T-50-3  power  supply 
operating  at  4  MHz.  The  power  supply  delivers  up  to  120  kVA  of  line  power  to  the  oscillator  plates 
with  a  maximum  of  12  kV  DC  and  7.5  A.  The  oscillator  plates  have  a  maximum  RF  power  output 
of  50  kW.  A  schematic  drawing  of  the  plasma  torch  head  is  shown  in  Figure  1.  The  feed  gas  is  in¬ 
jected  at  the  bottom  of  a  quartz  tube  (inner  diameter  7.6  cm,  thickness  3  mm)  surrounded  by  a  coax¬ 
ial  5-turn  copper  induction  coil  (mean  diameter  8.6  cm)  traversed  by  a  radio  frequency  current.  The 
outer  Teflon  body  acts  as  an  electrical  insulator  and  electromagnetic  screen.  The  coil  is  cooled  with 
de-ionized  water  to  prevent  arcing  between  its  turns.  The  rf  current  produces  an  oscillating  axial 
magnetic  field  that  forces  the  free-electrons  to  spin  in  a  radial  plane  and  thereby  generates  eddy  cur¬ 
rents.  The  oscillating  magnetic  field  induces  an  oscillating  electric  field  and  the  induced  electro¬ 
magnetic  force  accelerates  the  free  electrons,  which  then  transfer  their  kinetic  energy  to  heavy  par¬ 
ticles  through  collisions. 
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Figure  1.  Schematic  cross-section  of  torch  head  with  7  cm  diameter  nozzle. 


The  plasma  torch  can  operate  with  a  variety  of  gases  (air,  hydrogen,  nitrogen,  methane, 
argon,  or  mixtures  thereof).  For  the  baseline  experiments  reported  here,  the  feed  gas  was  primar¬ 
ily  air  with  a  small  amount  of  hydrogen  (less  than  2%  mole  fraction)  for  purposes  of  electron 
number  density  measurements.  The  feed  gas  can  be  injected  in  axial,  radial  or  swirl  modes 
through  the  manifold  located  at  the  bottom  of  the  torch.  Axial  injection  provides  bulk  movement 
to  the  gas  during  the  start-up  phase.  In  normal  operation,  only  swirl  and  radial  injectors  are  used. 
As  will  be  seen  below,  the  swirl-to-radial  feed  ratio  has  a  large  impact  on  the  temperature  and 
concentration  profiles  of  the  plasmas  produced  by  the  torch. 

The  plasma  generated  in  the  coil  region  expands  into  ambient  air  through  a  converging 
copper  nozzle.  Different  nozzle  diameters  from  a  few  millimeters  to  7  cm  can  be  mounted  on  the 
torch;  for  the  present  work,  a  7  cm  nozzle  was  used.  At  the  nozzle  exit  plane  the  maximum  axial 
velocity  is  estimated  to  be  ~10  m/s,  the  maximum  temperature  ~7000  K,  the  density  p  s  5.04x 
10'2  kgm*3  and  the  dynamic  viscosity  p  =  1.6X10-4  kg.nr’s'1.  Based  on  the  nozzle  diameter  of 
7  cm,  the  Reynolds  number  at  the  nozzle  exit  is  about  220.  The  plasma  jet  is  therefore  laminar  at 
locations  of  1  and  5  cm  downstream  of  the  nozzle  exit  where  the  measurements  presented  below 
were  made.  Downstream  of  the  nozzle  exit  at  a  distance  equal  to  a  few  nozzle  diameters,  the 
plasma  plume  becomes  turbulent  as  a  result  of  mixing  with  ambient  air. 

The  temperature  and  electron  number  density  measurements  presented  in  the  following 
section  were  made  by  means  of  optical  emission  spectroscopy.  The  experimental  set-up,  shown 
in  Figure  2,  includes  a  SPEX  model  750M,  3/4  meter,  scanning  monochromator  fitted  with  a 
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1200  lines/mm  grating  blazed  at  500  nm  and  a  backthinned,  UV-coated  SPEX  Model  TE-2000 
Spectrum  One  thermoelectrically  cooled  CCD  camera.  The  CCD  chip  measures  30x12  mm  and 
contains  2000x800  square  pixels  of  dimension  15x15  |im.  Absolute  intensity  calibrations  were 
obtained  with  an  Optronics  model  OL550  radiance  standard  traceable  to  NIST  standards. 
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Figure  2.  Experimental  Set-up  for  Emission  Diagnostics. 


3.2.  Plasma  Characterization 

Figure  3  shows  photographs  of  the  plasma  plume  downstream  of  the  nozzle  exit  for  three 
different  conditions  of  swirl/radial  ratio.  In  the  “low  swirl  case”,  the  flow  rates  were  67  slpm 
(standard  liter  per  minute)  in  radial  mode  and  33  slpm  in  swirl  mode.  The  “medium”  and  “high” 
swirl  cases  correspond  to  radial/swirl  flow  rates  of  67/50  and  67/67,  respectively.  In  all  three 
cases,  the  plate  power  was  approximately  41.2  kW  and  a  small  quantity  of  hydrogen  (~2.3  slpm) 
was  premixed  prior  to  injection  into  the  torch.  To  a  good  approximation  the  flow  injected  into 
the  torch  is  thermodynamically  equivalent  to  humid  air  with  ~2.3  slpm  of  water  vapor. 


Figure  3.  Air  Plasma  Plume  for  three  Conditions  of  Radial/Swirl  Flow. 


As  can  be  seen  from  Figure  3,  the  swirl/radial  injection  ratio  has  a  noticeable  influence  on 
the  physical  aspect  of  the  plasma.  The  length  of  the  plume  is  approximately  35  cm,  20  cm,  and 
10  cm  for  the  low  (67/33),  medium  (67/50)  and  high  (67/67)  swirl  cases,  respectively.  The 
plasma  luminosity  exhibits  a  strong  radial  gradient  in  the  low  swirl  case,  but  in  contrast  it  is  al¬ 
most  radially  uniform  in  the  high  swirl  case  (it  is  not  possible  to  observe  radial  variations  of  the 
luminosity  in  Figure  3  as  the  photographs  are  intensity-saturated).  It  thus  appears  that  the  plasma 
properties  (temperature,  electron  number  density)  are  likely  to  be  more  uniform  radially  in  the 
case  with  highest  swirl  injection. 

Measurements  were  made  of  temperature  and  electron  number  density  radial  profiles  at 
locations  1  and  5  cm  downstream  of  the  nozzle  exit.  Temperature  profiles  were  determined  from 
the  absolute  intensity  of  the  atomic  line  of  oxygen  at  777.3  nm.  The  temperature  profiles  meas¬ 
ured  at  the  location  1  and  5  cm  downstream  of  the  nozzle  exit  for  a  plate  power  of  41.2  kW  are 
shown  in  Figure  4  and  Figure  5  for  both  the  low  and  high  swirl  cases.  As  can  be  seen  from  these 
figures,  radial  profiles  are  flatter  in  the  high  swirl  case  (67/67)  than  in  the  low  swirl  case  (67/33), 
in  accordance  with  the  visual  aspect  of  the  plasma  plume. 
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Figure  4.  1  cm  downstream  of  the  nozzle  exit.  Measured  temperature  profiles  from  Abel-inverted  absolute 
intensity  profiles  of  the  atomic  oxygen  triplet  at  777.3  nm.  Plate  power  =  41.2  kW.  Gas:  Air+2.3  slpm  H2. 


Figure  5.  5  cm  downstream  of  the  nozzle  exit.  Measured  temperature  profiles  from  Abel-inverted  absolute 
intensity  profiles  of  the  atomic  oxygen  triplet  at  777.3  nm.  Plate  power  =  41.2  kW.  Gas:  Air+2.3  slpm  H2. 

In  previous  studies  conducted  in  our  laboratory,  we  had  found  that  air  plasmas  generated 
by  the  torch  under  similar  conditions  of  temperature  and  velocity  were  close  to  Local  Thermody¬ 
namic  Equilibrium  (LTE)  [1],  This  observation  reflects  the  fact  that  characteristic  chemical  re¬ 
laxation  times  are  typically  10  times  longer  than  the  characteristic  flow  time  for  a  7  cm  diameter 
nozzle  (smaller  diameter  nozzles  produce  faster  flowing  plasmas  for  which  characteristic  flow 
times  may  become  smaller  than  chemical  times).  It  is  thus  expected  that  the  relatively  slow 
plasmas  flowing  through  the  7  cm  diameter  nozzle  are  close  to  LTE  both  at  1  and  5  cm  down¬ 
stream  of  the  nozzle  exit.  Under  LTE  conditions,  electron  number  densities  can  be  determined 
from  the  knowledge  of  the  plasma  temperature  using  chemical  equilibrium  relations  (Saha  equa¬ 
tion).  Figure  6  and  Figure  7  show  the  equilibrium  electron  number  density  profiles  correspond¬ 
ing  to  the  temperature  profiles  of  Figure  4  and  Figure  5.  In  order  to  verify  the  LTE,  we  also  made 
direct  electron  number  density  measurements  from  the  Stark-broadened  atomic  hydrogen  Balmer 
(3  line  at  486  nm. 
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Figure  6.  1  cm  downstream  of  the  nozzle  exit.  Measured  electron  number  density  profiles  from  Abel-inverted 
Hp  lineshapes  and  equilibrium  electron  number  density  profiles  based  on  the  temperature  profiles  of  Figure  4. 

Plate  power  =  41.2  kW.  Gas:  Air+2.3  slpm  H2. 


Figure  7.  5  cm  downstream  of  the  nozzle  exit.  Measured  electron  number  density  profiles  from  Abel-inverted 
Hp  lineshapes  and  equilibrium  electron  number  density  profiles  based  on  the  temperature  profiles  of  Figure  5. 

Plate  power  =  41.2  kW.  Gas:  Air+2.3  slpm  H2. 


In  air  plasmas,  the  Hp  line  sits  on  top  of  an  intense  emission  background  that  is  mainly 
comprised  of  bands  of  molecular  nitrogen  (second  positive  system).  In  order  to  extract  the  Hp 
lineshape,  spectral  measurements  were  made  both  with  the  mixture  of  air/hydrogen  (spectrum 
labeled  “Hp  +  background”  in  Figure  8)  and  with  pure  air  (spectrum  labeled  “Background”). 
Spectra  measured  at  several  lateral  locations  along  chords  of  the  plasma  were  then  Abel-inverted 
to  provide  local  emission  spectra  as  a  function  of  the  radial  location.  At  each  radial  location,  the 
Hp  lineshapes  were  recovered  by  subtracting  the  background  from  the  total  signal.  The  Hp  line- 
shapes  were  then  fitted  with  Voigt  profiles,  which  represent  the  convolution  of  several  broaden¬ 
ing  mechanisms  including  pressure  (van  der  Waals,  resonance),  Doppler,  instrumental,  and  Stark 
broadening  (see  Figure  9).  Pressure  and  Doppler  broadening  widths  only  depend  on  the  pressure 
and  temperature  of  the  gas.  Instrumental  broadening  was  minimized  by  using  a  very  small  en- 
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trance  slit  on  the  monochromator  (30  pm).  Radial  electron  number  density  profiles  were  deter¬ 
mined  with  the  aid  of  the  curves  of  Figure  9.  Results  are  shown  in  Figure  6  and  Figure  7. 


Wavelength  [A] 

Figure  8.  Hp  lineshape  extraction  procedure.  1  cm  downstream  of  the  nozzle  exit.  Low  swirl  case  (67/33  Ra¬ 
dial/Swirl).  Plate  power  =  41.2  kW.  Gas:  Air+2.3  slpm  H2.  The  two  spectra  in  the  figure  are  those  obtained 
after  Abel-inversion  at  r  =  10  mm  from  the  plasma  centerline. 


Figure  9.  Hp  line  broadening  in  atmospheric  pressure,  equilibrium  air.  Instrumental  broadening  is  well  ap¬ 
proximated  by  a  Gaussian  of  half  width  at  half  maximum  of  0.14  A.  1  cm  downstream  of  the  nozzle  exit.  Low 
swirl  case  (67/33  Radial/Swirl).  Plate  power  =  41.2  kW.  Gas:  Air+2.3  slpm  H2.  The  two  spectra  in  the  figure 
are  those  obtained  by  the  Abel-inversion  procedure  at  radial  location  r  =  10  mm. 


For  the  high  swirl  case  shown  in  Figure  7,  the  Hp  line  intensity  was  so  weak  relative  to  the 
nitrogen  background  (see  Figure  10)  that  it  was  not  possible  to  obtain  a  reliable  series  of  Abel- 
inverted  Hp  lineshapes.  For  this  case,  the  lines  were  not  Abel-inverted.  Nevertheless,  since  the 
plasma  temperature  does  not  vary  significantly  over  the  central  part  of  the  plasma,  the  electron 
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number  density  determined  from  the  Hp  lineshape  measured  along  the  diameter  of  the  plasma 
provides  an  estimate  of  the  average  electron  number  density  in  the  central  region.  As  can  be  seen 
from  Figure  7,  the  measured  line-of-sight  averaged  electron  density  agrees  well  with  the  expected 
equilibrium  value  in  the  central  region  of  the  plasma. 


9695  9700  9705  9710  9715 

Wavelength  [A] 

Figure  10.  Hp  lineshape  extraction  procedure.  Line-of-sight  Hp  lineshape  at  location  5  cm  downstream  of  the 
nozzle  exit.  High  swirl  case  (67/67  Radial/Swirl).  Plate  power  =  41.2  kW.  Gas:  Air+2.3  slpm  H2.  Here  the  Hp 
lineshape  was  measured  in  second  order  so  as  to  reduce  instrumental  broadening  to  approximately  0.07  A. 

The  foregoing  measurements  show  that  the  RF  plasma  torch  can  generate  steady-state  air 
plasmas  with  electron  number  densities  greater  than  1013  cm'3  over  volumes  with  dimensions 
greater  than  5  cm  in  all  directions.  The  measurements  also  indicate  that  the  shape  of  electron 
number  density  profiles  can  be  controlled  by  modifying  the  ratio  of  radial-to-swirl  injection.  The 
measurements  presented  here  were  obtained  with  a  mixture  of  air  and  hydrogen.  Measurements 
were  also  made  for  dry  air,  in  which  case  the  electron  number  density  could  only  be  determined 
by  assuming  chemical  equilibrium  at  the  local  temperature  measured  from  the  oxygen  triplet  at 
777.3  nm.  Results  of  this  series  of  experiments  are  shown  in  Figure  11  and  Figure  12  where  it 
can  be  seen  that  similar  conclusions  can  be  drawn  as  for  the  case  of  humid  air. 


17 


Figure  11.  Dry  air.  Temperature  profiles  at  location  5  cm  downstream  of  the  nozzle  exit  as  a  function  of  ra- 

dial/swirl  injection  ratio.  Plate  power  =  40.1  kW. 


Figure  12.  Dry  air.  LTE  electron  number  density  profiles  at  location  5  cm  downstream  of  the  nozzle  exit  as  a 
function  of  radial/swirl  injection  ratio.  Plate  power  =  40.1  kW. 


3.3.  Power  Budget 

The  main  objective  of  the  Air  Plasma  Rampart  program  is  to  significantly  reduce  the 
power  budget  required  to  generate  large  volume  air  plasmas.  Budget  requirements  can  be  de¬ 
fined  in  two  ways,  either  as  the  total  wall  plug  power  (which  depends  on  the  efficiency  of  the 
specific  device  utilized  to  generate  the  plasma)  or  as  the  net  power  deposited  into  the  plasma.  In 
this  work,  both  measurements  were  made  and  are  described  next. 

The  total  wall  power  was  determined  by  measuring  the  current  in  each  power  lead  of  the 
440  V,  triphase  power  supply,  by  means  of  a  Fluke  Model  33  ampmeter.  The  average  measured 
rms  current  was  ~72  A  (70,  72,  and  74  A  in  each  phase,  respectively).  The  total  wall  power  is 
then  given  by  the  following  expression: 
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Kail  =^3  Imsvms  cos<|>  s  48  kW 

where  cos<|),  the  power  factor,  is  equal  to  ~0.9  for  this  type  of  power  supply  according  to  the  manu¬ 
facturer’s  specifications,  and  where  Vrms  =  440  V.  As  the  plasma  probed  in  the  present  experiments 
was  approximately  cylindrical  of  diameter  6  cm  and  length  5  cm  (~140  cm3),  the  total  volumetric 
power  is  ~340  MW/m3. 


Additional  measurements  were  made  to  determine  the  torch  efficiency  and  the  power  depos¬ 
ited  into  the  plasma.  All  power  losses  can  be  determined  by  calorimetry  measurements  on  the  cool¬ 
ing  water  circuit.  To  this  end,  the  cooling  circuit  was  instrumented  with  a  combination  of  thermo¬ 
couples  at  the  inlet  and  outlet  of  the  generator  and  turbine  flowmeters  in  the  flow  lines.  The  total 
power  removed  by  the  cooling  water  was  measured  to  be: 

P cooling  water  ~  1®  Cp  AT  =  33  kW 


The  power  deposited  into  the  plasma  is  given  by: 

P plasma  ~  P wall  ~  P cooling  water  =  15  kW 

The  torch  efficiency,  defined  as  77  =  PpiaSma  /  Pwaii  is  approximately  equal  to  3 1  %.  Thus,  the  mini- 

'I 

mum  power  required  to  sustain  the  thermal  plasma  volume  is  ~105  MW/m  . 


The  total  volume  of  plasma  generated  in  the  low  swirl  case  (67/33)  depicted  in  Figure  3  is 
actually  much  larger  than  the  probed  volume  of  ~140  cm3.  We  estimate  that  the  total  volume  with 
electron  number  densities  greater  than  1013  cm'3  was  ~1000  cm3.  This  estimate  includes  the  volume 
of  plasma  generated  inside  the  torch  head  and  the  volume  extending  10  cm  downstream  of  the  noz¬ 
zle  exit.  The  total  volumetric  wall  power  is  then  ~48  MW/m3,  and  the  minimum  volumetric  power 
for  an  ideal  generator  is  ~15  MW/m3. 

3.4.  Conclusions  on  Baseline  Experiments 

Steady-state  air  plasmas  with  electron  number  densities  greater  than  1013  cm'3  and  volumes 
with  dimensions  greater  than  5  cm  in  all  directions  were  generated  in  both  dry  and  humid  air.  The 
shape  of  the  electron  number  density  profile  can  be  controlled  by  modifying  the  ratio  of  radial-to- 
swirl  injection.  The  total  power  required  to  sustain  open  volumes  of  thermal  air  plasma  is  measured 
to  be  about  340  MW/m3,  or  105/ri  MW/m3,  where  r|  represents  the  efficiency  of  the  specific  device 
used  to  produce  the  plasma.  It  is  important  to  emphasize  again  that  the  plasmas  generated  in  this 
work  are  thermal  plasmas,  i.e.  in  a  state  of  LTE,  and  accordingly  that  the  plasma  temperature  was 
always  greater  than  -4500  K.  Efforts  to  reduce  the  plasma  temperature  require  the  use  of  nonequi¬ 
librium  plasmas.  This  motivates  the  work  presented  in  the  rest  of  this  report. 
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4.  DC  glow  discharges  in  atmospheric  pressure  air 


4.1.  Introduction 

This  section  presents  experimental  and  numerical  investigations  to  determine  whether  and  to 
what  extent  the  electron  number  density  can  be  increased  in  air  plasmas  by  means  of  DC  dis¬ 
charges.  The  strategy  is  to  elevate  the  electron  temperature,  Te,  relative  to  the  gas  temperature,  Tg, 
with  an  applied  DC  electric  field. 

Section  4.2  describes  our  numerical  investigations  of  two-temperature  air  plasma  chemical 
kinetics.  We  present  first  the  two-temperature  kinetic  mechanism  that  we  have  developed  to  predict 
electron  number  density  in  air  at  a  given  gas  temperature,  as  a  function  of  the  electron  temperature. 
We  have  paid  close  attention  to  the  influence  of  the  electron  temperature  on  the  rate  coefficients, 
because  collisions  with  energetic  electrons  can  affect  the  vibrational  population  distribution  of 
molecules,  thereby  the  rates  of  ionization  and  dissociation. 

Section  4.3  discusses  the  implications  of  this  analysis  for  the  generation  of  nonequilibrium 
air  plasmas  by  means  of  electrical  discharges.  We  first  determine  in  Section  4.3.1  the  relation  be¬ 
tween  electron  number  density  and  current  density,  and  between  electron  temperature  and  electric 
field.  This  is  accomplished  with  Ohm’s  law  and  the  electron  energy  equation,  as  discussed  in  Sec¬ 
tion  4.3.2.  A  key  quantity  in  the  electron  energy  equation  is  the  rate  of  electron  energy  lost  by  ine¬ 
lastic  collisions.  To  predict  inelastic  losses  in  air  plasmas,  we  have  developed  a  detailed  colli- 
sional-radiative  model.  This  model  is  presented  in  Section  4.3.3. 

Section  4.4  describes  our  experiments  with  DC  glow  discharges  in  air.  We  demonstrate  that 
stable  diffuse  glow  discharges  with  electron  densities  of  up  to  ~2xl012  cm'3  can  be  sustained  in 
flowing  preheated  atmospheric  pressure  air.  The  electrical  characteristics  and  thermodynamic  pa¬ 
rameters  of  the  glow  discharges  are  measured  and  compared  determined. 

Section  4.5  compares  the  measured  electrical  characteristics  of  DC  glow  discharges  in  air 
with  those  obtained  with  the  two-temperature  and  collisional-radiative  model.  This  comparison 
validates  the  two-temperature  model  theoretical  predictions.  In  addition,  it  enables  us  to  establish 
the  power  requirements  of  DC  discharges  in  air  plasmas.  This  fundamental  understanding  forms 
the  basis  for  the  power  budget  reduction  strategy  presented  in  Section  5  of  this  report. 
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4.2.  Two-temperature  kinetic  simulations 

This  section  presents  results  of  numerical  investigations  to  determine  whether  and  to  what 
extent  electron  number  densities  can  be  increased  in  air  plasmas  by  elevating  the  electron  tempera¬ 
ture,  Te,  relative  to  the  gas  temperature,  Tg.  The  two-temperature  kinetic  mechanism  and  rates  used 
for  this  work  are  presented  in  Section  4.2.1.  In  Section  4.2.2,  the  kinetic  model  is  used  to  predict 
the  temporal  evolution  and  steady-state  species  concentrations  in  an  atmospheric  pressure  air 
plasma  with  constant  gas  temperature  of  2000  K  and  with  electron  temperatures  varied  from 
4000  K  to  18,000  K.  In  Section  4.2.3,  the  key  reactions  controlling  ionization  and  recombination 
processes  are  identified.  An  analytical  model  based  on  the  set  of  controlling  reactions  is  then  used 
to  predict  steady-state  species  concentrations  in  two-temperature  air.  As  will  be  seen,  the  analytical 
model  not  only  reproduces  the  CHEMKIN  solution  but  also  predicts  an  additional  range  of  steady- 
state  electron  number  densities. 


4.2.1.  Two-temperature  kinetic  model 


The  rate  coefficients  required  for  the  two-temperature  kinetic  model  depend  on  the  relative 
velocities  of  collision  partners  (related  to  Tg  for  reactions  between  heavy  particles  and  to  Te  for 
electron-impact  reactions)  and  on  the  population  distributions  over  internal  energy  levels  of  atoms 
and  molecules.  Thus,  these  rate  coefficients  correspond  to  the  weighted  average  of  elementary 
rates  over  internal  energy  states  of  atoms  and  molecules.  This  forms  the  basis  of  the  Weighted  Rate 
Coefficient  (WRC)  method  described  in  Refs  [2-5].  With  this  method,  for  instance,  the  total  rate 
coefficient  of  ionization  of  atomic  nitrogen  is  calculated  as: 


K7  =7^-Z^exP 
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and  the  rate  coefficient  of  electron-impact  dissociation  of  N2  as: 
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where  Qn  and  Qm  are  the  total  partition  functions  of  N  and  N2,  Ey  and  gy  the  energy  and  degener¬ 
acy  of  electronic  state  Y  of  N  or  N2,  (7 j(v)  and  Fyv(J)  the  vibrational  and  rotational  energies  of  level 
v  and  J,  and  K'yon  and  Ky*vJ  the  elementary  ionization  and  dissociation  rate  coefficients,  respec¬ 
tively. 


As  can  be  seen  from  Eqns  (1)  and  (2),  we  assume  that  the  internal  energy  levels  of  atoms 
and  molecules  are  populated  according  to  Boltzmann  distributions  at  the  electronic  temperature  Tei, 
the  vibrational  temperature  Tv,  and  the  rotational  temperature  Tr.  Elementary  rate  coefficients  are 
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calculated  from  cross-section  data  assuming  Maxwellian  velocity  distribution  functions  for  elec¬ 
trons  and  heavy-particles  at  Te  and  Tg,  respectively.  We  further  assume  that  Tei  =  Te  and  Tr  =  Tg. 
The  remaining  parameter,  Tv,  can  only  be  determined  in  the  general  case  by  solution  of  the  master 
equation  for  all  vibrational  levels  by  means  of  a  collisional-radiative  (CR)  model  that  incorporates 
vibrationally-specific  state-to-state  kinetics.  We  have  recently  developed  such  a  model  for  nitrogen 
plasmas  [2,  5]  that  provides  insight  into  the  relation  between  Tv  and  Tg  and  Te  in  atmospheric  pres¬ 
sure  plasmas.  The  nitrogen  CR  model  accounts  for  electron  and  heavy-particle  impact  ionization 
(atoms  and  molecules)  and  dissociation  (molecules),  electron-impact  vibrational  excitation,  V-T 
and  V-V  transfer,  radiation,  and  predissociation.  Through  comparisons  between  the  results  of  the 
CR  model  and  of  a  two-temperature  kinetic  model  of  nitrogen  that  assumed  either  Tv  =  Tg  or  Tv  = 
Te,  we  have  shown  [5,  6]  for  the  case  of  a  nitrogen  plasma  at  Tg  =  2000  K  that  the  steady-state  spe¬ 
cies  concentrations  determined  with  the  two-temperature  kinetic  model  are  in  close  agreement  with 
the  CR  model  predictions  if  one  assumes  (1)  that  Tv  =  Tg  for  electron  temperatures  Te  <  9500  K  and 
electron  number  densities  ne  <  ~10u  cm"3,  and  (2)  that  Tv  =  Te  or  Tv  =  Tg  for  Te  >  9500  K  and  ne  > 
-1015  cm'3  (in  the  latter  range,  best  agreement  is  obtained  with  Tv  =  Te  but  assuming  Tv  =  Tg  leads 
to  electron  number  densities  that  are  underestimated  by  at  worse  a  factor  of  5).  It  should  be  noted 
that  the  often-used  assumption  Tv  =  Te  produces  steady-state  electron  number  densities  that  are  sev¬ 
eral  orders  of  magnitude  greater  than  those  obtained  with  the  CR  model  for  electron  temperatures 
Te  <  9500  K  and  electron  number  densities  ne  <  ~10n  cm'3.  We  extend  these  results  to  atmospheric 
pressure  air  by  calculating  all  WRC  rate  coefficients  with  the  assumption  Tv  =  Tg. 

The  full  11 -species  (O2,  N2,  NO,  O,  N,  02+,  N2+,  NO+,  N+,  0+,  and  electrons),  40-reaction 
mechanism  and  rate  coefficients  for  the  case  Tg  =  2000  K  are  summarized  in  Table  1.  Electron  at¬ 
tachment  reactions  can  be  neglected  in  atmospheric  pressure  air  at  temperature  >1500  K  because 
the  equilibrium  concentrations  of  OC  or  O-  are  negligibly  small  relative  to  the  concentration  of 
electrons  above  -1500  K  (Figure  13).  For  reactions  between  nitrogen  species,  the  rate  coefficients 
are  taken  from  Yu  et  al.  [6].  This  set  is  supplemented  by  two-temperature  rate  coefficients  deter¬ 
mined  using  the  WRC  method  for  electron  impact  dissociation  and  ionization  of  O2  and  NO.  For 
electron-impact  ionization  of  O,  we  adopt  the  two-temperature  rate  of  Lieberman  and  Lichtenberg 

[7] ,  Rate  coefficients  for  O'  +  N2  <=>  NO+  +  N  and  0+  +  O2  <=>  02+  +  O  are  taken  from  Hierl  et  al. 

[8] ,  and  the  rate  coefficient  of  the  charge  transfer  reaction  between  0+  and  NO  is  calculated  using 
the  experimental  cross-section  reported  by  Dotan  and  Viggiano  [9].  The  remaining  reactions  in¬ 
volve  collisions  between  heavy  particles  and  thus  mostly  depend  on  the  gas  kinetic  temperature  (as 
we  assume  Tr  =  Tv  =  Tg).  For  these  reactions,  the  rate  coefficients  of  Park  [10,  1 1]  are  employed. 
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The  two-temperature  kinetic  calculations  presented  in  the  rest  of  this  section  were  made 
with  the  CHEMKIN  solver  [12]  modified  [13]  so  as  to  allow  a  different  temperature  (Te)  to  be 
specified  for  the  rates  of  particular  reactions.  The  extended  code  functions  in  similar  manner  as 
CHEMKIN.  For  thermal  reactions,  reverse  rates  are  computed  from  equilibrium  thermodynamic 
functions  (detailed  balance).  For  reaction  rates  that  depend  on  Te,  a  flag  is  entered  in  the 
CHEMKIN  reaction  mechanism  file  after  the  particular  reaction  considered  and  the  electron  tem¬ 
perature  is  specified.  When  reverse  rates  also  depend  on  Te,  the  reverse  rate  coefficient  must  be 
entered  explicitly.  Otherwise  the  reverse  rate  is  always  calculated  at  the  gas  temperature.  Reverse 
rates  with  a  dependence  on  Te  were  determined  with  the  WRC  model  by  detailed  balance. 


Table  1.  Two-temperature  kinetic  model  of  air  plasmas.  The  temperature  entering  the  Arrhenius-type 
expressions  is  either  the  gas  (Tg)  or  the  electron  (Te)  temperature,  as  indicated  in  columns  kf  (forward  rate)  and 
kr  (reverse  rate).  The  present  mechanism  if  valid  at  gas  temperatures  greater  than  1500  K. 


Reaction 

Temperature  dependence 
Kf  kr 

Rate  coefficient,  k  =  A  Tb  exp(-E/RT) 

A  (mole,cm,s)  b  E/R  (K) 

Ref. 

02  Dissociation/Recombination 

1.  02  +  02  =  20  +  02 

Tg 

Tg 

2.00E+21 

-1.5 

59500 

a 

2.  02  +  N0  =  0  +  0  +  N0 

Tg 

Tg 

2.00E+21 

-1.5 

59500 

a 

3.  02  +  N2  =  0  +  0  +  N2 

Tg 

Tg 

2.00E+21 

-1.5 

59500 

a 

4.  02  +  0  =  0  +  0  +  0 

Tg 

Tg 

1.00E+22 

-1.5 

59500 

a 

5.  02  +  N  =  0  +  0  +  N 

Tg 

Tg 

1.00E+22 

-1.5 

59500 

a 

6f.  02  +  E=>0  +  0  +  E 

Te 

2.85E+17 

-0.6 

59500 

b 

6b.  O  +  O  +  E  — >  O2  +  E 

Te 

4.03E+18 

-0.4 

0 

b 

NO  Dissociation/Recombination 

7.  no  +  o2=n  +  o  +  o2 

Tg 

Tg 

5.00E+15 

0.0 

75500 

a 

8.  NO  +  NO  =  N  +  0  +  NO 

Tg 

Tg 

1.10E+17 

0.0 

75500 

a 

9.  N0  +  N2  =  N+0  +  N2 

Tg 

Tg 

5.00E+15 

0.0 

75500 

a 

10.  NO  +  0  =  N  +  0  +  0 

Tg 

Tg 

1.10E+17 

0.0 

75500 

a 

11.  NO  +  N  =  N  +  0  +  N 

Tg 

Tg 

1.10E+17 

0.0 

75500 

a 

12f.  NO  +  E=>N  +  0  +  E 

Te 

3.54E+16 

-0.2 

75500 

b 

12b.  N  +  0  +  E  =>  NO  +  E 

Te 

8.42E+21 

-1.1 

0 

b 

N2  Dissociation/Recombination 

13.  N2  +  02  =  N  +  N  +  02 

Tg 

Tg 

7.00E+21 

-1.6 

113200 

a 

14.  N2  +  NO-N  +  N  +  NO 

Tg 

Tg 

7.00E+21 

-1.6 

113200 

a 

15.  N2  +  N2=N  +  N  +  N2 

Tg 

Tg 

7.00E+21 

-1.6 

113200 

a 

16.  N2  +  0  =  N  +  N  +  0 

Tg 

Tg 

3.00E+22 

-1.6 

113200 

a 

17.  N2  +  N  =  N  +  N  +  N 

Tg 

Tg 

3.00E+22 

-1.6 

113200 

a 

18f.  N2  +  E=>N  +  N  +  E 

Te 

1.18E+18 

-0.7 

113200 

b 

18b.  N  +  N  +  E=>N2  +  E 

Te 

1.36E+23 

-1.3 

0 

b 

Zeldovich  reactions 

19.  N2  +  0  =  N0  +  N 

Tg 

Tg 

6.40E+17 

-1.0 

38400 

a 

20.  N0  +  0  =  02  +  N 

Tg 

Tg 

8.40E+12 

0.0 

moo 

a 

Associative  ionization/Dissociative  recombination 

2  If.  N  +  0  ->  NO+  +  E  Tg 

8.80E+08 

1.0 

31900 

a 

21b.  NO+  +  E  =>  N  +  0 

Te 

9.00E+18 

-0.7 

0 

c 

22  f.  N  +  N=>N2+  +  E 

Tg 

6.00E+07 

1.5 

67500 

b 

22b.  N2+  +  E=>N  +  N 

Te 

1.53E+18 

-0.5 

0 

b 

23  f.  0  +  0=>02+  +  E 

Tg 

7.10E+02 

2.7 

80600 

a 

23b.  02+  +  E=>0  +  0 

Te 

1.50E+18 

-0.5 

0 

c 

Electron  impact  ionization/Three-body  recombination 

24f.  0  +  E=>0+  +  E  +  E  Te 

7.74E+12 

0.7 

157760 

d 

24b.  0+  +  E  +  E=>0  +  E 

Te 

2.19E+21 

-0.8 

0 

e 

25f.  N  +  E=>N+  +  E  +  E 

Te 

5.06E+19 

0.0 

168200 

b 

25b.  N+  +  E  +  E=>N  +  E 

Te 

5.75E+26 

-1.3 

0 

b 

26f.  02  +  E  =>  02+  +  E  +  E 

Te 

5.03E+12 

0.5 

146160 

b 
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26b. 

02+  +  E  +  E  =>  02  +  E 

Te 

8.49E+23 

-1.9 

0 

b 

27f. 

N2  +  E  =>  N2+  +  E  +  E 

Te 

2.70E+17 

-0.3 

181000 

b 

27b. 

n2+  +  e  +  e=>n2  +  e 

Te 

2.05E+21 

-0.8 

0 

b 

28f. 

NO  +  E  ->  NO+  +  E  +  E 

Te 

2.20E+16 

-0.3 

107400 

b 

28b. 

NO+  +  E  +  E  =>  NO  +  E 

Te 

2.06E+25 

-2.0 

0 

b 

Charge  exchange/Charge  transfer 

29  f. 

n+  +  n2=>n2+  +  n 

Tg 

4.60E+11 

0.5 

12200 

b 

29b. 

N2+  +  N=>N2  +  Nt 

Tg=2000K 

1.93E+13 

0.0 

0 

b 

30. 

NO+  +  O  =  N+  +  02 

Tg 

Tg 

1.00E+12 

0.5 

77200 

a 

31. 

NO  +  0+  =  N+  +  02 

Tg 

Tg 

1.40E+05 

1.9 

15300 

a 

32. 

0+  +  N2  =  NO+  +  N 

Tg 

Tg 

4.40E+13 

0.0 

5664 

f 

33. 

o++n2=n2+  +  o 

Tg 

Tg 

9.00E+1 1 

0.4 

22800 

a 

34. 

no++n=n2+  +  o 

Tg 

Tg 

7.20E+13 

0.0 

35500 

a 

35. 

o2+  +  n=n+  +  o2 

Tg 

Tg 

8.70E+13 

0.1 

28600 

c 

36. 

o2++n2=n2+  +  o2 

Tg 

Tg 

9.90E+12 

0.0 

40700 

c 

37. 

no+  +  o2  =  o2+  +  no 

Tg 

Tg 

2.40E+13 

0.4 

32600 

c 

38. 

no+  +  o  =  o2+  +  n 

Tg 

Tg 

7.20E+12 

0.3 

48600 

c 

39. 

0+  +  02  =  02+  +  0 

Tg 

Tg 

3.26E+13 

0.0 

2064 

f 

40. 

0+  +  N0=N0+  +  0 

Tg 

Tg 

2.42E+13 

0.0 

902 

g 

a.  Park  [1 1] 

b.  WRC  [4,  5],  These  rates  were  calculated  at  Tg=Tv=2000  K.  The  present  fitting  formulas  are  valid  for  6000  <  Tc  <  20,000  K 

c.  Park  [10] 

d.  Lieberman  [7] 

e.  Detailed  Balance 

f.  Hierletal.  [8] 

g.  Dotan  and  Viggiano  [9] 


Temperature  (K) 

Figure  13.  Ion  concentrations  relative  to  the  electron  concentration  in  equilibrium  air  (P  =  1  atm). 


4.2.2.  Results 

We  consider  first  the  case  of  an  air  plasma  taken  to  be  in  equilibrium  (Tg°  =  Te°  =  2000  K,  P 
=  1  atm)  at  time  zero  when  an  elevated  electron  temperature  is  instantaneously  prescribed,  in  an 
idealized  way  modeling  an  electrical  glow  discharge  in  a  reactor  section.  In  the  example  shown  in 
Figure  14,  the  gas  temperature  is  held  constant  at  2000  K  and  the  electron  temperature  is  increased 
to  13,000  K  at  time  zero.  The  time  evolution  of  species  concentrations  computed  with  the  two- 
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shown  in  Figure  14,  the  gas  temperature  is  held  constant  at  2000  K  and  the  electron  temperature 
is  increased  to  13,000  K  at  time  zero.  The  time  evolution  of  species  concentrations  computed 
with  the  two-temperature  CHEMKIN  solver  is  shown  in  Figure  14.  The  electron  number  density 
rises  from  its  initially  low  value  of  3.3x1 06  cm"3  to  a  steady-state  value  of  ~4xl012  cm3  in  about 
0.1  ms.  The  dissociation  fraction  of  oxygen  atoms  increases  from  ~0.03%  at  time  zero  to  ~1%  at 
steady-state.  NO+  is  the  dominant  ion  at  all  times. 

Additional  calculations  were  made  for  various  electron  temperatures  while  keeping  the 
gas  temperature  constant  at  2000  K.  The  predicted  temporal  evolutions  of  the  electron  number 
density  are  shown  in  Figure  15.  Practically  no  increase  in  the  electron  number  density  is  ob¬ 
served  for  electron  temperatures  below  a  threshold  value  of  Te  s  6000  K,  which  corresponds  to 
the  temperature  where  electron-impact  ionization  reactions  begin  to  dominate  over  heavy  parti¬ 
cle-impact  dissociation.  As  the  electron  temperature  is  further  increased,  the  steady-state  electron 
concentration  increases  significantly,  with  a  very  abrupt  change  at  Te  =  16800  K. 


Time  (ms) 


Figure  14.  Temporal  evolution  of  species  concentrations  in  atmospheric  pressure  air  at  constant  gas  tempera¬ 
ture  (2000  K)  and  constant  electron  temperature  (13000  K). 
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Figure  15.  Temporal  evolution  of  the  electron  number  density  in  a  two-temperature  air  plasma.  Initial  condi¬ 
tions  are  equilibrium  air  at  2000  K  (ne<t=#>  =  3.3xl06  cm'3). 

Figure  16  shows  the  steady-state  electron  number  densities  predicted  with  CHEMKIN  as 
a  function  of  the  electron  temperature.  The  abrupt  change  in  electron  number  density  values  oc¬ 
curs  for  an  electron  temperature  of  — 1 6,800  K  where  the  predicted  steady-state  electron  number 
density  suddenly  increases  from  -l.OxlO14  to  -lxlO18  cm'3  over  a  few  Kelvin.  It  is  interesting  to 
examine  the  reverse  case  where  steady-state  electron  concentrations  are  calculated  for  an  initial 
composition  given  by  the  steady-state  solution  corresponding  to  Tg  =  2000  K  and  Te  =  20,000  K 
(corresponding  to  ne(t=0)  =  1.7xl018  cm'3).  As  can  be  seen  from  Figure  17,  in  this  case  the  pre¬ 
dicted  steady-state  electron  number  densities  start  by  decreasing  along  the  same  curve  as  in 
Figure  16,  but  instead  of  the  abrupt  decrease  at  16,800  K,  continue  their  slow  decrease  until  the 
electron  temperature  reaches  ~1 4,300  K.  When  the  electron  temperature  is  further  decreased  be¬ 
low  -14,300  K,  the  steady-state  electron  number  density  abruptly  decreases  to  the  level  of  the 
curve  of  Figure  16.  Thus  a  hysteresis  occurs  as  the  electron  temperature  is  increased  and  de¬ 
creased  in  a  cyclical  fashion. 
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8000  10000  12000  14000  16000  18000  20000  22000 

Electron  Temperature,  Te  (K) 

Figure  16.  Steady-state  electron  number  density  predicted  by  CHEMKIN  for  air  at  Tg  =  2000  K,  as  a  function 
ofTe.  For  each  steady-state  calculation,  initial  conditions  are  equilibrium  air  at  2000  K. 
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8000  10000  12000  14000  16000  18000  20000  22000 

Electron  Temperature,  Te  (K) 

Figure  17.  Steady-state  electron  number  density  predicted  by  CHEMKIN  for  air  at  Tg  =  2000  K,  as  a  function 
of  Te.  For  each  steady-state  calculation,  the  initial  condition  corresponds  to  the  steady-state  composition  pre- 
dieted  by  CHEMKIN  at  Tg  =  2000  K  and  Te  =  20,000  K  (ne(t=#)  =  1.7xl018  cm’3). 


4.2.3.  Analysis  of  the  ionization  mechanisms 


Through  detailed  examinations  of  the  reactions  and  rates,  we  found  that  the  behavior  in 
each  of  the  two  regions  A  and  B  of  the  curve  in  Figure  16  can  be  explained  in  terms  of  the  fol¬ 
lowing  simplified  reaction  mechanisms: 
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4.2.3 .a.  Ionization  mechanism  in  Region  A 

In  region  A,  the  initial  rapid  electron  concentration  rise  (see  Figure  14  for  the  case  Te  = 
13,000  K)  is  the  result  of  electron-impact  ionization  of  N2  and  O2  via  three-body  reactions: 

O2  +  e  02+  +  e  +  e 
N2  +  e  =>  N2+  +  e  +  e 

and  via  electron-impact  dissociation  of  O2  followed  by  electron-impact  ionization  of  O: 

02  +  e=i>0  +  0  +  e 
O  +  e  =>  0+  +  e  +  e 

The  charged  species  produced  by  these  processes  undergo  rapid  charge  transfer  to  NO+,  via 

o++n2^  no++n 

N2+  +  O2  =>  N2  +  02+ 

o2++no=>  no+  +  o2 

The  main  path  for  electron  recombination  is  the  two-body  dissociative  recombination  reaction: 

NO+  +  e  =>  N  +  O 

When  the  concentration  of  NO+  becomes  sufficiently  large,  the  rate  of  dissociative  recombination 
balances  the  rate  of  electron  production  and  the  plasma  reaches  steady-state.  Thus,  in  Region  A, 
the  termination  step  for  the  ionization  process  is  the  two-body  recombination  of  a  molecular  ion. 


4.2,3.b.  Ionization  mechanism  in  Region  B 

An  example  of  the  temporal  evolution  of  species  concentrations  in  Region  B  is  shown  in 
Figure  18  where  the  electron  temperature  is  fixed  at  Te=l  8,000  K.  As  in  region  A,  the  electron 
number  density  initially  increases  by  electron-impact  ionization  of  N2  and  O2  via  : 

O2  +  e  =>  02+  +  e  +  e 
N2  +  e  =>  N2+  +  e  +  e 

and  via  electron-impact  dissociation  of  O2  followed  by  electron-impact  ionization  of  O: 

02  +  e=>0  +  0  +  e 
0  +  e=>0+  +  e  +  e 

The  difference  with  Region  A  is  that  the  charge  transfer  reactions  are  not  fast  enough  to 
produce  NO+  at  a  high  enough  rate.  This  is  because  these  reactions  are  controlled  by  the  gas  tem¬ 
perature,  whereas  electron  impact  ionization  reactions  are  controlled  by  Te.  The  critical  electron 
temperature  that  defines  the  limit  between  Regions  A  and  B  corresponds  approximately  to  the  elec¬ 
tron  temperature  for  which  the  rate  of  the  transfer  reaction  Cf  +  N2  =>  NO+  +  N  is  comparable  with 
the  rate  of  avalanche  ionization  by  electron  impact.  Above  this  critical  electron  temperature,  the 
avalanche  ionization  process  continues  until  all  molecular  species  are  dissociated.  Eventually  the 
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Eventually  the  rates  of  three-body  electron  recombination  reactions  balance  the  rate  of  ionization, 
and  steady-state  is  reached. 

It  is  noted  that  in  Region  A,  electron  impact  dissociation  of  N2  (or  NO)  is  negligible  be¬ 
cause  the  dissociation  energy  of  N2  (9.76  eV)  is  much  larger  than  that  of  02  (5.11  eV),  and  the 
concentration  of  NO  is  small  relative  to  the  concentration  of  02.  It  is  only  above  the  critical  tem¬ 
perature  that  electron  impact  dissociation  of  N2  starts  having  a  noticeable  effect. 


Figure  18.  Temporal  evolution  of  species  concentrations  in  atmospheric  pressure  air  at  constant  gas  tempera¬ 
ture  (2000  K)  and  constant  electron  temperature  (18,000  K). 


4.2. 3 .c.  Analytical  solution 

The  kinetics  in  Regions  A  and  B  can  be  described  with  a  simplified  subset  of  reactions 
that  takes  into  account  the  dominant  channels  discussed  in  the  foregoing  section.  With  this  sim¬ 
plified  mechanism,  the  steady-state  concentrations  of  major  species  are  obtained  by  solving  the 
following  system  of  equations: 

•  Steady-state  for  e_: 

02  +  e  =>  02+  +  e  +  e 
N2  +  e  =>  N2+  +  e  +  e 
0  +  e  <=>  0+  +  e  +  e 
NO+  +  e=>N  +  0 

•  Steady-state  for  02: 

02  +  e=>0  +  0  +  e 

0  +  0  +  M=>02  +  M,  with  M  =  02,  N2,  O 

•  Steady-state  for  NO+: 
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N0+  +  e  =>  N  +  O 

o+  +  n2^no++n 
o2++no=>no++o2 

•  Steady-state  for  0+: 

O  +  eo  0+  +  e  +  e 
0+  +  N2z^N0++N 

0+  +  o2  =>  o2+  +  o 

•  Steady-state  for  02+: 

02  +  e  =>  02+  +  e  +  e 

n2+  +  o2  =>  o2+ + n2 
o2++no=>no++o2 

•  Steady-state  for  N2+: 

N2  +  e  =^>  N2+  +  e  +  e 
N2+  +  02  02+  +  N2 

In  writing  the  corresponding  steady-state  relations,  we  make  the  approximation  that  the 
change  of  plasma  volume  due  to  the  increase  of  the  total  number  of  moles  at  constant  T  and  P  is 
negligible,  and  that  the  concentration  of  N2  remains  constant  and  equal  to  its  initial  value.  Further¬ 
more,  we  consider  that  the  dominant  neutral  species  in  the  plasma  are  02,  O,  and  N2. 

By  elimination  of  no,  no2,  no+,  no2+,  nN2+  and  nNo+,  we  obtain  a  tenth-degree  polynomial  in 
ne  with  coefficients  that  only  depend  on  Tg,  Te,  no2(0),  nN2(0)  and  the  rate  coefficients.  The  roots  of 
this  polynomial  were  extracted  with  the  mathematical  package  MATLAB.  The  roots  that  give 
negative  values  for  no  are  omitted  from  the  solution.  The  remaining  roots  of  the  polynomial  are 
plotted  in  Figure  19  along  with  the  CHEMKIN  predictions  corresponding  to  the  full  mechanism  of 
Table  1.  As  can  be  seen  from  Figure  19,  the  approximate  solution  obtained  with  the  simplified 
mechanism  is  in  very  good  agreement  with  the  CHEMKIN  predictions  in  regions  A  and  B.  The 
small  discrepancy  in  Region  B  is  due  to  the  neglecting  in  our  simplified  model  of  nitrogen  dissocia¬ 
tion.  Furthermore,  the  tenth-degree  polynomial  exhibits  an  extra  solution  that  could  not  be  attained 
with  CHEMKIN  (Region  C).  The  limits  of  Region  C  are  the  turning  points  labeled  (a)  and  (p)  in 
Figure  19.  If  we  initialize  CHEMKIN  with  the  plasma  composition  corresponding  to  a  point  in  Re¬ 
gion  C  of  Figure  19,  CHEMKIN  produces  a  new  steady-state  electron  number  density  located  on 
either  the  lower  (Region  A)  or  upper  (Region  B)  limb  of  the  steady-state  curves.  Thus  Region  C  of 
Figure  19  cannot  be  obtained  by  fixing  the  electron  temperature. 
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Figure  19.  Steady-state  electron  number  densities  predicted  by  CHEMKIN  and  by  the  tenth-order  polynomial 
analytical  solution. 

Figure  20  shows  the  concentrations  of  dominant  species  as  predicted  by  the  analytical 
model.  We  see  that  the  concentration  of  NO+  increases  up  to  turning  point  (a),  and  then  stays 
approximately  constant  as  charge  transfer  becomes  slower  than  oxygen  ionization.  The  concen¬ 
tration  of  oxygen  atoms  steadily  increases  throughout  Regions  (A)  and  (C).  Beyond  turning  point 
(p),  molecular  oxygen  is  nearly  fully  dissociated.  The  electron  concentration  is  approximately 
equal  to  the  concentration  of  NO+  in  Region  (A)  and  to  the  concentration  of  0+  in  Regions  (B) 
and  (C). 


Figure  20.  Steady-state  species  concentrations  at  Tg  =  2000  K,  as  predicted  by  the  analytical  solution. 
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4.3.  Predicted  Electric  Discharge  Characteristics 

The  electron  number  density  and  electron  temperature  can  be  related  to  the  current  density 
and  electric  field,  respectively,  by  means  of  Ohm’s  law  and  the  electron  energy  equation.  These 
relations  provide  guidance  for  the  design  of  nonequilibrium  DC  discharges,  as  well  as  an  estimate 
of  the  power  requirements  of  such  discharges. 


4.3.1.  Ohm’s  law 


Ohm’s  law  relates  the  current  density  j  to  the  electric  field  E: 

]  =  <?E  =  — 4=, —  E 


m„ 


veh 


(3) 


where  a  is  the  electrical  conductivity  of  the  plasma,  me  the  electron  mass,  and  veh  the  average  fre¬ 
quency  of  collisions  between  electrons  and  heavy  particles.  The  total  collision  frequency  is  the  sum 
of  the  collision  frequencies  of  electrons  with  neutrals,  n ,  and  ions,  /: 

Kh=Kn+Ki  (4) 


The  electron-neutral  collision  frequency  ven  can  be  expressed  in  terms  of  the  number  densities  of 
neutral  species,  n„,  the  electron  velocity,  ge  =  ^8kTe  / 7ime ,  and  the  energy-averaged  momentum 
transfer  cross-sections  Qm  as: 

Kn  =  nn  ge  Qen 


The  energy-averaged  electron-neutral  momentum  transfer  collision  cross-sections  Qeh  are  calculated 
by: 


Qeh  =tif2exp[-7- 1 Qeh(e)de 


3  T: 


(5) 


v  w 


where  Te  is  the  electron  temperature,  s  is  the  electron  impact  energy,  and  Qeh(s)  is  the  momentum 
collision  cross-section.  For  N2,  O2  and  O,  these  cross-sections  have  been  taken  from  Brown  [14], 
Shkarofsky  et  al.  [15],  and  Tsang  et  al.  [16].  The  resulting  average  energy  cross-sections  Qeh  are 
presented  in  Figure  21. 
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Figure  21:  Energy-averaged  momentum  transfer  cross  sections  for  collisions  between  electrons  and  N2,  Oi,  O. 

In  Region  A,  where  the  dominant  neutral  species  are  N2  and  O2,  and  where  ions  have  negligible 
concentrations,  the  total  collision  frequency  is  well  approximated  by  the  electron-neutral  collision 
frequency. 

In  Region  B,  where  the  plasma  is  almost  fully  ionized,  the  total  collision  cross  section  is  ap¬ 
proximately  equal  to  the  electron-ion  collision  frequency,  which  can  be  expressed  as  [17]: 

Vei  =  3.64x10-6h,^^-  sec-1,  (6) 

where  A  represents  the  ratio  of  the  Debye  length  to  the  impact  parameter  for  90°  scattering  and  is 
approximately  equal  to  2.5  for  electron  number  densities  of  1018  cm'3  (Mitchner  and  Kruger 
[17]). 


4.3.2.  Electron  energy  equation 


4.3.2.a.  Introduction 

For  a  stationary  plasma  in  a  DC  electric  field,  the  electron  energy  equation  can  be  written: 

(7) 


OE1  =  n,lh (T.  -Tt )^v„  +  +  R 


mL 


In  Eqn.  (7),  k  is  the  Boltzmann  constant  and  Th  the  kinetic  temperatures  of  heavy  species  (as¬ 
sumed  to  be  the  same  for  all  heavy  species  and  thus  equal  to  Tg).  The  term  on  the  left  hand  side 
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second  terms  on  the  RHS  represent  the  volumetric  power  lost  by  free  electrons  through  elastic  and 
inelastic  collisions,  respectively,  and  the  last  term  on  the  RHS  stands  for  volumetric  radiative 
losses. 


In  Region  A  of  the  S-shaped  curve,  inelastic  energy  losses  dominate  the  elastic  and  radiation 
losses  by  at  least  two  orders  of  magnitude.  In  Region  B,  however,  inelastic  losses  are  negligible 
relative  to  elastic  and  radiation  losses.  In  the  rest  of  this  analysis,  we  limit  ourselves  to  the  lower 
limb  of  the  S-shaped  curve  (Region  A),  and  therefore  we  neglect  radiation  losses  in  the  rest  of  this 
work. 


The  basis  for  calculations  of  inelastic  losses  in  atmospheric  pressure  air  plasmas  is  summa¬ 
rized  below.  Electrons  lose  energy  through  the  following  inelastic  processes:  vibrational  excitation 
of  molecular  species  (VE  transfer),  ionization,  electronic  excitation,  dissociation  of  molecules,  elec¬ 
tronic  excitation  and  ionization  of  atoms.  At  electron  temperatures  below  about  17,000  K,  the  main 
channel  for  inelastic  electron  energy  loss  in  air  is  via  electron  impact  vibrational  excitation  of  nitro¬ 
gen.  This  process  is  particularly  important  for  the  ground  state  of  molecular  nitrogen  because  vi¬ 
brational  excitation  by  electron  impact  of  this  state  occurs  via  resonant  transitions  to  the  ground 
state  of  the  unstable  negative  ion  N2-: 

N2  (X,v")  +  e  ->  N2  (X,v)  N2  (X,v’)  +  e 

The  net  rate  of  energy  lost  by  this  process  is: 

8vV=(X,.4N2(X,v")][e]-Xv,v.[N2(X,v')][e])AEvV 

where  Xv’v  and  Kv\”  are  the  rate  coefficients  of  electron-impact  vibrational  excitation  and  de¬ 
excitation,  and  AEv’v  stands  for  the  difference  of  energy  between  the  two  vibrational  levels.  The 
net  rate  of  inelastic  losses  is  a  function  of  the  gas  and  electron  temperatures,  the  electron  number 
density,  and  the  vibrational  population  distribution  of  the  ground  state  of  N2.  In  the  limiting  case 
where  the  vibrational  level  populations  follow  a  Boltzmann  distribution  at  the  electron  temperature, 
the  net  rate  of  energy  loss  by  VE  transfer  is  equal  to  zero  because  the  energy  lost  by  VE  excitation 
reactions  is  exactly  balanced  by  the  energy  gained  from  de-excitation.  In  the  other  limiting  case 
where  the  vibrational  levels  follow  a  Boltzmann  distribution  at  the  gas  temperature,  the  rate  of  exci¬ 
tation  is  much  larger  than  the  rate  of  de-excitation.  In  the  general  case,  the  vibrational  population 
distribution  is  intermediate  between  the  two  previous  cases.  The  vibrational  population  distribution 
is  then  governed  by  the  relative  importance  of  the  rates  of  vibrational  excitation  by  electron  impact, 
and  the  rates  of  de-excitation  which  are  mostly  determined  by  collisions  with  heavy  species.  The 
dominant  de-excitation  processes  are  vibrational-vibrational  transfer  between  two  N2  molecules  (V- 
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V  transfer),  vibrational-vibrational  transfer  between  one  N2  molecule  and  other  molecular  species 
such  as  O2  and  NO  (V-V’  transfer),  and  vibrational-translation  (V-T)  relaxation  by  collisions  of  N2 
with  other  heavy  species  (N2,  O2,  NO,  N  and  O).  Here  the  main  relaxation  processes  are  V-T  re¬ 
laxation  by  O  and  N2.  To  calculate  inelastic  energy  losses  in  air,  we  must  therefore  predict  the  vi¬ 
brational  population  distribution  of  the  nitrogen  ground  state  by  taking  into  account  the  aforemen¬ 
tioned  processes.  This  is  a  complex  calculation  that  requires  the  use  of  a  vibrationally  specific  Col- 
lisional-Radiative  (CR)  model.  We  have  developed  such  a  model  for  pure  nitrogen  plasmas  [2,  3,  5, 
18]  and  have  recently  extended  it  to  air  plasmas  [19]. 


4.3.2.b.  Rate  coefficients  controlling  the  vibrational  distribution  of  N?  levels 

The  rate  coefficients  for  VE  transfer  are  calculated  [5]  using  the  cross  section  calculation 
method  of  Kazansky  and  Yelets  [20].  This  method  reproduces  available  experimental  cross-sections 
relative  to  the  low-lying  vibrational  levels  within  about  10%. 

For  electron  temperatures  up  to  17000K  (turning  point  a  in  the  S-shaped  curve),  the  domi¬ 
nant  N2  V-T,  V-V  and  V-V’  relaxation  processes  are: 

N2(Y,v  +  \)  +  M  <->  N2(Y,v)  +  M  (8) 

N2(X,vt+l)  +  AB(X,v2)^  N2 (X, v, )  +  AB(X, v2  + 1)  (9) 

In  Eqn.  (8),  M  represents  the  heavy  particle  collision  partner  M=N2,  O2 ,  NO  or  O  and  in  Eqn.  (9), 
AB  is  the  diatomic  molecule  TV),  O2  or  NO. 

To  our  knowledge,  no  experimental  data  have  been  reported  for  the  V-T  relaxation  of  N2  by 
collisions  with  N.  However  Kozlov  et  al.  [21]  experimentally  determined  an  upper  limit  value  of 
the  V-T  relaxation  rate  for  transition  v=l—>  v=0.  They  showed  that  for  temperatures  between 
2500K  and  4500K,  this  rate  is  about  one  order  of  magnitude  lower  than  the  rate  of  V-T  relaxation 
of  N2  by  collisions  with  O  atoms.  Since  for  our  conditions  the  concentration  of  N  atoms  is  at  least 
four  orders  of  magnitude  lower  than  the  concentration  of  O  atoms,  we  neglect  the  V-T  relaxation  of 
N2  by  collisions  with  N  atoms. 

Most  measured  vibrational  relaxation  rate  coefficients  are  for  transitions  between  vibra¬ 
tional  levels  v=0  and  1.  When  available,  experimental  rates  have  been  preferred  over  theoretical 
ones.  The  existing  experimental  rates  have  been  compared  and  critically  selected.  Rates  for  transi¬ 
tions  between  higher  levels  have  been  calculated  using  scaling  functions  derived  from  SSH  theory, 
which  is  a  reasonable  approximation  when  the  gas  temperature  is  below  3000K.  The  reverse  rates 
have  been  determined  by  application  of  the  detailed  balance  method. 
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For  each  V-T  transfer  process,  the  rates  ki,o  corresponding  to  transition  v=l—>v=0  have 
been  calculated  with  the  following  analytical  expression: 

f  n  r  Y  ,  _  ,  , 

--7JT+—  1-Dexp  --a  (10) 
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where  ki,o  is  expressed  in  cm3 s'1,  and  Tg  and  Ew  (energy  of  the  N2  transition  v=l->v=0)  are  ex¬ 
pressed  in  Kelvin.  The  parameters  A,  B,  C,  D,  m  and  n  are  listed  in  Table  2. 


V-T  relaxation  of  N2  by  collision  with  O  atoms  has  been  studied  experimentally  between 
300K  and  4500K  by  Breshears  and  Bird  [22],  McNeal  et  al.  [23]  and  Eckstrom  [24],  These  three 
sets  of  experimental  data  can  be  well  fitted  by  the  Landau-Teller  expression  of  Taylor  [25],  which 
reproduces  the  experimental  results  within  a  factor  less  than  two  for  temperatures  lower  than 
3000K.  The  parameters  of  Taylor’s  expression  are  reported  in  Table  2. 


V-T  relaxation  of  N2  by  collision  with  N2  molecules  has  been  investigated  by  Lukasik  and 
Young  [26],  Henderson  [27],  Millikan  and  White  [28],  Hurle  [29,  30],  Appleton  [31],  Kovacs  and 
Mack  [32]  and  Kozlov  et  al.  [21].  The  analytical  expression  proposed  by  Nikitin  [33]  agrees  within 
a  factor  less  than  two  with  the  experimental  values  for  temperatures  between  1500K  and  3000K. 
The  parameters  of  Nikitin’s  expression  are  reported  in  Table  2. 

No  experimental  data  have  been  found  for  V-T  relaxation  of  V?  by  collision  with  O2  and  NO 
molecules.  Since  these  collision  partners  have  about  the  same  mass  as  N2,  we  have  used  the  same 
parameters  A,  B,  C,  D,  m  and  n  to  describe  ki,o. 

Table  2:  parameters  for  the  V-T  rates  k]>0 

Mb  A  B  C  D  m  n 

O  1.07xl0'10  69.9  0  0  0  0 

N2,02,N0  7.8xlQ-12  218  690  1  1  1 


The  rate  coefficients  kv+  /,v  for  transitions  v  +  1  — >  v  between  upper  vibrational  levels  have 
been  calculated  using  appropriate  scaling  laws  from  the  measured  kw  rates: 

*v+1.v=V?(v  + 1)  (H) 

Using  SSH  theory  [34]  and  some  approximations  for  the  Morse  oscillator  model  [35],  G(v+1)  can 
be  expressed  as: 
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G  (v+l)D 


(v  +  lXl-xJ^^,.,) 


l-xe(v+l)  F  (^10 ) 
where  xe  is  the  anharmonicity  of  the  N2  molecule,  andyv+  i,v  is  given  by: 

yv+l=0.32Ev+,rL^ 


(12) 


(13) 


where  Fv+j,v  is  the  energy  of  the  v+ 1  v  transition  in  K,  L  is  the  characteristic  parameter  of  the 
short  range  repulsive  potential  in  A,  /u  is  the  reduced  mass  of  the  two  colliding  particles  in  a.u., 
and  Tg  is  the  gas  temperature  in  K.  We  have  taken  L=0.25  A  for  all  V-T  processes. 


The  function  F  in  Eqn.  (12)  is  given  by  [34]: 
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(14) 


The  calculated  forward  and  reverse  rates  are  plotted  in  Figure  22.  The  rates  for  Y-T  relaxa¬ 
tion  of  N2(v)  by  collision  with  O  atoms  are  between  two  and  three  orders  of  magnitude  higher  than 
the  rates  of  V-T  relaxation  of  N2OO  by  collision  with  N2. 
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Figure  22.  Rate  coefficient  for  V-T  relaxation:  N2  (X,v)  +M  —>N2  (X,v—1)  +  M,  with  M=0,  N2  at  Tg=2000K 


The  rates  k°J  corresponding  to  the  V-V  and  V-V’  processes  (9)  with  v/=0  and  V2=  0,  are 


calculated  using  the  following  expression: 


kio  -  ATg  exP 
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rrt  1/3 

Lg  J 


(15) 


where  the  gas  temperature  Tg  is  expressed  in  K  and  the  parameters  A,  B,  n  are  listed  in  Table  3. 


Table  3.  Parameters  for  the  V-V  and  V-V’  rates  k°0‘ 


V-V  or  V 

-V’  proc- 

A  ! 

B 

jt;: 

n2-n2 

1.27E-17 

0 

1.483 

n2-o2 

L23E-14 

104 

1 

n2-no 

4.22E-10 

86.35 

0 

The  N2-N2  V-V  exchange  process  has  been  studied  experimentally  at  300K  by  [36-41]  and 
theoretically  by  semiclassical  trajectory  methods  [42,  43]  and  SSH  theory  [44,  45],  These  results 
are  plotted  in  Figure  23. 
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Figure  23.  Rate  coefficients  for  V-V  relaxation:  N2  (X,l)  +  N2  (X, 0)  —>N2  (X,0)  +N2  (X,l) 

There  is  a  disagreement  of  about  one  decade  among  the  different  results,  and  the  reason  for 
the  discrepancy  remains  unclear.  The  theoretical  models  are  not  exact  and  contain  assumptions  that 
can  easily  change  the  calculated  rates  by  up  to  a  factor  of  5.  On  the  other  hand,  measuring  V-V 
probabilities  for  nitrogen  is  difficult  because  of  the  very  weak  absorption  vibrational  spectrum.  Po- 
dobedov  [36],  Valyanskii  et  al  [37,  38],  Suchkov  and  Shebeko  [39]  and  Gordeev  and  Shakatov  [41] 
monitored  the  relaxation  of  nitrogen  vibrational  levels  using  Raman  methods.  They  obtained  the 
rate  constant  by  fitting  their  experimental  results  with  solutions  of  the  coupled  kinetic  equa¬ 
tions  describing  the  relaxation  process.  The  uncertainties  of  the  rate  constants  obtained  depend  on 
the  accuracy  of  all  parameters  implemented  in  the  kinetic  equations.  However  accurate  V-V  values 
for  CO  are  available  and  have  been  satisfactorily  modeled  theoretically  by  Jeffers  and  Kelley  [46]. 
CO  and  N2  are  similar  in  that  they  have  the  same  number  of  electrons,  the  same  molecular  weight 
and  nearly  the  same  fundamental  vibrational  frequency.  The  most  significant  difference  between 
V-V  exchange  in  CO  and  N2  is  the  long-range  dipole-dipole  interaction  in  the  former  system.  If 
this  interaction  is  absent,  as  it  is  for  N2,  the  short-range  repulsive  interaction  dominates  the  V-V 
process.  The  contribution  of  the  short-range  repulsive  interaction  to  the  V-V  transfer  probabilities 
should  be  much  the  same  in  N2  as  in  CO.  We  have  done  a  calculation  for  N2  based  on  the  short- 
range  repulsive  interaction  model  presented  in  [46].  This  calculation  leads  to  rate  constants  of 
about  6xl0'14  cm3/s  at  300  K  and  about  lxlO'12  cm3/s  at  2000  K.  These  rate  constants  has  been  fit¬ 
ted  with  an  analytical  expression  of  the  form  given  by  Eqn.  (15),  with  the  parameters  reported  in 
Table  3. 
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The  N2-O2  V-V’  energy  exchange  has  been  studied  experimentally  for  temperatures  be¬ 
tween  100  and  4500  K  by  White  and  Millikan  [47,  48],  Bauer  and  Roesler  [49],  Breashears  and 
Bird  [22],  Kozlov  et  al.  [21],  and  Wilson  et  al.  [50].  The  different  sets  of  experimental  data  are  in 
good  agreement  with  one  another  and  are  well  fitted  by  the  analytical  expression  of  Kozlov  [21] 
which  reproduces  the  experimental  results  within  a  factor  of  2  for  temperatures  between  1000  and 
3000  K.  The  parameters  of  this  expression  are  reported  in  Table  3. 


The  N2-NO  V-V’  energy  exchange  has  been  studied  experimentally  for  temperatures  be¬ 
tween  300  and  3500  K  by  Callear  [51],  Taylor  et  al.  [52],  Stephenson  [53]  and  Doyennette  and 
Margotin-Maclou  [54].  The  analytical  expression  of  Taylor  [25],  with  parameters  reported  in  Table 
3,  reproduces  very  well  the  low  temperature  (300-400  K)  experimental  data  of  Refs.  [51,  53,  54].  It 
also  reproduces  the  high  temperature  (1500-3400  K)  data  of  Ref.  [25]  within  a  factor  of  three. 


The  rate  coefficients  for  exothermic  transitions  between  upper  vibrational  levels 


have  been  calculated  using  the  relation: 

%tl=CG(v1+l,vJ+l)  (16) 

where  G(vi+  1,V2+ 1)  is  an  appropriate  function  which  can  be  expressed  using  SSH  theory  [34]  and 
some  approximations  for  the  Morse  oscillator  model  as: 


G(v,  +1,v2+1)D 


1-*«i(v,+1)  l-*,2(v2+l)  F  ( Ti°o ) 


(17) 


where  xei  and  xe2  are  the  anharmonicities  of  the  two  molecules  involved.  F(y)  is  given  by  Eqn.  (14) 
with  defined  as: 


where  Ev  are  the  energies  of  the  initial  and  final  levels  in  K.  We  have  taken  L=0.25  A  for  all  V-V 
and  V-V’  processes.  Note  that  j/*;]2*1  must  be  always  positive  since  we  are  considering  the  reaction 
in  the  exothermic  direction. 


The  calculated  forward  and  reverse  rates  are  plotted  in  Figure  24  as  a  function  of  the  vibra¬ 
tional  number  v/  for  V2=  0.  For  the  N2-O2  process,  the  rates  increase  up  to  the  resonance  point  at 
v/= 27,  and  decrease  after  this  value.  We  observe  the  same  behavior  for  the  N2-NO  process  but  the 
resonance  appears  at  a  lower  value  of  v/  (v/=16)  because  the  spacing  between  NO  levels  is  larger 
than  between  O2  levels.  For  the  N2-N2  process,  the  rates  increase  until  v/=J  and  then  decrease  be¬ 
cause  of  the  increasing  vibrational  energy  gap  between  the  two  N2  molecules. 
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Figure  24.  Rate  coefficients  for  V-V  and  V-V’  exchange  :  N2  (X,v,+I)  +AB(X,0)  —>  N2  (X,v,)  +AB(X,1),  with 

AB  =  N2,02  and  NO,  at  Tg=2000K. 

4.3.2.C.  Vibrational  distribution  of  N?  ground  state 

The  V-T,  V-V  and  V-V’  rates  of  the  foregoing  section  were  implemented  in  the  model  and 
the  vibrational  distribution  of  the  N2  ground  and  excited  electronic  states  was  determined  by  solving 
a  system  of  kinetic  equations  at  steady  state  in  which  the  vibrational  levels  of  the  N2  ground  and 
excited  electronic  states  are  the  unknowns. 

The  total  concentration  of  N2  was  determined  with  the  two-temperature  kinetic  model  [13] 
and  fixed  by  replacing  the  vibrational  level  v=0  of  ground  electronic  state  by  the  mass  conservation 
equation. 

The  total  populations  of  the  other  species  were  fixed  and  determined  with  the  two- 
temperature  kinetic  model,  and  their  internal  distribution  was  calculated  according  to  a  Boltzmann 
distribution  at  the  vibrational  temperature  Tv=Tg  and  at  the  electronic  temperature  Tei=Te.  We  now 
present  our  calculations  of  the  vibrational  distribution  of  the  N2  ground  state  at  Tg=2000K  and  for 
different  electron  temperatures. 

For  electron  temperatures  Te  lower  than  6000K,  the  vibrational  distribution  is  very  close  to  a 
Boltzmann  distribution  at  the  gas  temperature  Tg=2000K.  Figure  25  and  Figure  26  show  the  calcu- 
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lated  vibrational  distributions  for  a  gas  temperature  of  2000K  and  an  electron  temperature  of  9000K 
and  16000K  respectively.  The  Boltzmann  distributions  at  Tv=Tg  and  Tv=Te  are  also  shown  on  these 
figures. 


For  Te=9000K,  the  vibrational  excitation  introduced  by  VE  transfer  is  mainly  redistributed 
via  V-T  relaxation  of  N2  by  collision  with  N2,  and  via  N2-N2  V-V  exchange.  The  N2-O2  and  N2-NO 
V-V’  processes  do  not  significantly  affect  the  populations  of  N2  levels.  We  checked  that  this  con¬ 
clusion  remains  valid  if  we  assume  a  different  internal  distribution  for  the  O2  and  NO  molecules. 
Indeed,  the  rates  of  N2-NO  exchange  are  faster  than  those  of  N2-N2  exchange  above  v=3  (see  Figure 
24)  but  the  total  concentration  of  NO  is  two  orders  of  magnitude  lower  than  the  concentration  of  N2 
and  the  rates  for  N2-O2  V-V  exchange  are  fast  for  v>20  but  the  population  of  those  levels  is  mainly 
governed  by  V-T  transfer  processes.  The  vibrational  distribution  at  Te=9000K  lies  between  the 
Boltzmann  distributions  at  Tv=Tg  and  Tv=Te,  but  remains  closer  to  a  distribution  at  Tv=Tg. 


Figure  25.  N2(X,v)  vibrational  distribution  function  at  TS=2000K  and  Te=9000K,  P=latm 

For  Te=T6000K,  almost  25%  of  the  O2  molecules  are  dissociated  and  the  vibrational  excita¬ 
tion  is  mainly  redistributed  by  V-T  relaxation  of  N2  by  collision  with  O  atoms  and  by  N2-N2  V-V 
exchange.  At  this  electron  temperature,  the  vibrational  distribution  of  the  first  15  levels  is  close  to 
the  Boltzmann  distribution  at  Tv=Te. 
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Figure  26.  N2(X,v)  vibrational  distribution  function  at  TS=2000K  and  Te=16000K,  P=latm. 


4.3.2.d.  Inelastic  electron  energy  losses  in  air  plasmas 

Electron  inelastic  energy  losses  can  now  be  calculated  by  summing  the  contributions  of  all 
electron  impact  collisional  processes: 


Qm 


processes 


(19) 


where  Ej-E,  represents  the  internal  energy  gained  by  heavy  species  during  the  collision  (E/must  be 
greater  than  E I)  and  dnf  /  dt  is  the  net  volumetric  rate  of  production  of  heavy  species  in  the  final 
energy  level  f  In  an  atmospheric  pressure  air  plasma  characterized  by  a  gas  temperature  between 
1000  and  3000  K  and  electron  temperatures  up  to  17000  K,  the  dominant  contribution  to  electron 
inelastic  energy  losses  is  the  electron-impact  vibrational  excitation  of  N2  ground  state.  The  electron 
impact  vibrational  excitation  cross-sections  of  O2  and  NO  ground  states  are  two  orders  of  magni¬ 
tude  lower  than  those  of  N2,  and  therefore  the  contribution  of  these  molecules  is  negligible. 


The  total  rate  of  energy  loss  can  be  expressed  as: 

Qe-V  ~  £?v,v2 


(20) 


Vj  V2>Vj 


where  v/  and  V2  are  the  initial  and  final  vibrational  levels  of  the  transition,  and  where  the  elementary 
rate  Qm  is  written  as: 

6v,v2  =(km  [^(X.v,)]  — ^  [N2{X,v2)])neAEVi 


(21) 
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In  Eqn.  (21),  ne  is  the  concentration  of  electrons,  km  and  km  are  the  excitation  and  deexcitation 
rate  coefficients  and  A E  stands  for  the  difference  of  energy  between  the  two  vibrational  levels  V2 
and  vy. 


0  depends  strongly  on  the  N2  ground  state  internal  distribution.  Vibrational  population 
distributions  calculated  with  the  method  presented  in  the  foregoing  section  are  used  in  Eqn.  (21)  to 
determine  the  electron  inelastic  energy  losses.  The  predicted  inelastic  power  losses  are  shown  in 
Figure  27.  At  low  electron  temperatures  and  densities,  the  vibrational  levels  of  N2  ground  state  are 
close  to  a  Boltzmann  distribution  at  Tv=Tg.  The  excited  vibrational  levels  have  low  population. 
Therefore,  the  power  lost  by  e-V  excitation  is  not  balanced  by  the  power  regained  from  V-e  supere¬ 
lastic  de-excitation.  As  the  electron  temperature  increases,  the  electron  density  also  increases  and 
eventually  the  vibrational  population  distribution  tends  toward  Tv=Te .  The  net  power  losses  do  not 
increase  as  rapidly  because  of  the  increased  importance  of  superelastic  collisions. 


Figure  27.  Predicted  inelastic  electron  power  losses  in  atmospheric  pressure  air  at  2000  K. 


It  is  sometimes  convenient  to  define  an  electron  “energy  loss  factor”  as  the  ratio  of  total 
(elastic  +  inelastic)  energy  losses  to  the  elastic  energy  losses: 

f  _  Qi'l  Qinel_ 

e~  a, 


(22) 
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§  -  Qd  -  Odd  (22) 

Qe, 

where  Qel  is  the  volumetric  power  lost  by  free  electrons  through  elastic  collisions,  and  Qe,  is  the 
sum  of  contributions  of  collisions  between  electrons  and  heavy  species  h=N2, 02  and  O: 

Qel=^k(Te-Th)^Veh  (23) 

*  mh 

In  Eqn.  (23),  k  is  the  Boltzmann  constant,  me  and  mh  are  the  masses  of  electron  and  heavy  species 
respectively,  7*  is  the  kinetic  temperature  of  the  heavy  species  (equal  to  Tg),  and  veh  represents 
the  average  frequency  of  collisions  between  the  electrons  and  heavy  particle  h.  veh  can  be  ex¬ 
pressed  in  terms  of  the  number  density  of  neutral  species  nh,  the  electron  velocity 
ge  =  ^8kTe  /  nme  and  the  average  elastic  collision  cross-section  QeJh : 


o  1 — . - ■ - * - * - 1 - * - - - - - • - 

4000  9000  14000 

electron  temperature  Te  (K) 

Figure  28.  Energy  loss  factor  Se  at  Tg=1800K,  2000,  and  2900K,  as  a  function  of  Te 

Figure  28  shows  the  calculated  electron  energy  loss  factor  as  a  function  of  the  electron 
temperature  for  two  values  of  the  gas  temperature,  Tg  =1800  and  2900  K.  As  can  be  seen  from 
this  figure,  the  inelastic  loss  factor  is  a  relatively  weak  function  of  the  gas  temperature.  It  in¬ 
creases  up  to  Te=8000  K  as  the  net  rate  of  production  of  N2  molecules  in  vibrational  level  v2>vj 
increases  with  Te,  and  then  decreases  due  to  the  transition  Tv=Tg  to  Tv=Te.  When  Tv  becomes 
close  to  Te,  the  forward  and  reverse  rates  are  practically  balanced  and  the  net  rate  of  energy  lost 
by  VE  transfer  approaches  zero. 
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4.3.2.e.  Predicted  DC  Discharge  Characteristics  in  Atmospheric  Pressure  Air 

The  results  of  the  previous  subsections  enable  us  to  convert  the  “S-shaped”  curve  of  ne  vs. 
Te  into  electric  field  vs.  current  density  discharge  characteristics.  This  result  is  obtained  by  com¬ 
bining  Ohm’s  law  and  the  electron  energy  equation.  The  latter  incorporates  the  results  of  the  col- 
lisional-radiative  model  to  account  for  non-elastic  energy  losses  from  the  free  electrons  to  the  mo¬ 
lecular  species.  The  predicted  discharge  characteristics  for  atmospheric  pressure  air  at  2000  K 
are  shown  in  Figure  29.  These  discharge  characteristics  exhibit  variations  that  reflect  both  the  S- 
shaped  dependence  of  electron  number  density  versus  Te,  and  the  dependence  of  the  inelastic  en¬ 
ergy  loss  factor  on  the  electron  temperature  and  number  density.  We  have  used  these  predicted 
characteristics  as  a  starting  point  to  design  the  DC  glow  discharge  experiments  presented  in  Sec¬ 
tion  4.4.  If  these  predictions  are  correct,  the  production  of  1013  electron/cm3  requires  an  electric 
field  of  ~1.35kV/cm,  and  a  current  density  of  ~10.4A/cm2.  Thus  the  power  required  to  produce 
1013  electrons/cc  in  air  at  2000  K  is  approximately  14kW/cm3. 


Figure  29.  Predicted  discharge  characteristics  for  atmospheric  pressure  air  at  2000  K. 
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density  versus  electric  field,  direct-current  (DC)  glow  discharged  have  been  produced  in  flowing 
low  temperature,  atmospheric  pressure  air  plasmas.  The  flow  velocity  is  around  400  m/s,  and  the 
gas  temperature  is  varied  between  1800  and  2900  K.  These  experiments  show  that  it  is  feasible  to 
create  stable  diffuse  glow  discharges  with  electron  number  densities  in  excess  of  1012  cm-3  in  at¬ 
mospheric  pressure  air  plasmas.  Electrical  characteristics  were  measured  and  the  thermodynamic 
parameters  of  the  discharge  were  obtained  by  spectroscopic  measurements.  The  measured  gas  tem¬ 
perature  is  not  noticeably  affected  by  whether  or  not  the  DC  discharge  is  applied.  The  discharge 
area  was  determined  from  spatially  resolved  optical  measurements  of  plasma  emission  during  dis¬ 
charge  excitation.  The  measured  discharge  characteristics  are  compared  in  Section  4.5  with  the 
predicted  electrical  characteristics. 


4.4.2.  DC  discharge  experimental  setup 

The  ionization  process  in  the  discharge  region  is  accompanied  by  energy  transfer  to  the  gas 
through  collisions  between  electrons  and  heavy  particles.  Electrons  lose  more  than  99.9%  of  the 
energy  gained  from  the  electric  field  to  molecular  N2  through  vibrational  excitation,  and  the  vibra- 
tionally  excited  N2  transfers  energy  to  translational  modes  through  vibrational  relaxation.  Thus  the 
degree  of  gas  heating  (ATg)  is  a  function  of  the  volumetric  power,  jE,  deposited  into  the  plasma  by 
the  discharge  and  the  competition  of  the  vibrational  relaxation  time  and  the  residence  time  r  of  the 
plasma  in  the  discharge  region.  To  limit  gas  heating  to  acceptable  levels  for  given  volumetric 
power,  it  is  desirable  to  flow  the  plasma  at  high  velocity  through  the  discharge  region. 

The  experimental  set-up  is  shown  schematically  in  Figure  30.  Atmospheric  pressure  air  is 
heated  with  a  50  kW  radio-frequency  inductively-coupled  plasma  torch  operating  at  a  frequency  of 
4  MHz.  A  2  cm  exit  diameter  nozzle  is  mounted  at  the  exit  of  the  torch  head.  The  flow  rate  in¬ 
jected  in  the  torch  was  approximately  96  standard  liters  per  minute  (slpm)  (64  slpm  radial  and 
32  slpm  swirl)  and  the  plate  power  settings  were  8.9  kV  x4.1  A,  with  approximately  14  kW  of 
power  deposited  into  the  plasma.  Under  these  conditions,  the  temperature  of  the  plasma  at  the  exit 
of  the  2  cm  diameter  nozzle  is  about  5000  K  and  its  velocity  is  ~  100  m/s. 
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power  deposited  into  the  plasma.  Under  these  conditions,  the  temperature  of  the  plasma  at  the 
exit  of  the  2  cm  diameter  nozzle  is  about  5000  K  and  its  velocity  is  ~  100  m/s. 
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Figure  30.  Set-up  (not  to  scale)  for  discharge  experiments  showing  the  torch  head,  the  injection  ring,  the  2-cm 
diameter,  18-cm  long  water-cooled  quartz  mixing  test-section,  the  2  — >  1  cm  converging  nozzle,  electrodes, 
voltage  pins,  and  electrical  circuit. 


The  plasma  then  enters  a  quartz  test-section  where  it  is  cooled  to  the  desired  temperature 
by  mixing  with  an  adjustable  amount  of  cold  air  injected  into  the  plasma  stream  through  a  radial 
mixing  ring.  The  quartz  test-section  length  of  18  cm  ensures  that  the  flow  residence  time  (ap¬ 
proximately  1.6  ms  here)  is  greater  than  the  characteristic  time  for  chemical  and  thermal  equili¬ 
bration  of  the  plasma  (<  1ms).  Thus  at  the  exit  of  the  quartz  test-section  the  air  flow  is  close  to 
Local  Thermodynamic  Equilibrium  (LTE)  conditions.  Finally,  a  1  cm  exit  diameter  converging 
water-cooled  copper  nozzle  is  mounted  at  the  exit  of  the  mixing  test-section.  This  nozzle  is  used 
to  control  the  velocity,  hence  the  residence  time,  of  the  flow  within  the  discharge  region.  Two- 
dimensional  computational  fluid  dynamics  (CFD)  calculations  performed  at  the  University  of 
Minnesota  show  that  the  axial  velocity  at  the  entrance  of  the  discharge  region  is  approximately 
445  m/s  [55],  The  discharge  itself  is  produced  between  two  platinum  pin  electrodes  of  0.5  mm 
diameter  held  along  the  axis  of  the  air  stream  by  two  water-cooled  1/16”  stainless-steel  tubes 
placed  crosswise  to  the  plasma  flow.  The  bottom  electrode  is  mounted  on  the  copper  nozzle  and 
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The  discharge  was  driven  by  a  Del  Electronics  Model  RHSV10-2500R  power  supply  with 
reversible  polarities,  capable  of  operation  in  control  current  or  control  voltage  mode,  with  current 
and  voltage  outputs  in  the  ranges  of  0  -  250  mA  and  0-10  kV,  respectively.  For  the  present  ex¬ 
periments,  the  cathode  (bottom  electrode)  was  biased  to  negative  potentials  with  respect  to  ground. 

The  electric  field  within  the  discharge  region  is  measured  from  the  potential  on  a  high  purity 
platinum  wire  (0.02  inch  diameter)  that  extends  to  the  center  of  the  discharge  region.  The  platinum 
wire  is  held  by  a  small  ceramic  tube  installed  on  a  two-way  (horizontal  and  vertical)  translation 
stage.  Horizontal  translation  moves  the  pin  into  the  discharge  region  for  electric  field  measure¬ 
ments,  and  out  of  the  discharge  during  spectral  emission  measurements.  Vertical  translation  moves 
the  pin  along  the  discharge  axis  to  determine  the  electric  field  from  potential  measurements.  Al¬ 
though  pure  platinum  melts  at  ~  2045  K,  radiation  cooling  prevents  melting  of  the  platinum  wires 
for  plasma  temperatures  up  to  at  least  3000  K.  The  voltage  measurements  reported  here  were  made 
with  a  Tektronix  Model  P6015A  high  voltage  (20  kV  DC,  40  kV  peak  pulse)  probe  and  a  Hewlett 
Packard  Model  545 10 A,  digitizing  oscilloscope.  The  current  was  measured  from  the  voltage  drop 
across  the  12  kQ  ballast  resistor  of  the  DC  circuit. 

The  setup  for  optical  emission  spectroscopy  diagnostics  includes  a  SPEX  model  750M,  3/4 
meter  monochromator  fitted  with  a  1200  lines/mm  grating  blazed  at  200  nm  and  a  backthinned 
Spectrum  One  thermoelectrically  cooled  charge-coupled  device  (CCD)  camera.  The  30  x  12  mm 
CCD  chip  contains  2000  x  800  pixels  of  dimension  15x15  pm.  The  dispersion  of  the  optical  sys¬ 
tem  is  -  1.1  nm/mm.  The  monochromator  entrance  slit  width  was  set  at  200  pm,  and  26  columns  of 
800  pixels  were  binned  to  produce  an  equivalent  exit  slit  width  of  390  pm.  The  spatial  resolution 
was  ~0.5  mm  as  determined  by  the  monochromator  entrance  slit  width  and  the  magnification  of  the 
optical  train  (2.5  for  two  lenses  of  focal  length  50  and  20  cm).  Absolute  spectral  intensity  calibra¬ 
tions  were  obtained  with  an  Optronics  model  OL550  tungsten  filament  lamp  and  a  1  kW  argon  ar- 
cjet,  with  radiance  calibrations  traceable  to  National  Institute  of  Standards  and  Technology  (NIST) 
standards. 


4.4.3.  Spectroscopic  Measurements 

4.4.3 . a.  Measurements  without  DC  discharge  applied 

The  gas  temperature  (rotational  temperature)  without  DC  discharge  applied  was  measured 
by  emission  spectroscopy  of  the  OH  (A  -^X)  transition.  The  OH  (X  -»  20  transition  is  one  of  most 
intense  emission  systems  in  low  temperature  (T  <  4000  K)  air  plasmas  containing  even  a  small 
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4.4.3.  Spectroscopic  Measurements 

4.4.3.  a.  Measurements  without  DC  discharge  applied 

The  gas  temperature  (rotational  temperature)  without  DC  discharge  applied  was  measured 
by  emission  spectroscopy  of  the  OH  (A  — >  X)  transition.  The  OH  (A  — >  X)  transition  is  one  of 
most  intense  emission  systems  in  low  temperature  ( T  <  4000  K)  air  plasmas  containing  even  a 
small  amount  (~  1%)  of  H2  or  H20.  In  the  present  experiments,  the  water  content  of  the  air  in¬ 
jected  into  the  torch  was  sufficient  to  produce  intense  OH  radiation.  Rotational  temperatures  can 
be  obtained  from  the  relative  intensities  of  two  groups  of  rotational  lines  ,  corresponding  to  the  R 
(Ri  and  R2)  and  P  (Pi  and  P2)  branches  of  the  (0,0)  vibrational  band,  that  form  distinct  peaks  at 
306.8  and  309.2  nm,  respectively.  Figure  31  shows  a  simulation  of  the  OH  A  — >  X  bands  ob¬ 
tained  with  the  NEQAIR2  radiation  code  [1]  which  incorporates  an  accurate  model  for  this 
transition  [56].  One  advantage  of  this  rotational  temperature  measurement  technique  is  that  it 
does  not  require  the  measured  signal  to  be  corrected  for  the  spectral  intensity  response  of  the 
detection  system  because  the  response  of  the  detection  system  used  is  nearly  constant  over  the 
less  than  3  nm  spectral  range  of  interest.  The  relative  intensity  of  the  two  peaks  is  fairly  sensitive 
to  the  rotational  temperature,  as  can  be  seen  from  Figure  32. 


Wavelength  (nm) 

Figure  31.  NEQAIR2  simulations  of  the  OH  A  -»  X  emission  spectrum  as  a  function  of  temperature  in  LTE 

air,  normalized  to  the  peak  of  the  P  branch. 
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Figure  32.  Intensity  ratio  of  the  P/R  band  head  as  a  function  of  the  rotational  temperature  (NEQAIR2  simula¬ 
tions  for  an  instrumental  slit  function  of  of  HWHM=0.21  nm). 

Line-of-sight  OH  emission  spectra  were  recorded  with  the  discharge  off.  The  amount  of 
cold  air  mixing  was  adjusted  to  vary  the  temperature  of  the  preheated  air.  The  measured  OH 
spectra  were  later  fitted  with  NEQAIR2.  As  shown  in  Figure  33,  the  gas  temperature  can  be  var¬ 
ied  from  1800  to  2900  K  by  adjusting  the  amount  of  cold  air  mixing. 
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Figure  33.  Measured  OH  A  —>X  emission  spectra  without  discharge  applied  as  a  function  of  the  amount  of 

cold  air  mixing. 


Absolute  spectral  measurements  were  made  of  plasma  emission  between  200  and  400  nm 
in  order  to  verify  that  the  preheated  air  is  chemical  or  thermal  equilibrium.  Comparisons  are 
shown  in  Figure  34  with  spectral  simulations  made  with  the  NEQAIR2  code,  assuming  LTE 
conditions,  and  using  the  measured  temperature  profile  and  corresponding  equilibrium 
concentrations  for  the  case  of  145  slpm  mixing  cold  air.  As  can  be  seen  from  Figure  34,  the 
numerical  LTE  simulations  reproduce  the  experimental  spectrum  both  in  shape  and  absolute 
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duce  the  experimental  spectrum  both  in  shape  and  absolute  intensity,  thus  indicating  that  the  pre¬ 
heated  air  (with  no  discharge  applied)  is  close  to  LTE.  Thus  the  electron  number  density  should  be 
close  to  the  equilibrium  value  when  the  discharge  is  not  applied.  At  2000  K  for  instance,  this  equi¬ 
librium  value  is  ~3xl06  cm'3. 


Wavelength  (nm) 

Figure  34.  Measured  emission  spectra  without  discharge  applied,  and  comparisons  with  NEQAIR2  spectral  simu¬ 
lations  assuming  LTE. 


4.4.3 .b.  Measurements  with  DC  discharge  applied 

Emission  spectra  were  also  measured  with  the  discharge  applied  (discharge  current  of 
150  mA),  and  results  are  shown  in  Figure  35.  As  can  be  seen  in  the  figure,  a  factor  of  ~  104  en¬ 
hancement  of  the  emission  due  to  NO  gamma  (A  — >  X)  and  a  factor  of  ~  105  enhancement  of  the 
emission  due  to  N2  (C  — »  B )  bands  were  observed.  Figure  35  also  shows  that  the  N2  second  positive 
system  (C  -»  B)  bands  overlap  the  OH  (A  —>  X)  feature  around  308  nm.  This  overlap  precludes  the 
accurate  measurements  of  the  rotational  temperature  from  the  OH  (A  —>X)  transition.  Therefore  the 
gas  temperature  (rotational  temperature)  was  measured  by  means  of  emission  spectroscopy  of  the 
(0,0)  band  of  the  N2  second  positive  system  -  the  N2  (C  — >  B)  transition.  Line-of-sight  N?  emission 
spectra  were  recorded  along  lateral  chords  of  the  plasma.  The  spectra  were  fitted  with  NEQAIR2  to 
obtain  a  rotational  temperature  Tr.  Figure  36  shows  a  typical  experimental  spectrum  and  NEQAIR2 
fit  for  150  mA  discharge  current  and  145  slpm  of  cold  air  mixing. 
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Figure  35.  Line-of-sight  emission  spectra  measured  at  a  discharge  current  /  =  150  mA. 


Figure  36.  Emission  Spectrum  of  N2  second  positive(C-»l?)  bands  and  comparisons  with  NEQAIR2  spectral 

simulations. 

The  vibrational  structure  of  the  N2  second  positive  (C  — >  B)  system  can  be  used  to  determine 
the  vibrational  temperature  of  the  C  state.  As  shown  in  Figure  36,  the  (0,2),  (1,3),  (2,4)  and  (3,5) 
vibrational  sequences  were  utilized  here.  By  matching  NEQAIR2  spectra  (calculated  with 
Tr  =  Tr, measured)  to  the  shape  of  the  measured  (0,2),  (1,3),  (2,4)  and  (3,5)  bands  spectra,  the  N2  C  state 
vibrational  temperature  Tv,meaSured  was  obtained. 


The  absolute  intensity  of  the  spectrum  can  then  be  used  to  determine  the  population  of  the 
N2  C  electronic  state.  As  discussed  in  Ref.  [1],  the  integrated  intensity  of  a  spectral  line  is  defined 
as 


/[w/cm3  /sr]=— 


4  % 


(25) 


53 


where  Au/  (s'1)  is  the  Einstein  A  coefficient  representing  the  probability  of  a  radiative  transition  from 
level  u  to  level  /,  A Eul  is  the  energy  difference  between  the  level  u  and  /,  and  nu  is  the  population  of 
electronic  level  u.  The  intensity  of  the  transition  gives  a  measurement  of  the  number  density  of  the 
upper  electronic  energy  level,  which  can  then  be  used  to  determine  an  electronic  temperature  7W  de¬ 
fined  by  the  Boltzmann  relation: 

(26) 

where  £u  is  the  electronic  energy  of  level  u;  kB  is  Boltzmann’s  constant;  no  is  the  population  of  the 
ground  state;  gu  and  go  are  the  degeneracies  of  the  upper  and  ground  level,  respectively. 


—  =  — expf^^- 

n0  So  ^  klFel 


The  absolute  (0,0)  bandhead  intensity  lmodei  can  be  obtained  from  NEQAIR2  by  assuming 
T,  -  Tr, measured  Tv  =  Tv, measured,  and  Te\  =  Tr.  Thus  from  equations  (25)  and  (26)  we  obtain  the  meas¬ 
ured  N2  C  state  electronic  temperature  Tei,m(o  via 


Additional  discharge  experiments  were  conducted  with  different  gas  temperatures.  Figure 
37  shows  the  measured  N2  second  positive  system  spectra  and  their  rotational  and  vibrational  tem¬ 
peratures  at  a  discharge  current  of  150  mA,  as  a  function  of  the  amount  of  mixing  air.  It  can  be 
seen  in  the  figure  that  both  the  rotational  and  vibrational  temperatures  are  lower  with  higher  amount 
of  mixing  cold  air.  Figure  38  shows  the  measured  spectrum  as  a  function  of  the  discharge  current 
for  145  slpm  of  mixing  air.  As  can  be  seen  from  the  figure,  the  rotational  temperature  remains  the 
same  at  all  currents,  but  the  vibrational  temperature  increases  with  increasing  discharge  current. 


Figure  37.  Measured  N2  second  positive(C  — >  B)  bands  with  discharge  on  as  a  function  of  the  amount  of  cold  air 

mixing  (/=  150  mA). 
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Figure  38.  Measured  N2  second  positive(C  — >  B)  bands  with  discharge  on  as  a  function  of  discharge  current  for 

the  case  of  145  slpm  cold  air  mixing. 


Radial  rotational  temperature  profiles  with  and  without  the  discharge  applied  were  measured 
along  chords  of  the  plasma  from  Abel-inverted  N2  second  positive  system  emission  spectra  and  OH 
emission  spectrum,  respectively.  Figure  39  shows  the  measured  radial  temperature  profiles  at  a  dis¬ 
tance  of  1 .5  cm  downstream  of  the  cathode(i.e.  along  midway  between  the  two  electrodes).  As  can 
be  seen  from  the  figure,  the  applied  discharge  does  not  noticeably  increase  the  rotational  tempera¬ 
ture  of  the  plasma  at  this  location.  Figure  40  shows  the  radial  N2  C  state  electronic  and  vibrational 
temperature  profiles.  On  the  axis  of  the  discharge,  the  electronic  temperature  of  the  N2  C  state 
reaches  about  5000  K,  and  the  vibrational  temperature  is  about  3000  K. 


Figure  39.  Rotational  temperature  profiles  with  and  without  the  applied  DC  discharge  at  1.5  cm  downstream  of 
the  bottom  electrode.  The  temperature  profile  without  the  discharge  was  measured  from  rotational  lines  of  the 
OH  (A  -4  A)  transition.  With  the  discharge  applied,  the  rotational  temperature  was  measured  from  lines  of  the  N2 
(C  -4  B )  transition  in  the  ultraviolet. 
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Figure  40.  Electronic,  vibrational  and  rotational  temperature  profiles  of  N2  C  state  with  an  applied  discharge 

current  of  I- 150  mA. 


Figure  41a  shows  a  photograph  of  the  air  plasma  plume  at  a  temperature  of  approximately 
2200  K  in  the  region  between  the  two  electrodes  without  the  discharge  applied.  Figure  41b  shows 
the  same  region  when  a  DC  discharge  of  5.2  kV  and  200  mA  is  applied  between  the  two  electrodes. 
In  these  experiments,  the  distance  between  electrodes  is  3.5  cm.  The  bright  region  in  Figure  41b 
corresponds  to  the  discharge-excited  plasma.  Thus  the  plasma  plume  without  discharge  applied  ap¬ 
pears  to  be  homogeneous  over  a  larger  diameter  than  the  plasma  plume  with  the  discharge  applied. 
However,  Figure  39  showed  that  the  gas  temperature  profile  is  practically  the  same  as  in  the  dis¬ 
charge  applied  case.  The  increased  brightness  in  Figure  41b  is  due  to  the  emission  of  excited  elec¬ 
tronic  states  of  molecular  NO  and  N2  (see  Figure  40).  Thus  the  applied  discharge  increases  excited 
state  populations  without  significantly  increasing  the  gas  temperature. 


(a)  (b) 

Figure  41.  (a)  Air  plasma  at  2000  K  without  electrical  discharge,  (b)  Air  plasma  at  2000  K  with  applied  dis¬ 
charge  (1.4  kV/cm,  200  mA).  Interelectrode  distance  =  3.5  cm.  The  measured  electron  number  density  in  the 
bright  discharge  region  is  around  1012  cm'3. 
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4.4.4.  Current  Density  Measurements 


The  current  density  at  the  center  of  the  plasma  was  determined  by  dividing  the  measured 
current  by  the  effective  discharge  area  A*,  i.e.  j(r  =  0)  =  HA* .  The  effective  discharge  area  is  ob¬ 
tained  from  the  following  relation, 

R 

[. j(r  =  0)x/]=  |2 7trj(r)dr ,  (28) 

o 


where  j(r)  is  the  local  current  density.  As  shown  in  Ref.  [13],  j(r)  is  approximately  proportional  to 
ne  (r).  Thus,  A*  can  be  calculated  as 


A*  = 


jne  {r)2nrdr 


(29) 


In  separate  discharge  experiments  conducted  with  a  nitrogen  plasma  [57],  the  electron  num¬ 
ber  density  profile  ne  was  measured  using  various  techniques  (from  Hp  Stark  broadening  and  N2  1st 
positive  emission  spectra)  to  calculate  A*  using  Equation  (29).  The  discharge  area  A*  was  also  es¬ 
timated  from  the  full  width  at  half  maximum  (FWHM)  of  the  N2  C-B  (0,0)  bandhead  intensity 
profile.  In  these  nitrogen  discharge  experiments,  the  effective  area  A*  obtained  with  Equation  (29) 
and  the  measured  ne  profile  was  found  to  be  equal  to  the  effective  area  obtained  from  the  N2  C  -  B 
(0,0)  emission  intensity  measurement.  Thus  for  the  present  air  plasma  discharge  experiments,  we 
estimate  the  effective  discharge  area  from  the  spatially  resolved  optical  measurements  of  N2  C  state 
emission.  Spectroscopic  measurements  of  N2  C  -  B  (0,0)  emission  with  the  applied  discharge  are 
shown  in  Figure  42.  It  can  be  seen  from  the  figure  that  the  diameter  (FWHM)  of  the  discharge  is 
approximately  3.2  mm.  This  diameter  was  monitored  and  found  to  be  constant  for  all  discharge 
currents  ranging  from  5  to  250  mA.  The  discharge  diameter  was  therefore  taken  to  be  3.2  mm  and 
assumed  constant  along  the  axis  of  the  discharge. 


Radius  (mm) 


Figure  42.  Spatial  extent  of  the  plasma  produced  by  the  discharge. 
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4.4.5.  Electric  Field  Measurements 

Electrode  and  pin  potentials  were  measured  as  a  function  of  the  applied  discharge  current 
which  was  varied  from  0  to  250  mA.  The  cathode  current  was  measured  from  the  voltage  drop 
across  the  12  kfi  ballast  resistor  placed  in  series  with  the  discharge  (see  Figure  30).  There  is  a 
small  difference  of  7  mA  between  the  anode  and  the  cathode  currents  that  was  found  to  be  due  to  a 
current  leak  through  the  water  cooling  circuit  of  the  anode.  All  results  reported  below  are  based  on 
the  measured  cathode  current,  which  is  not  affected  by  current  losses  to  the  cooling  circuit. 

Figure  43  shows  the  measured  pin  voltage  as  a  function  of  the  applied  current  along  the  axis 
of  the  discharge.  The  potential  varies  approximately  linearly  along  the  axis  of  the  discharge,  indi¬ 
cating  that  the  electric  field  is  approximately  uniform  in  the  discharge  region.  The  electric  field 
measurements  reported  here  were  determined  from  the  slope  of  a  linear  fit  of  the  pin  potentials.  In 
the  vicinity  of  the  cathodes,  voltage  falls  of  up  to  several  hundred  volts  were  observed.  These  val¬ 
ues  are  typical  of  the  cathode  fall  voltage  in  glow  discharges  [58]. 

The  total  voltage  across  the  discharge  was  also  measured  as  a  function  of  electrode  separa¬ 
tion,  by  translating  the  top  electrode  (anode)  vertically.  The  voltage-length  characteristic  for  a  dis¬ 
charge  current  of  150  mA  is  shown  in  Figure  44.  The  lowest  voltage  reading  as  the  electrodes  are 
brought  within  less  than  0.2  mm  from  one  another  provides  an  approximation  to  the  discharge  volt¬ 
age  at  zero  gap  length  [59,  60].  The  value  of  this  voltage  is  found  to  be  285  V  and  is  independent 
of  the  current  in  the  current  range  investigated  (10  -  250  mA).  This  value  agrees  with  the  cathode 
fall  voltage  reported  in  the  literature  [7,  61]  for  glow  discharges  in  air  with  a  platinum  cathode.  The 
voltage  gradient  in  the  positive  column,  given  by  the  slope  of  the  voltage-length  characteristic,  is 
constant  as  the  discharge  length  is  increased.  For  a  discharge  current  of  150  mA,  the  gradient  is 
about  1400  V/cm  (see  Figure  44)  and  is  consistent  with  the  electric  field  value  determined  from  the 
pin  measurements. 
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Distance  along  Discharge  Axis  (cm) 

Figure  43.  Measured  potentials  as  a  function  of  applied  current  in  the  discharge  section. 


Figure  44.  Voltage-length  characteristic  in  the  discharge  region. 


4.5.  Electrical  Characteristics  and  Power  Requirements  of  DC  Discharges  in  Air 

The  experimental  discharge  characteristics  presented  in  Section  4.4  for  plasma  temperatures 
ranging  from  1800  to  2900  K  are  shown  in  Figure  45.  They  are  also  compared  with  the  predicted 
characteristics  at  the  corresponding  gas  temperatures.  The  method  employed  to  predict  the  dis¬ 
charge  characteristics  was  discussed  in  Section  4.3.  As  can  be  seen  from  Figure  45,  good  agree¬ 
ment  is  obtained  between  the  measured  and  predicted  discharge  characteristics  over  a  range  of  ex¬ 
periments  spanning  over  three  orders  of  magnitude  in  current  density. 

Figure  45  also  shows  (dashed  curve)  the  resistive  characteristic  of  equilibrium  air  at  2900  K, 
given  by  the  relation: 
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E  = - - ,  (30) 

® equilibrium, 2900 K 

where  a  equilibrium,  2900  k,  the  electrical  conductivity  of  equilibrium  air  at  2900  K,  is  calculated  as: 


nequil,2900Ke2 


9  equilibrium, 2900K 


(31) 


where  neJiu,l>2900K  -  4xl010  cm'3  and  the  collision  frequency  ve_air  is  well  approximated  by 

r  \ 

P 


v  .  = 


ge(l.5xl0-'5[cm2]). 


(32) 


where  p  is  the  pressure  (1  atm),  Tg  =  2900  K  is  the  gas  temperature,  and  ge  =  yJSkTe  1  nme  is  the 
electron  thermal  velocity.  For  Tg  =  2900  K,  as  can  be  seen  from  Figure  45  the  predicted  E  vs.  j 
characteristic  is  close  to  the  resistive  equilibrium  characteristic  for  current  densities  below 
0.2  A/cm2.  In  this  current  density  range,  the  predicted  electron  temperature  remains  below  ap¬ 
proximately  8000  K  and  electron-impact  reactions  are  inefficient  in  ionizing  the  plasma.  Thus  the 
electron  number  density  increases  only  by  a  few  percent.  As  the  electron  temperature  increases,  the 
frequency  of  collisions  increases  with  resulting  in  a  decrease  of  the  electrical  conductivity  of 
the  plasma.  This  explains  why  the  E  vs.  j  characteristic  is  higher  than  the  resistive  equilibrium 
characteristic  for  j  below  ~0.2  A/cm2.  At  higher  values  of  the  current  density,  where  the  discharge 
produces  a  significant  increase  in  the  electron  density,  the  conductivity  increases  dramatically  and 
the  slope  of  the  E  vs  .j  characteristic  decreases.  Thus  the  region  to  the  right  of  the  resistive  equilib¬ 
rium  characteristic  is  where  the  discharge  increases  the  electron  number  density.  The  experimental 
data  at  Tg  between  1 800  K  and  2900  K  all  show  the  turning  trends  of  the  non-resistive  discharge 
characteristics.  We  note  however  that  the  predicted  resistive  part  appears  to  be  shifted  to  lower  cur¬ 
rent  densities  relative  to  the  experimental  curves.  This  difference  may  be  due  to  the  fact  that  the 
electron  number  density  was  slightly  above  the  equilibrium  value  in  the  incoming  air  stream.  We 
recall  that  the  “equilibrium”  air  was  produced  by  cooling  of  an  air  stream  initially  heated  to  high 
temperatures.  Slow  electron  recombination  could  therefore  explain  the  differences  at  low  current 
densities. 
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Current  Density,  j  (A/cm2) 


Figure  45.  Measured  (symbols)  and  predicted  (solid  lines)  electrical  discharge  characteristics  in  atmospheric 
pressure  air  plasmas  generated  by  DC  electric  discharges. 

Figure  45  also  shows  experimental  data  obtained  by  Stark  and  Schoenbach  [62]  in  an  at¬ 
mospheric  pressure  glow  discharge  in  air.  The  discharge  was  produced  between  a  microhollow 
cathode  and  a  positively  biased  electrode,  as  described  in  Ref.  [63].  The  gas  temperature  was 
measured  to  be  around  2000  K,  and  the  center  electron  number  density  is  reported  to  be  above  1012 
cm'3  [62,  64].  This  measurement  adds  further  support  to  the  kinetic  mechanism  predictions. 

We  conclude  this  section  with  the  power  required  to  produce  a  given  electron  density  in  air 
at  2000  K  by  means  of  DC  glow  discharges.  The  results  are  shown  in  Figure  46.  We  predict  that 
the  production  of  1013  electron/cm3  in  air  at  2000  K  requires  about  14kW/cm3.  The  corresponding 
electric  field  is  ~1.35kV/cm,  and  the  current  density  is  ~10.4A/cm2,  as  shown  earlier  in  Figure  29. 

This  level  of  Joule  heating  may  not  lead  to  significant  overall  gas  heating  in  small  scale  sta¬ 
tionary  DC  discharges  where  conduction  to  ambient  air  and  to  the  electrodes  is  high.  This  was  the 
case  for  instance  in  Gambling  and  Edels’s  experiments  [60]  where  the  positive  column  was  a  few 
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millimeters  in  length  and  approximately  0.2  mm2  in  area.  In  larger  volume  DC  discharges,  how¬ 
ever,  it  is  necessary  to  control  the  effect  of  Joule  heating  of  the  gas,  for  instance  by  flowing  the  gas 
through  the  discharge  at  high  velocities.  For  air  at  2000  K  flowing  through  a  1  cm  diameter  region 
of  length  3.5  cm  at  a  velocity  of  450  m/s,  the  residence  time  is  78  ps.  The  vibrational  relaxation 
times  x  reported  by  Park  [10]  indicate  that  the  fastest  vibrational  relaxation  rate  of  molecular  N2  is 
through  collisions  with  atomic  oxygen.  The  rate  constant  is  given  by  pot  =  10'6  (atm.s),  where  po  is 
the  partial  pressure  of  atomic  oxygen.  In  our  discharge  experiments,  the  atomic  oxygen  mole  frac¬ 
tion  is  less  than  1%,  according  to  our  two-temperature  kinetic  model  predictions.  Thus,  the  vibra¬ 
tional  relaxation  time  x  (>  100  ps)  is  larger  than  flow  time  (78  ps).  This  is  consistent  with  the  ob¬ 
servation  that  we  observed  little  gas  heating  in  the  experiments.  To  limit  gas  heating  to  acceptable 
levels  for  a  given  volumetric  power,  it  is  desirable  to  flow  the  plasma  at  high  velocity  through  the 
discharge  region. 


Figure  46.  Power  required  to  produce  an  elevated  electron  density  in  atmospheric  pressure  air  at  2000  K  by 

means  of  DC  discharges. 
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4.6.  DC  Discharges  in  Nitrogen 

Discharge  experiments  were  also  conducted  with  atmospheric  pressure  nitrogen  either  at 
room  temperature  or  preheated  to  temperatures  between  1800-2250  K  [57].  For  the  room  tempera¬ 
ture  experiments,  nitrogen  is  injected  between  the  electrodes  at  a  velocity  of  about  2  m/s.  The  gas 
heats  up  to  temperatures  between  2200  and  2800  K  as  a  result  of  discharge-induced  Joule  heating. 
For  the  experiments  with  preheated  nitrogen,  the  gas  is  injected  at  about  450  m/s.  Owing  to  the 
relatively  fast  flow  rate,  the  gas  temperature  remains  practically  constant  in  the  discharge  region. 
The  discharge  diameter  is  about  1.7  mm  in  the  room  temperature  nitrogen  experiments,  and  5  mm 
in  preheated  nitrogen.  The  difference  in  discharge  diameters  is  due  to  the  larger  radial  thermal  gra¬ 
dients  in  cold  nitrogen  than  in  the  preheated  flow.  In  all  cases  the  measured  cathode  fall  is  about 
300  V ,  which  is  typical  of  a  glow  discharge.  The  measured  discharge  characteristics  are  shown  in 
Figure  47,  along  with  the  characteristics  calculated  with  the  nitrogen  CR  model.  As  shown  in 
Figure  47,  the  measured  characteristics  support  the  results  of  the  chemistry  and  discharge  models. 

The  foregoing  calculations  were  performed  using  Maxwellian  distribution  functions  for  the 
translational  energies  of  the  electrons  and  heavy  particles.  To  investigate  the  limitations  of  this  sim¬ 
plification,  the  electron  Boltzmann  equation  was  solved  for  nitrogen  plasmas  in  conjunction  with  a 
somewhat  simplified  form  of  the  collisional-radiative  model  [57],  The  results,  shown  in  Figure  48, 
do  not  differ  significantly  from  those  with  the  Maxwellian  version. 


Current  density,  j  (A/cm2) 

Figure  47.  Measured  (symbols)  and  predicted  DC  discharge  characteristics  in  nitrogen  at  1  atm. 
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j  (A/cm2) 

Figure  48.  Predicted  discharge  characteristics  in  atmospheric  pressure  nitrogen  at  2000  K. 


4.7.  Conclusions 

Investigations  have  been  made  of  the  mechanisms  of  ionization  in  two-temperature  air 
plasmas  with  electron  temperatures  elevated  with  respect  to  the  gas  temperature.  Numerical  simu¬ 
lations  of  these  mechanisms  yield  the  notable  result  that  the  electron  number  density  exhibits  an  S- 
shape  dependence  on  the  electron  temperature  at  fixed  gas  temperature.  This  S-shaped  behavior  is 
caused  by  competing  ionization  and  charge  transfer  reactions.  The  characteristic  of  electric  field 
vs.  current  density  also  exhibits  a  non-monotonic  dependence. 

Discharge  experiments  were  conducted  in  air  at  atmospheric  pressure  and  temperatures 
ranging  from  1 800  to  3000  K.  In  these  experiments,  a  DC  electric  field  was  applied  to  flowing  air 
plasmas  with  electron  concentrations  initially  close  to  equilibrium.  These  experiments  have  shown 
that  it  is  possible  to  obtain  stable  diffuse  glow  discharges  in  atmospheric  pressure  air  with  electron 
number  densities  of  up  to  2.5x1 012  cm'3,  which  is  up  to  six  orders  of  magnitude  higher  than  in  the 
absence  of  the  discharge  The  value  of  2.5xl012  cm'3  corresponds  to  the  maximum  current  that  can 
be  drawn  from  the  250  mA  power  supply  used  in  these  experiments.  The  diffuse  discharges  are  ap¬ 
proximately  3.5  cm  in  length  and  3.2  mm  in  diameter.  No  significant  degree  of  gas  heating  was  no¬ 
ticed  as  the  measured  gas  temperature  remained  close  to  its  value  without  the  discharge  applied. 
Results  from  these  experiments  are  in  good  agreement  with  the  predicted  E  vs.  j  characteristics. 
Additional  comparisons  were  made  with  results  from  glow  discharge  experiments  in  atmospheric 
pressure  ambient  air  by  Gambling  and  Edels  [60]  and  Stark  and  Schoenbach  [62].  The  measure¬ 
ments  of  these  authors  are  also  consistent  with  the  predicted  E  vs  .j  characteristics.  As  these  meas- 
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urements  were  made  in  the  reactive  region  of  the  E  vs.  j  curve,  they  support  our  proposed  mecha¬ 
nism  of  ionization  for  two-temperature  air. 

As  the  power  budget  for  DC  electron  heating  is  higher  than  desired  for  the  practical  use  of 
air  plasmas  in  many  applications,  methods  to  reduce  the  power  budget  were  explored  in  our  labora¬ 
tory.  Based  on  the  predictions  of  our  chemical  kinetics  and  electrical  discharge  models,  we  have 
found  that  a  repetitively  pulsed  electron  heating  strategy  can  provide  power  budget  reductions  of 
several  orders  of  magnitude  with  respect  to  DC  electron  heating.  Repetitively  pulsed  discharges  are 
the  object  of  Section  5  of  this  report. 
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5.  Repetitively  Pulsed  Discharges  in  Air 


5.1.  Introduction 

As  the  power  required  to  sustain  elevated  electron  densities  with  DC  discharges  is  large,  a 
power  reduction  strategy  based  on  pulsed  electron  heating,  proposed  in  our  1999  MURI  Annual 
Progress  Report,  was  explored.  This  strategy  is  illustrated  in  Figure  49.  Short  voltage  pulses  are 
applied  to  increase  the  electron  number  density.  After  each  pulse,  ne  decreases  according  to  electron 
recombination  processes.  When  ne  reaches  the  minimum  desired  value,  a  second  pulse  is  applied. 
The  average  electron  density  obtained  with  this  method  depends  on  the  pulse  duration,  pulse  volt¬ 
age,  and  the  interval  between  pulses.  Since  power  is  not  expended  during  the  recombination  period, 
the  time  averaged  power  will  be  greatly  reduced  if  the  pulse  time  is  much  smaller  than  the  recombi¬ 
nation  time. 


time 

Figure  49.  Repetitively  pulsed  strategy. 


As  seen  in  Section  4,  DC  discharges  can  maintain  ne  >  1012  cm'3  in  atmospheric  pressure  air  with 
electric  fields  producing  an  electron  temperature  on  the  order  of  leV.  To  produce  the  same  average 
electron  density  with  short  (l-lOns)-pulsed  discharges,  a  higher  electron  temperature  of  about  3-5 
eV  is  required.  Although  the  corresponding  field  is  higher  than  for  a  DC  discharge,  the  ionization 
efficiency  is  larger  in  the  pulsed  case  than  in  the  DC  case  because  the  energy  lost  to  nitrogen  mole¬ 
cules,  per  electron  created,  is  smaller  at  Te=3-5  eV  than  at  1  eV.  This  coupled  with  the  finite  elec¬ 
tron  recombination  time  allows  the  power  budget  to  be  dramatically  reduced  with  pulsed  discharges. 


66 


Figure  50  shows  the  predicted  inelastic  energy  losses  of  electrons  by  collisions  with  N2,  per  unit 
number  density  of  N2  and  electrons.  The  losses  to  nitrogen  represent  the  main  fraction  of  the  losses 
in  air.  This  is  because  at  electron  temperatures  below  ~20,000  K  the  resonant  e-V  transfer  in 
ground  state  N2  is  by  far  the  dominant  loss  channel  (there  is  no  such  resonant  channel  for  02  or 
NO).  At  electron  temperatures  above  30,000  K,  the  inelastic  losses  are  dominated  by  electron- 
impact  electronic  excitation,  dissociation,  and  ionization,  and  the  total  losses  per  unit  number  den¬ 
sity  of  N2  and  electrons  are  about  the  same  as  the  total  losses  per  unit  number  density  of  02  and 
electrons.  Nitrogen  losses  dominate  at  Te  >  20,000  K  because  the  density  of  N2  is  much  higher  than 
the  density  of  02.  Figure  50  shows  that  the  power  into  ionization  represents  an  increasingly  large 
fraction  of  the  total  power  as  the  electron  temperature  is  increased. 


Figure  50.  Inelastic  power  losses  by  electron-impact  vibrational  excitation,  electronic  excitation,  dissociation,  and 
ionization  of  N2.  In  these  calculations,  the  vibrational  temperature  is  fixed  equal  to  the  gas  temperature  (Tv=Tg 
=2000  K),  and  the  electronic  temperature  of  internal  energy  levels  is  fixed  equal  to  the  electron  temperature. 
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Figure  51.  Schematic  of  the  combined  pulsed  and  DC  discharge  experiments. 

5.2.  Experiments  with  a  Single  Pulse 

To  test  the  pulsing  scheme,  experiments  were  conducted  [65]  in  atmospheric  pressure,  2000 
K  air  using  a  pulse  forming  line  capable  of  generating  a  10  ns  rectangular  pulse  with  peak  voltage 
up  to  16  kV.  To  experimentally  simulate  the  conditions  of  a  repetitively  pulsed  discharge,  the  initial 
elevated  electron  number  density  generated  by  the  “previous”  pulse  is  created  by  means  of  a  DC 
discharge  in  parallel  with  the  pulser.  The  circuit  schematic  is  shown  in  Figure  5 1 .  With  a  DC  volt¬ 
age  of  2  kV  and  current  of  150  mA,  the  initial  electron  density  is  6.5x10“  cm'3.  A  10  kV,  10  ns 
pulse  is  superimposed  to  further  increase  the  electron  density.  The  measured  discharge  diameter  of 
about  3  mm  is  comparable  with  the  diameter  of  the  DC  discharge  (Figure  52).  The  temporal  varia¬ 
tion  of  plasma  conductivity  was  measured  from  the  voltage  across  the  electrodes  and  the  current 
density  through  the  plasma.  The  electron  density  increases  from  6.5x10“  to  9x1 012  cm'3  during  the 
pulse,  then  decays  to  1012  cm'3  in  about  12  |xs  (Figure  53).  The  average  measured  electron  density 
over  the  12  |J.s  duration  is  2.8x1 012  cm'3. 

Figure  53  shows  a  comparison  of  the  measured  electron  number  density  with  the  predictions 
of  our  two-temperature  model.  The  predictions  agree  well  with  the  measured  electron  decay  time  of 
12  (is.  This  decay  time  is  consistent  with  the  dissociative  recombination  time  of  NO+ predicted  to 
be  8.7  jus  without  the  DC  background.  Thus  these  results  provide  validation  of  our  chemical  kinetic 
model  of  the  recombination  phase. 


68 


Radius  (mm) 

Figure  52.  Spatial  extent  of  the  plasma  produced  with  pulsed  and  DC  discharges. 


Time  (ns) 

Figure  53.  Temporal  electron  density  profile  in  the  10  ns  pulsed  discharge 


5.3.  Experiments  with  100  kHz  Repetitive  Discharge 

The  success  of  the  proof-of-concept  experiments  conducted  with  the  single  pulse  discharge 
led  us  to  investigate  the  generation  of  air  plasmas  with  a  repetitively  pulsed  discharge.  A  repetitive 
pulser  capable  of  generating  10  ns  pulses,  with  peak  voltages  of  3-12  kV  and  pulse  repetition  fre¬ 
quencies  up  to  100  kHz  was  acquired  from  Moose-Hill/FID  Technologies.  This  pulser  operates  with 
a  solid-state  opening  switch  or  Drift-Step  Recovery  Diode  (DSRD).  The  experimental  set-up  is 
shown  in  Figure  54  and  the  electric  circuit  in  Figure  55.  The  discharge  is  applied  to  preheated,  LTE 
air  at  atmospheric  pressure  and  about  2000  K.  The  DC  circuit  in  parallel  with  the  pulser  was  used 
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only  to  determine  the  electron  number  density  from  the  plasma  conductivity.  In  regular  operation, 
the  DC  circuit  is  disconnected  and  the  discharge  operated  with  the  pulser  only. 


Figure  55.  Repetitive  pulse  discharge  circuit  schematic 
(DC  circuit  applied  only  for  conductivity  measurements). 


Figure  56.  Photograph  of  10  ns,  100  kHz  repetitive  pulse  discharge  in  air  at  2000  K,  1  atm. 


Figure  57.  Time-resolved  images  of  10  ns  pulsed  discharge  in  air  at  2000  K,  1  atm. 


71 


A  photograph  of  the  repetitively  pulsed  discharge  in  operation  in  atmospheric  pressure  pre¬ 
heated  (2000  K)  air  is  shown  in  Figure  56.  The  diffuse  character  of  the  discharge  was  confirmed 
with  time-resolved  (1.5  ns  frames  every  2  ns)  measurements  of  plasma  emission  during  the  pulse 
(see  Figure  57).  These  measurements  were  made  with  a  high-speed  intensified  camera,  Roper  Sci¬ 
entific  PI-MAX 1024.  The  diameter  of  the  discharge  is  approximately  3.3  mm.  Additional  time- 
and  spectrally-resolved  measurements  of  emission  during  the  pulse  and  the  recombination  phase 
show  that  the  pulse  excites  the  C  state  of  N2  and  the  A  state  of  NO.  After  the  pulse,  emission  from 
the  C  state  of  N2  decays  to  a  constant  value  within  30  ns,  and  emission  from  the  A  state  of  NO 
shows  a  two-step  decay,  with  first  an  abrupt  decrease  by  over  four  orders  of  magnitude  from  the 
end  of  the  pulse  until  320  ns  after  the  pulse,  and  then  a  slower  decrease  by  one  order  of  magnitude 
until  the  next  pulse. 

Figure  58  shows  the  measured  temporal  variations  of  the  electron  density  during  three  cy¬ 
cles  of  the  pulsed  discharge.  The  electron  number  density  varies  from  7x1 011  to  1.7xl012  cm'3,  with 
an  average  value  of  about  10 12  cm"3.  The  power  deposited  into  the  plasma  by  the  repetitive  dis¬ 
charge  was  determined  from  the  pulse  current  (measured  with  a  Rogowski  coil),  the  voltage  be¬ 
tween  the  electrodes  (6  kY  peak)  minus  the  cathode  fall  voltage  (measured  to  be  1 525  V  by  varying 
the  gap  distance),  and  the  measured  discharge  diameter.  The  peak  pulse  current  was  240  mA.  The 
power  deposited  is  found  to  be  12  W/cm3,  consistent  with  the  theoretical  value  of  9  W/cm3  for  an 
optimized  pulsed  discharge  producing  1012  electrons/cm3.  It  is  lower,  by  a  factor  of  250,  than  the 
power  of  3000  W/cm3  required  to  sustain  1012  electrons/cm3  with  a  DC  discharge. 

More  details  about  these  experiments  and  modeling  can  be  found  in  Ref.  [66],  In  this  refer¬ 
ence,  a  study  was  made  of  the  effect  of  the  pulse  repetition  frequency.  Experiments  were  conducted 
with  repetition  frequencies  of  30  and  100  kHz.  In  both  case  the  power  requirements  were  close  to 
10  W/cm3  for  about  1012  electrons/cm3.  The  main  difference  between  the  plasmas  produced  is  the 
amplitude  of  electron  density  variations.  In  the  30  kHz  discharge,  the  amplitude  varies  by  about  a 
factor  of  10,  whereas  in  the  100  kHz  the  amplitude  varies  by  a  factor  of  two  only. 
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Figure  58.  Electron  number  density  measurements  in  the  repetitive  pulse  discharge  in  air  at  1  atm,  2000  K. 


The  results  of  our  research  on  DC  and  pulsed  electrical  discharges  are  summarized  in  Figure 
59,  which  shows  the  power  required  to  generate  elevated  electron  number  density  in  2000  K,  at¬ 
mospheric  pressure  air,  with  DC  and  pulsed  discharges.  The  experimental  point  represents  the 
measured  power  requirement  of  our  repetitively  pulsed  discharge  experiment.  Power  budget  reduc¬ 
tions  by  an  additional  factor  of  about  5  are  possible  with  repetitive  pulses  of  1  nanosecond  duration. 
Such  repetitive  pulsers  are  already  commercially  available.  Therefore,  power  budget  reductions  by  a 
factor  of  1000  relative  to  the  DC  case  at  1012  electrons/cm3  can  be  readily  obtained  with  a  repeti¬ 
tively  pulsed  technique. 
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Figure  59.  Power  budget  requirements  vs.  electron  number  density  for  DC  and  pulsed  discharges 

in  air  at  1  atm,  2000  K. 
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6.  Computational  Analysis  of  Diffuse  Discharges  in  Atmospheric  Pressure  Air 
6.1.  Introduction 

In  this  part  of  the  report,  we  use  computational  simulations  to  study  the  mechanisms  that  de¬ 
termine  the  size  and  stability  of  the  DC  discharges  described  in  Section  4  of  this  report.  The  CFD 
models  the  gas  flow  coupled  to  the  Joule  heating  of  the  electrons.  A  three-temperature  model  for 
the  plasma  energetics  is  used  in  conjunction  a  two-temperature  chemical  kinetics  model  for  11- 
species  air.  Finite-rate  coupling  between  the  energy  modes  and  all  of  the  relevant  transport  proc¬ 
esses  are  included.  We  use  this  code  to  simulate  the  Stanford  University  DC  and  pulsed  plasma 
torch  experiments.  We  find  that  turbulent  mixing  in  the  torch  jet  appears  to  affect  the  discharge 
behavior,  and  a  simplified  flow  model  produces  results  consistent  with  those  obtained  in  experi¬ 
ments. 


6.2.  Flow  Field  and  Discharge  Model 

The  flow  within  the  plasma  torch  experiment  test  section  is  described  by  the  Navier-Stokes 
equations  that  have  been  extended  to  include  the  effects  of  nonequilibrium  thermochemistry  and  the 
applied  electric  field.  Thus,  we  solve  a  separate  mass  conservation  equation  for  each  of  the  chemi¬ 
cal  species  to  allow  the  finite-rate  formation  and  destruction  of  all  relevant  chemical  species.  In  our 
simulations  of  the  air  plasmas,  we  use  an  11 -species  model  that  accounts  for  all  air  species  as  well 
as  their  ions  (N2,  O2,  NO,  N,  O,  N2+,  02+,  NO+,  N+,  0+,  e).  The  chemical  kinetics  of  the  multi¬ 
temperature  plasma  are  described  by  the  model  developed  at  Stanford  under  the  APRP.  Additional 
details  of  the  model  are  given  in  this  report  and  in  [65]  and  [13]. 

The  velocity  field  is  described  by  solving  conservation  equations  for  the  radial,  axial,  and 
swirl  direction  momenta.  Under  typical  torch  operating  conditions,  about  half  of  the  test  gas  is  in¬ 
jected  into  the  torch  with  a  component  of  velocity  in  the  swirl  direction,  resulting  in  a  medium  to 
high  degree  of  swirl  at  the  entrance  to  the  test-section.  The  swirl  was  thought  to  be  instrumental  for 
explaining  the  poor  agreement  of  the  previous  simulations  with  experimental  results  as  far  as  the 
radial  profiles  of  individual  variables  are  concerned.  The  effect  of  swirl  is  properly  accounted  for 
in  the  present  simulations. 

When  an  electric  field  is  applied  to  the  plasma,  the  rate  of  energy  transfer  between  the  dif¬ 
ferent  energy  modes  largely  determines  the  steady-state  electron  concentration  in  the  plasma.  Thus, 
a  complete  model  of  the  plasma  energetics  is  required.  We  solve  a  total  energy  conservation  equa¬ 
tion,  a  separate  mass-averaged  vibration-electronic  energy  conservation  equation,  and  an  electron 
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translational  energy  equation.  The  vibration-electronic  energy  of  a  given  species  is  defined  to  be 
the  difference  between  that  species'  internal  energy  computed  from  the  Gordon  and  McBride  data 
[67]  and  the  sum  of  its  translational-rotational  energy  and  heat  of  formation.  This  approach  is  taken 
by  Gnoffo  [68],  The  rate  of  electrical  energy  input  and  the  rate  of  energy  transfer  between  the 
modes  are  modeled  as  detailed  below. 

In  our  model,  the  DC  discharge  adds  energy  to  the  plasma  through  the  electrons  via  the 
Joule  heating  term.  As  a  result,  the  electron  temperature  increases,  causing  electron-impact  ioniza¬ 
tion  and  the  further  production  of  electrons,  thus  cooling  the  electron  gas.  At  the  same  time,  the 
electrons  transfer  their  energy  to  the  heavy  particles  through  collisions  and  inelastic  losses.  We 
represent  model  these  processes  with  the  following  energy  conservation  equations.  In  general,  the 
change  of  energy  =  divergence  of  energy  flow  +  energy  source  terms. 

j-2 

V  •(( E  +  p)u  +  q+u-z)+  — 

a 

-V  •  ( Evu  +  qv)+QT_v  +  Qv.e  +  £  wsevs 

s 

-V  •  (Ee  u  +  qe )+ i- —  Qh_e  -  Qv_e  +  weee 
a 

Note  that  the  total  energy,  E,  has  a  source  term  due  to  the  Joule  heating,  as  well  as  the  usual 
viscous  transport  and  pressure-work  terms.  The  vibrational  energy,  Ev,  has  the  source  terms  repre¬ 
senting  energy  transfer  to  vibration  from  the  V-T  relaxation  and  the  electron-vibration  energy  trans¬ 
fer.  In  addition,  the  last  term  on  the  right-hand  side  represents  the  addition  and  loss  of  vibrational 
energy  due  to  chemical  reactions.  The  electron  energy  equation  has  the  same  form,  but  also  in¬ 
cludes  the  Joule  heating  term.  The  heat  flow  vectors  q ,  qv ,  and  qe ,  contain  heat  conduction  and 
heat  carried  by  mass  diffusion.  The  semi-empirical  formula  of  Blottner  et  al.  [69]  and  curve  fits 
provided  by  Gupta  et  al.  [70]  are  used  to  obtain  the  energy  transport  properties.  Unfortunately,  for 
flows  with  bulk  phenomena  such  as  turbulent  mixing,  a  general  multi-component  approach  for 
transport  properties  is  not  available.  This  limitation  will  be  discussed  later  in  the  report. 

The  energy  transfer  terms  Qr-V,  Qh-e,  and  Qv.e,  the  translation-vibration  energy  transfer  rate, 
the  heavy-particle-electron  energy  transfer  rate,  and  the  vibration-electron  energy  transfer  rate,  re¬ 
spectively,  are  taken  from  this  report,  as  well  as  [13]  and  [71].  The  diffusional  mass  fluxes  carry 
energy  and  are  calculated  per  species  using  the  explicit  expression  derived  by  Ramshaw  [72]  in  his 
self-consistent  effective  binary  diffusion  approximation.  The  driving  forces  in  this  formulation  are 
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dt 

dt 
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the  species  concentration  gradients  and  the  body  forces  acting  on  the  charged  species  due  to  the 
ambipolar  electric  field.  Ambipolar  diffusion  is  included  because  concentration  gradients  of 
charged  particles  exist  in  regions  of  no  current  flow  and  their  diffusion  is  curbed  by  the  presence  of 
neighboring  charged  particles.  The  first  expansion  of  the  binary  diffusion  coefficient,  utilizing  the 
weighted  averages  of  the  pair  cross  sections,  is  used. 

The  discharge  is  assumed  to  be  axisymmetric,  with  no  Lorentz  force  and  a  coincident  axial 
electric  current.  As  shown  in  the  energy  conservation  equations,  the  electric  field  influences  the 
total  energy  and  the  electron  energy  through  the  Joule  heating  source  term.  The  current  carried  by 
the  plasma  is  controlled  experimentally,  and  is  thus  a  prescribed  parameter.  Therefore,  the  current 
must  be  conserved  at  each  axial  location  between  the  electrodes.  A  flux  of  electrons  is  introduced 
at  the  cathode  location.  It  is  calculated  from  the  total  current  in  the  discharge  and  the  effective  area 
of  the  cathode.  Sheath  details  are  not  resolved,  and  the  current  is  considered  to  be  entirely  carried 
by  the  electrons.  Between  the  electrodes  the  charged  particles  are  acted  upon  by  the  electric  field. 
The  flux  of  electrons  passing  through  a  unit  area  per  unit  time  is 

re  =  PX 

where  ue  is  the  mean  drift  velocity  and  is  related  to  the  kinetic  energy  acquired  by  the  electrons  un¬ 
der  the  influence  of  the  electric  field.  The  mean  drift  velocity  of  the  electrons  depends  on  the  elec¬ 
tric  field  and  the  electron  mobility,  /ue 

The  electron  mobility  is  related  to  the  collision  frequency,  ve,  between  electrons  and  heavy 
particles. 

e 


For  the  experimental  conditions,  the  ions  have  negligible  concentrations  and  the  collision 
frequency  is  approximated  by 

K=nV,heQen 

where  n  is  the  number  density  and  v,/,e  is  the  thermal  velocity  of  the  electrons, 


and  the  average  cross-section  Qen  is  the  average  momentum  transfer  cross  section  taken  to  be 
10'15  cm"3. 

The  most  critical  issue  of  the  discharge  modeling  is  how  to  enforce  current  conservation  be¬ 
tween  the  electrodes.  It  is  possible  to  solve  some  form  of  a  Poisson  equation  for  the  electrical  po¬ 
tential  as  a  function  of  the  local  plasma  properties.  However,  this  approach  has  a  number  of  draw- 
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backs,  including  many  numerical  difficulties,  expense  of  solving  the  Poisson  equation,  and  diffi¬ 
culty  with  enforcing  total  current  conservation  at  each  axial  location  between  the  electrodes.  Also, 
because  the  region  of  the  discharge  is  well  defined,  a  simpler  more  direct  approach  may  be  used. 
The  basic  concept  is  to  directly  enforce  total  current  conservation  at  each  axial  location  between  the 
electrodes. 

To  formulate  the  discharge  model,  we  start  with  the  current  density,/,  from 

j  ~  }e  =  ~ene% 

where  again  ue  is  the  mean  drift  velocity.  If  we  assume  that  there  is  no  variation  in  the  drift  veloc¬ 
ity  in  the  radial  direction,  we  can  then  write  the  total  current,  /,  as 

/  =  J  jdA  =  -eue  (x)  jnedA 

A  A 

where  the  integral  is  taken  in  the  radial  direction,  and  again  we  assume  that  the  mean  drift  velocity 
is  only  a  function  of  the  axial  distance,  x.  We  can  then  compute  ue(x )  from  the  above  equation. 
We  then  obtain  the  local  value  of  j  for  use  in  the  Joule  heating  source  term,  f/a.  The  electrical 
conductivity  is  defined  using  the  calculated  particle  mobility  and  electron  number  density 

<7  =  fttent 

This  approach  is  essentially  equivalent  to  the  classical  channel  model  discussed  in  Ref.  [73]. 
The  radial  extent  of  the  discharge  is  then  determined  by  the  local  electron  number  density  in  rela¬ 
tion  to  the  total  electron  number  density  at  each  axial  location.  Thus,  the  rate  of  ionization,  energy 
transfer  between  modes,  and  the  rate  of  radial  diffusion  of  the  electrons  reach  a  balance  to  define 
the  extent  of  the  elevated  electron  concentration  region. 

This  discharge  model  will  be  accurate  in  regions  where  the  equipotential  lines  are  not 
curved,  which  is  a  good  approximation  for  most  of  the  discharge.  However,  near  the  cathode  the 
discharge  radius  visually  appears  to  expand  rapidly,  and  the  potential  lines  tend  to  wrap  around  the 
cathode.  But  this  is  a  localized  effect  (confined  to  about  a  discharge  diameter)  and  does  not  appear 
to  have  an  important  influence  on  the  predicted  discharge  characteristics. 


6.3.  Numerical  Method 

The  conservation  equations  discussed  above  and  discussed  in  more  detail  in  [65]  and  [74], 
may  be  written  as 

dU  dF  IdrG  _ 

- + — + - =  W 

dt  dx  r  dr 

where  U  is  the  vector  of  conserved  variables,  F  and  G  are  the  flux  vectors,  and  W  is  the  source  vec¬ 
tor.  The  equations  may  be  solved  using  many  different  computational  fluid  dynamics  methods.  We 
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choose  to  use  a  finite-volume  formulation  and  evaluate  the  fluxes  with  an  upwind  modified  Steg- 
ger-Warming  flux  vector  splitting  approach  [75]. 

Under  the  conditions  of  the  DC  discharge  experiments,  the  energy  relaxation  processes  are 
very  fast  relative  to  the  fluid  motion  time  scales  and  the  chemical  kinetic  processes.  To  handle  this 
large  disparity  in  characteristic  time  scales,  we  would  usually  use  an  implicit  time  integration 
method  [75].  However,  for  this  problem  a  complete  linearization  of  the  problem  is  itself  very  ex¬ 
pensive.  We  solve  17  conservation  equations,  and  the  cost  of  evaluating  the  Jacobians  and  invert¬ 
ing  the  system  scales  with  the  square  of  the  number  of  equations.  Therefore,  we  linearize  only 
those  terms  that  are  relatively  fast,  which  results  in  a  simple  and  inexpensive  semi-implicit  method 
that  very  substantially  reduces  the  cost  of  the  calculations. 


The  relatively  fast  terms  are  the  internal  energy  relaxation  and  the  Joule  heating  terms  in  the 
source  terms  for  the  three  energy  equations.  Therefore,  we  split  the  source  vector,  W,  into  those 
“fast”  terms,  Wfasl,  and  all  of  the  other  terms,  Wsi0w-  The  conservation  equations  are  then  written  as 


dU  dF  1  drG 

- I - 1 - : 

dt  dx  r  dr 


■wfasl+wsl0W 


We  then  linearize  W/ast  in  time 


w;::=w;s,+c;astsun+o(At2) 

where  Cfast  is  the  Jacobian  of  Wfas,  with  respect  to  U  and  SU"  =  Un+'  -  U” .  Because  of  the  form  of 
Wfas,,  Cfast  is  a  simple  matrix  that  can  be  inverted  analytically.  Then  the  solution  is  integrated  in 
time  using 

.  f  f  r)P"  1  r)rCn 

U”+'  =(i-Atc;asiy  At(w^,+w;;0j-At 


This  approach  increases  the  stable  time  step  by  a  factor  of  at  least  50  compared  to  an  ex¬ 
plicit  Euler  method  at  essentially  no  additional  computational  time  per  time  step.  This  results  in  a 
very  large  reduction  in  the  computer  time  required  to  obtain  a  steady-state  solution  and  enables  the 
simulation  of  the  Stanford  plasma  experiments.  This  numerical  method  has  also  been  used  to  simu¬ 
late  the  pulsed  discharge  experiments  performed  at  Stanford  University  [65].  We  find  that  the  time 
step  must  be  reduced  during  the  duration  of  the  pulse,  but  then  the  time  step  may  be  increased 
significantly  as  the  plasma  cools  between  pulses. 


A  two-block  grid  is  used  to  facilitate  the  implementation  of  the  boundary  conditions.  The 
first  grid  block  represents  a  portion  of  the  plasma  torch,  the  nozzle  and  the  test  section,  and  the  sec¬ 
ond  grid  block  acts  as  a  large  constant-pressure  exhaust  reservoir  at  one  atmosphere.  A  schematic 
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of  the  computational  grid  is  shown  in  Figure  60.  Grid  block  1  has  35  points  in  the  axial  direction 
and  35  points  in  the  radial  direction,  with  exponential  clustering  near  the  test  section  surface.  The 
second  block  is  110x75,  with  an  exact  matching  of  the  radial  spacing  at  the  block  interface. 


Figure  60.  Schematic  of  two-block  grid  used  to  model  the  Stanford  discharge  experiments. 

The  inflow  boundary  conditions  are  set  by  varying  the  assumed  inflow  static  pressure  until 
we  obtain  the  experimental  gas  mass  flow  rate  of  4.9  g/s.  The  inflow  is  assumed  to  be  in  LTE  at  the 
measured  temperature  profile.  This  results  in  a  consistent  representation  of  the  inflow  conditions. 
The  boundary  conditions  along  the  test-section  surface  are  straight-forward.  The  velocity  is  zero  at 
the  surface,  the  temperature  is  specified,  and  the  normal-direction  pressure  gradient  is  zero.  We  as¬ 
sume  that  the  metallic  surface  is  highly  catalytic  to  ion  recombination.  Otherwise,  the  surface  is  as¬ 
sumed  to  be  non-catalytic  to  recombination  for  neutrals. 

The  computation  is  initialized  as  follows:  first,  the  inflow  conditions  are  specified  as  above. 
Then  the  test-section  and  reservoir  are  all  initialized  at  atmospheric  pressure,  and  at  each  axial  loca¬ 
tion  the  temperature  profiles  and  chemical  concentration  profiles  are  set  identical  to  the  inflow 
boundary  profiles.  Once  a  converged  solution  is  obtained  for  the  flow  in  LTE,  the  discharge  is  ig¬ 
nited  by  injecting  a  flux  of  electrons  at  the  cathode  and  applying  the  Joule  heating  source  term  to  the 
energy  equations.  Then  a  steady-state  solution  for  the  DC  discharge  is  obtained. 
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6.4.  Results 

Let  us  first  consider  the  flow  of  the  torch  jet  into  the  atmosphere  without  electrical  energy 
addition.  Figure  61  plots  contours  of  velocity  magnitude  and  density,  overlaid  with  instantaneous 
streamlines.  Clearly  the  jet  is  unstable,  with  the  shear  layer  rolling  up  into  vortices.  This  results  in 
an  unsteady  flow  field,  with  the  vortical  structures  propagating  downstream  at  approximately  the 
shear  layer  mean  convective  speed.  Simulations  have  been  carried  out  to  many  characteristic  times, 
and  these  structures  persist  in  the  calculations. 

Previous  results  [65,  74,  76]  did  not  demonstrate  this  unsteadiness,  but  those  calculations 
used  an  artificially  high  reservoir  temperature  of  800  K,  rather  than  the  300  K  used  here.  Thus,  the 
present  results  have  a  reservoir  density  that  is  nearly  three  times  larger  and  a  significantly  lower 
viscosity  in  the  shear  layer.  These  differences  greatly  enhance  the  jet  instability.  Additionally,  the 
roll-up  of  vortical  structures  in  a  shear  layer  is  greatly  affected  by  the  dissipation  of  the  numerical 
method.  In  our  results,  we  have  taken  care  to  reduce  the  numerical  error  by  using  a  second-order 
accurate  with  low  dissipation  and  by  clustering  the  grid  at  the  jet  interface.  These  high-quality  nu¬ 
merics  increases  the  propensity  of  the  jet  to  form  vortical  structures. 

The  flow  field  described  by  the  Navier-Stokes  equations  assumes  that  the  flow  from  the 
nozzle  to  the  open  atmosphere  is  laminar.  The  flow  at  the  end  of  the  nozzle  section  has  a  Reynolds 
number  of  about  104,  suggesting  that  the  flow  in  the  nozzle  is  likely  transitional  or  turbulent.  The 
cold  gas  mixing  upstream  of  the  nozzle  also  enhances  the  turbulent  motion  in  the  nozzle.  Thus  it  is 
highly  likely  that  the  jet  flow  is  transitional  or  turbulent.  The  effect  of  turbulence  is  to  increase  the 
mixing  across  the  jet  shear  layer,  which  tends  to  enhance  its  stability  and  reduce  the  formation  of 
the  organized  vortical  structures  seen  in  Figure  61. 

With  the  background  flow  described  above,  we  now  include  electrical  energy  addition.  For 
computational  reasons,  we  establish  the  discharge  by  propagating  the  electron  flux  from  the  cathode 
downstream  to  the  anode  before  “igniting”  the  discharge.  This  provides  an  initial  channel  of  elec¬ 
trons  to  carry  the  current.  These  electrons  absorb  the  electrical  energy,  are  boosted  to  high  tem¬ 
perature,  and  rapidly  produce  more  electrons  until  a  balance  between  the  energy  addition  and  losses 
is  established. 

Once  the  discharge  is  ignited,  it  is  impossible  to  obtain  a  steady-state  solution  for  the  dis¬ 
charge.  Figure  62  plots  contours  of  vibrational  temperature  and  electron  number  density  6.4  psec 
after  the  energy  addition  is  started.  Note  the  non-uniformity  in  the  discharge,  with  substantial 
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tions  in  ne  along  the  axis.  In  particular,  the  elevated  region  of  ne  corresponds  to  the  narrowing  of  the 
jet  by  the  vortical  structures.  The  discharge  is  not  robust,  and  as  the  cycle  of  vortical  motion  occurs 
in  the  shear  layer,  the  discharge  is  pinched  and  transforms  into  an  arc. 

It  should  be  noted  that  the  time  scale  of  the  vortical  motion  is  much  larger  than  that  of  the 
discharge  formation.  The  vortices  are  convected  downstream  at  the  mean  convective  speed,  which 
for  this  case  is  about  200  m/s.  Thus,  the  time  required  for  a  vortex  to  propagate  between  the  elec¬ 
trodes  is  about  200  (is.  This  is  obviously  much  larger  than  the  time  required  for  the  discharge  to 
collapse  into  an  arc.  Therefore,  it  is  highly  unlikely  that  the  modifications  to  the  flow  field  by  the 
electrical  energy  addition  would  prevent  arc  formation. 

In  spite  of  these  numerical  results,  the  experimentally  observed  discharge  is  stable,  and  the 
obvious  question  is:  Why  does  arcing  happen  in  the  computations,  but  not  in  the  experiments? 


Figure  61.  Contours  of  velocity  magnitude  and  density  at  the  experimental  conditions  without  the  discharge. 
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Figure  62.  Contours  of  vibrational  temperature  and  electron  number  density,  6.4  fisec  after  discharge  initiation. 

One  important  mechanism  to  stabilize  the  discharge  is  radial  heat  conduction.  If  the  flow  is 
transitional  or  turbulent,  the  heat  flux  must  have  a  component  due  to  turbulent  motion.  There  is  no 
adequate  model  for  predicting  turbulence  in  gas  discharges  or  plasmas.  Huang  [77]  introduced  a 
two-fluid  turbulence  model  for  thermal  plasma  sprays  that  qualitatively  reproduces  the  large-scale 
transport.  He  showed  that  pockets  of  heated  gas  move  away  from  the  shear  layer,  convectively 
transporting  large  amounts  of  heat.  This  mechanism  would  also  inhibit  the  formation  of  large-scale 
vortical  structures. 

Unfortunately,  without  additional  information  about  the  properties  of  the  experimental  jet,  it 
is  very  difficult  to  make  predictive  simulations  for  turbulent  flows.  We  could  include  a  turbulence 
model,  but  this  would  involve  a  tedious  and  unrewarding  game  of  parameter  adjustment  (tuning)  to 
get  reasonable  results.  We  have  therefore  resorted  to  an  approximate  approach  to  better  understand 
the  role  of  mixing  in  the  discharge  stability. 


6.5  Approximate  Discharge  Model 

The  jet  shown  in  Figure  61  and  Figure  62  grows  primarily  by  mixing  across  the  shear  layer. 
This  mixing  is  slow,  with  the  jet  size  remaining  (at  least  visually)  constant  for  the  3.5  cm  length  of 
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the  discharge.  We  can  exploit  this  slow  growth  rate  of  the  jet,  and  develop  a  simple  and  cost- 
effective  jet/discharge  model. 

If  we  assume  that  the  growth  of  the  jet  is  simply  due  to  diffusive  processes  (viscous  dissipa¬ 
tion,  thermal  conduction,  and  mass  diffusion),  then  the  flow  gradients  that  drive  these  processes  are 
negligible  in  the  axial,  x,  direction.  Then  the  above  conservation  equation  becomes 

dU  1  8rG 

dt  r  dr 

and  the  two-dimensional  steady-state  equation  is  converted  to  a  one-dimensional  time-dependent 
problem  that  is  much  easier  to  solve.  Within  the  limitations  of  our  assumptions,  this  approach  is 
useful  to  study  the  discharge  growth  mechanisms. 

To  implement  this  model,  we  start  at  the  cathode  plane  where  a  known  flux  of  electrons  is 
imposed  on  the  inflow  velocity  and  temperature  profile.  We  then  integrate  the  above  equation  in 
time.  Since  the  axial  velocity,  u,  varies  only  with  r,  we  can  construct  the  entire  flow  field  by  simply 
using  x=ut. 

Figure  63  plots  contours  of  temperature  and  electron  number  density  in  the  flow  field  com¬ 
puted  with  the  approximate  method  at  the  experimental  conditions.  Note  that  because  we  have  ne¬ 
glected  the  streamwise  gradients,  there  is  no  mechanism  to  promote  instability  of  the  jet.  Interest¬ 
ingly,  if  the  calculation  is  run  without  the  diffusive  transport  terms  included,  the  discharge  immedi¬ 
ately  collapses  to  an  arc.  The  jet  grows  very  slowly  because  the  thermal  diffusivity  is  extremely 
low  in  the  cold  quiescent  gas. 

Figure  64  and  Figure  65  plot  the  discharge  properties  along  the  centerline  and  at  one  axial 
location.  A  close  examination  of  the  plot  shows  that  the  electron  temperature  increases  to  very  high 
values  near  the  cathode  (at  x=0)  where  the  electrical  energy  is  added  to  a  small  number  of  electrons. 
The  high  Te  causes  the  rapid  production  of  electrons,  which  then  lowers  the  electron  temperature, 
and  a  balance  between  the  electron  density  and  temperature  is  reached.  This  figure  shows  that  the 
heavy  particle  and  vibrational  temperatures  do  not  increase  significantly,  but  the  coupling  between 
these  modes  and  the  electrons  is  responsible  for  limiting  the  electron  temperature.  The  electron 
temperature  is  approximately  10,000  K,  which  is  somewhat  lower  than  the  12,000  K  inferred  from 
the  experiments.  The  vibrational  and  translational  temperatures  closely  match  the  measured  values. 
The  predicted  electron  number  density  for  most  of  the  discharge  is  about  4xlOn  cm'3,  which  is  con¬ 
siderably  below  the  value  of  2xl012  cm"3  measured  in  the  experiments.  The  lower  value  of  ne  is 
consistent  with  the  lower  value  of  Te  predicted  by  the  simulations  because  the  electron  concentra- 
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is  very  sensitive  to  the  electron  temperature.  The  reason  for  this  difference  is  unknown  at  the  pre¬ 
sent  time,  though  the  electron  number  density  depends  on  the  details  of  the  ionization  and  energy 
transfer  rates.  In  addition,  the  turbulent  spreading  of  the  jet  would  likely  affect  the  discharge  prop¬ 
erties  by  cooling  the  edges  of  the  jet  more  efficiently,  resulting  in  more  of  the  current  being  carried 
near  the  jet  centerline.  This  would  result  in  an  increased  centerline  electron  temperature  and  con¬ 
centration. 


Density 


Figure  63.  Contours  of  temperature  and  electron  number  density  computed  using  the  approximate  method  at 

experimental  conditions. 
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Figure  64.  Centerline  temperatures  and  electron  number  density  computed  using  the  approximate  method  at  ex¬ 
perimental  conditions. 
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Figure  65.  Cross-section  of  the  computed  discharge  2  cm  from  the  cathode  at  experimental  conditions. 
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Figure  66.  Electrical  conductivity  2  cm  from  the  cathode  at  experimental  conditions. 
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Figure  67.  Variation  of  centerline  electron  number  density  with  jet  speed. 
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Figure  68.  Radial  variation  of  electron  number  density  with  jet  speed  2  cm  from  the  cathode. 
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Figure  69.  Radial  variation  of  discharge  properties.  Solid  lines:  nominal  transport  properties;  dashed  lines: 

transport  properties  increased  by  a  factor  of  10. 


The  discharge  diameter  was  measured  to  be  approximately  3.2  mm.  Figure  66  plots  the 
electrical  conductivity  at  a  location  2  cm  from  the  cathode.  Although  the  edge  of  the  discharge  is 
not  precisely  defined,  this  plot  shows  that  the  computations  produce  a  conducting  channel  of  ap¬ 
proximately  the  same  dimension. 
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Figure  67  and  Figure  68  summarize  how  the  centerline  and  radial  electron  number  density 
depend  on  the  jet  speed.  Here,  we  have  increased  and  decreased  the  jet  speed  by  a  factor  of  two 
from  the  measured  440  m/s.  The  higher  jet  speed  reduces  the  rate  of  spreading  of  the  jet,  and  there¬ 
fore  increases  the  centerline  electron  number  density.  However,  this  increase  is  minimal  and  is  a 
result  of  the  very  slow  mixing  at  the  jet  edge.  The  radial  plot  is  consistent  with  the  centerline  re¬ 
sults,  showing  an  increased  spreading  at  low  jet  speed.  Again,  the  effect  is  small. 

Based  on  the  preceding  arguments,  we  expect  turbulent  mixing  to  increase  the  spreading 
rate  of  the  jet.  We  therefore  performed  additional  simulations  with  the  transport  properties  in¬ 
creased  by  a  factor  of  10  to  mimic  increased  mixing.  Figure  69  compares  the  radial  variation  of  the 
temperatures  and  electron  number  density.  The  jet  radius  increases  only  minimally.  Clearly,  this 
simple  approach  of  increasing  the  molecular  transport  properties  does  not  produce  substantial  in¬ 
creased  jet  spreading.  This  is  because  turbulent  mixing  is  a  large-scale  convective  phenomenon.  A 
more  complete  characterization  of  the  jet  is  required  to  incorporate  turbulent  mixing  into  the  model. 


6.6.  Conclusions 

We  have  developed  a  computational  methodology  to  analyze  Stanford  University  DC  and 
pulsed  discharge  experiments.  A  complete  axisymmetric  simulation  including  an  11 -species  finite- 
rate  air  kinetics  model  and  a  three-temperature  energy  formulation  was  used.  A  three-temperature 
model  is  used  to  allow  separate  heavy-particle,  vibrational  and  electron  translation  temperatures  of 
the  plasma.  The  self-consistent  effective  binary  diffusion  model  is  used  for  the  molecular  transport 
properties.  The  discharge  is  modeled  as  a  conducting  channel,  allowing  a  direct  enforcement  of 
total  current  conservation. 

The  governing  equations  are  solved  with  a  semi-implicit  finite-volume  method.  The  “fast” 
source  terms  are  linearized  in  the  time  integration,  resulting  in  a  method  that  at  least  50  times  faster 
than  an  explicit  method.  A  fully  implicit  method  is  about  100  times  more  expensive  per  time  step, 
and  does  not  significantly  increase  the  stable  time  step  size.  Thus,  for  this  problem  that  is  domi¬ 
nated  by  very  fast  (order  of  nanosecond)  time-scales,  our  semi-implicit  approach  is  optimal  and  en¬ 
ables  simulations  of  this  type.  The  convective  fluxes  are  represented  by  standard  second-order  ac¬ 
curate  upwind  fluxes,  and  the  viscous  fluxes  use  a  central  difference  formulation. 

The  plasma  torch  simulations  using  this  CFD  model  produced  an  unstable  jet  with  large- 
scale  vortical  structures.  Contrary  to  the  experiments,  the  non-uniformities  pinch  the  discharge  and 
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transition  to  an  arc  occurs.  We  postulate  that  turbulent  mixing  in  the  experiments  stabilizes  the  jet 
and  prevents  arcing. 

To  study  the  discharge  spreading,  we  reformulated  the  model,  freezing  the  gradients  in  the 
axial  direction.  This  prevents  formation  of  large  vortical  structures,  and  we  obtain  a  stable  diffuse 
discharge.  The  size  of  the  computed  electrically  conducting  channel  matches  the  experiments,  and 
the  electron  number  densities  are  fairly  well  predicted  considering  that  they  are  highly  sensitive  to 
the  electric  field  or  equivalently  to  the  electron  temperature.  The  vibrational  and  translational  tem¬ 
peratures  also  agree  well. 

Using  the  approximate  formulation,  we  varied  the  jet  speed  and  transport  properties  to  study 
how  the  jet  spreads.  We  find  that  the  very  low  thermal  diffusivity  in  the  cold  external  gas  makes 
the  jet  spreading  essentially  insensitive  to  these  changes.  The  lack  of  turbulent  mixing  in  our  model 
is  a  major  limitation  and  future  work  should  address  this  issue. 

The  effect  of  turbulent  mixing  appears  to  be  to  stabilize  the  jet  under  the  high  flow  rate 
cases  simulated.  However,  we  expect  that  at  lower  flow  rates  where  the  jet  does  not  roll  up  into 
vortical  structures,  the  discharge  will  be  stable.  Under  these  conditions,  the  flow  is  very  closely 
represented  by  the  approximate  discharge  model,  which  produces  a  stable  diffuse  discharge.  Thus 
although  our  simulations  result  in  arc  formation  at  the  experimental  conditions,  this  does  not  imply 
that  our  full  axisymmetric  simulations  are  unstable  at  all  flow  conditions. 

The  modeling  of  turbulent  mixing  in  these  flows  is  very  challenging.  We  are  investigating 
the  approach  of  Huang  as  one  means  of  increasing  the  heat  transfer  across  the  jet  interface.  How¬ 
ever,  this  approach  is  largely  empirical  and  will  not  result  in  a  truly  predictive  capability  to  simulate 
this  type  of  flow.  Other  turbulence  modeling  approaches  for  free  jets  can  be  used  that  may  result  in 
a  more  accurate  and  general  model.  This  should  be  supported  by  a  systematic  study  of  varying  the 
torch  flow  properties  in  experiments  to  better  characterize  the  role  of  turbulence  in  these  flows. 
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7.  Measurements  of  Ion  Concentration  by  Cavity  Ring-Down  Spectroscopy 

7.1.  Introduction 

Measurements  of  ion  and/or  electron  number  density  are  needed  to  characterize  experiments 
and  validate  models  for  the  atmospheric  pressure  air  and  nitrogen  plasmas  presented  in  this  report. 
The  ionization  fraction  in  these  atmospheric  pressure  plasmas  is  of  the  order  of  10'5,  so  that  sensi¬ 
tive  measurement  techniques  are  required.  Physical  probes  tend  to  disturb  the  system,  and  tech¬ 
niques  such  as  Thomson  scattering  and  interferometry  do  not  readily  provide  results  with  high  spa¬ 
tial  resolution.  Optical  techniques  that  measure  ion  concentrations  are  widely  used.  Of  these, 
emission  provides  information  only  on  excited  species,  fluorescence  suffers  from  quenching  effects 
and  optical  interference  that  complicate  interpretation,  and  absorption  often  lacks  sensitivity.  Cav¬ 
ity  ring-down  spectroscopy  (CRDS),  on  the  other  hand,  is  a  sensitive  line-of-sight  laser  absorption 
technique  that  has  been  used  to  measure  species  concentrations  in  low-pressure  plasmas  [78-82].  In 
particular,  the  N2+  ion  has  been  studied  in  low-pressure  hollow  cathode  sources  [83,  84].  We  have 
selected  CRDS  as  our  primary  diagnostic  tool  to  characterize  ion  concentrations.  In  this  section,  we 
describe  the  use  of  CRDS  to  measure  ion  concentrations  in  atmospheric  pressure  discharges.  By 
implementing  CRDS  in  its  “standard”  form,  we  perform  spatially  resolved  ion  concentration  meas¬ 
urements.  We  also  develop  a  temporally  resolved  variant  of  CRDS,  which  we  used  to  study  ion 
recombination  in  pulsed  plasmas.  Measurements  have  been  performed  in  both  air  and  nitrogen 
plasmas.  In  nitrogen  plasmas,  CRDS  measurements  of  the  dominant  N2+  ion  enable  characteriza¬ 
tion  of  plasma  ionization.  In  the  air  plasmas,  CRDS  measurements  of  the  NO+  ion  were  performed. 

An  overview  of  the  CRDS  technique  is  provided  in  Section  7.2,  including  a  discussion  of 
temporally  resolved  CRDS.  Section  7.3  presents  the  experimental  schemes  used  for  spatially  re¬ 
solved  ion  concentration  measurements  of  the  N2+  ion  in  DC  discharges,  as  well  as  temporally  re¬ 
solved  N2+  ion  concentration  measurements  in  pulsed  discharges.  Measurement  results,  and  discus¬ 
sion,  are  provided.  To  aid  in  interpreting  results,  a  collisional  radiative  (CR)  model  is  used  to  com¬ 
pute  population  fractions  and  to  relate  the  measured  ion  concentrations  to  electron  number  densi¬ 
ties.  The  inferred  electron  number  density  profiles  are  compared  with  electrical  measurements,  and 
the  nonequilibrium  nature  of  the  plasma  is  discussed.  Section  7.4  discusses  CRDS  measurements 
of  the  NO+  ion  in  air  plasmas.  The  experimental  scheme  and  a  discussion  of  results  are  presented. 
Conclusions  are  provided  in  Section  7.5. 
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7.2.  Cavity  ring-down  spectroscopy 

Cavity  ring-down  spectroscopy  has  become  a  widely  used  method  in  absorption  spectros¬ 
copy  owing  primarily  to  its  high  sensitivity.  Detailed  reviews  of  the  technique  may  be  found  in 
References  [85,  86].  Essentially,  a  laser  beam  is  coupled  into  a  high-finesse  optical  cavity  contain¬ 
ing  a  sample,  where  it  passes  many  times  between  the  mirrors.  As  the  light  bounces  back  and  forth 
inside  the  cavity,  its  intensity  decays  (rings  down)  owing  to  sample  absorption,  particle  scattering 
loss  (generally  negligible),  and  mirror  transmission  loss.  A  photodetector  is  used  to  measure  the 
ring-down  signal,  which  is  fitted  to  yield  the  sample  loss.  The  technique  affords  high  sensitivity 
owing  to  a  combination  of  long  effective  path  length  and  insensitivity  to  laser  energy  fluctuations, 
and  therefore  is  well  suited  to  the  detection  of  trace  species  in  plasmas.  Under  appropriate  condi¬ 
tions,  the  laser  lineshape  may  be  neglected,  and  the  ring-down  signal  S(t )  decays  exponentially  [87, 
88]  as: 


S(t )  =  S0  exp[-f/r] 
l/r=  y  [/^(vJ  +  Cl-*)] 


(33) 


where  r  is  the  1/e  time  of  the  decay  (termed  the  ring-down  time),  c  is  the  speed  of  light,  /  is  the  cav¬ 
ity  length,  labs  is  the  absorber  column  length,  k(y)  is  the  absorption  coefficient,  vL  is  the  laser  fre¬ 
quency,  and  1-R  is  the  effective  mirror  loss  (including  scattering  and  all  cavity  losses).  Generally, 
the  measured  ring-down  signal  is  fit  with  an  exponential,  and  the  ring-down  time  t  is  extracted. 
Combining  x  with  the  ring-down  time  to  measured  with  the  laser  detuned  from  the  absorption  fea¬ 
ture  allows  a  determination  of  the  sample  absorbance,  and  hence  absorption  coefficient: 


Abs  =  labsk{vL)  ■ 


(34) 


We  minimize  any  potential  laser  lineshape  dependence  by  tuning  the  laser  frequency  across 
an  absorption  line,  and  measuring  the  frequency-integrated  absorption  coefficient  [88]. 


The  prior  discussion  of  CRDS  has  implicitly  assumed  that  the  sample  concentration  (and  as¬ 
sociated  absorption  loss)  is  independent  of  time,  as  would  be  the  case  in  the  DC  discharge.  How¬ 
ever,  a  more  complex  approach  is  required  to  measure  time  dependent  concentrations  in  pulsed  dis¬ 
charges,  since  the  ion  concentration  varies  over  the  duration  of  the  optical  ring-down  (decay  of  light 
in  the  cavity).  It  might  be  tempting  to  consider  using  lower  reflectivity  mirrors  with  shorter  ring- 
down  times  so  that  the  losses  may  be  treated  as  constant  over  the  ring-down,  but  the  sensitivity  of 
such  an  approach  is  inferior  [87].  Although  a  number  of  kinetics  studies  have  been  performed  with 
CRDS,  nearly  all  of  these  experiments  study  processes  that  are  slow  compared  to  experimental 
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ring-down  times.  An  exception  is  the  work  of  Brown  et  al.  [89],  who  perform  gas-phase  measure¬ 
ments  in  cases  where  the  populations  do  change  over  the  duration  of  the  ring-down.  We  follow  a 
related  approach  to  measure  ion  recombination  in  a  plasma  over  time-scales  comparable  to  the  ring- 
down  time  (microseconds).  For  the  case  of  a  time-dependent  absorption,  the  ring-down  signal  S(t) 
may  be  written  as  [89]: 


S(t)  =  S0exp 


c 

T 


\Kv,t)iabfi  +0  ~R)t 
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(35) 


where  the  absorption  coefficient  now  has  a  time  dependence.  Rearranging  equation  (35)  leads  to  an 
expression  for  the  absorbance  as  a  function  of  time: 


Abs{t)  =  k(v,t  )labs 


l_d_ 
c  dt 


In 


rS{t)' 


V  °0  J 


(1-R) 


(36) 


The  derivative  (local  slope)  of  the  logarithm  of  the  ring-down  signal  is  proportional  to  the 
loss  (sample  plus  empty  cavity)  at  that  time.  To  directly  obtain  time-dependent  concentrations, 
Brown  et  al.  analyzed  their  data  with  this  method.  We  choose  to  follow  an  alternative  approach  in 
which  the  ring-down  signal  is  divided  into  a  series  of  time-windows,  each  of  which  is  fit  to  an  ex¬ 
ponential  decay  (ring-down  time).  We  believe  that  this  approach  is  less  noisy  because  it  avoids  dif¬ 
ferentiation. 


7.3.  N?+  Measurements 


7.3.1.  Experimental  Setup 

7.3.1  .a.  Atmospheric  Pressure  Discharge 

We  have  developed  a  compact  atmospheric  pressure  plasma  source  for  diagnostic  develop¬ 
ment.  The  discharge  may  be  operated  with  both  nitrogen  and  air.  A  photograph  of  the  nitrogen  dis¬ 
charge  with  a  schematic  representation  of  the  ring-down  cavity  is  shown  in  Figure  70.  Nitrogen  is 
injected  through  a  flow  straightener  and  passes  through  the  discharge  region  with  a  velocity  of 
about  20  cm/s.  The  discharge  is  formed  between  a  pair  of  platinum  pins  (separation  0.85  cm)  that 
are  vertically  mounted  on  water-cooled  stainless-steel  tubes.  The  discharge  is  maintained  by  a  DC 
current  supply  (imax=250  mA)  in  a  ballasted  circuit  (Rb=9.35  kQ).  The  pins  are  brought  together  to 
ignite  the  discharge,  and  are  then  separated  using  a  translation  stage.  The  position  of  the  discharge 
is  observed  to  be  stable  and  reproducible.  The  discharge  is  contained  within  a  Plexiglas  cylinder 
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through  the  top  to  avoid  accumulation  of  undesirable  by-products  of  the  discharge  (such  as  ozone  or 
oxides  of  nitrogen),  and  enable  passage  of  the  laser  beam  through  the  discharge.  A  second  transla¬ 
tion  stage  is  used  to  displace  the  entire  discharge  cylinder  relative  to  the  optical  axis  in  order  to  ob¬ 
tain  spatial  profiles. 


to  data 
acquisition 


Figure  70  Photograph  of  the  atmospheric  pressure  nitrogen  discharge  and  schematic  diagram  of  the  ring-down 
cavity.  Electrode  separation:  0.85  cm.  Discharge  current:  187  mA. 


To  explore  the  repetitively  pulsed  approaches,  we  connect  a  high-voltage  pulser  in  parallel  to 
the  DC  discharge  circuit.  The  pulser  is  capacitively  coupled  to  the  discharge  so  that  it  is  isolated 
from  the  DC  supply.  The  DC  field  serves  to  give  a  baseline  of  ionization,  and  to  heat  the  gas.  We 
operate  the  high-voltage  pulser  (pulse  width-1 0ns,  pulse  voltage-8  kV)  at  10  Hz  so  that  it  may  be 
synchronized  relative  to  the  laser.  At  this  repetition  rate  the  plasma  equilibrates  between  high- 
voltage  pulses  so  that  the  behavior  during  and  following  each  pulse  is  not  affected  by  the  presence 
of  other  pulses. 

7 . 3 . 1  ,b .  CRDS  Measurements. 

We  study  the  N2+  ion  by  probing  the  (0,0)  band  of  its  first  negative  system  (B  2Zu+  -  X  2Ig+) 
in  the  vicinity  of  391  nm.  We  select  this  spectral  feature  because  it  is  comparatively  strong  and  op¬ 
tically  accessible.  The  optical  layout  is  shown  in  Figure  71 .  An  OPO  system  (doubled  idler)  is  used 
as  the  light  source  (repetition  rate  =  10  Hz,  pulse  width  -  7  ns,  pulse  energy  -  3  mJ, 
linewidth  -  0.14  cm'1).  The  output  from  the  OPO  passes  through  a  Glan-Taylor  polarizer  to  attenu¬ 
ate  the  energy,  and  several  beam  shaping  optics  to  approximately  mode-match  the  beam  into  the 
cavity.  We  use  a  linear  cavity  of  75-cm  length  with  50-cm  radius-of-curvature  (ROC)  mirrors.  The 
selection  of  cavity  geometry  is  discussed  below.  Typically,  about  100  jiJ  per  pulse  is  incident  on  the 
back  face  of  the  entrance  ring-down  mirror.  The  ring-down  signal  is  collected  behind  the  output 
mirror  with  a  fast  photomultiplier  tube  (Hammamatsu-R1104),  which  we  filter  against  the  pump 
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mirror  with  a  fast  photomultiplier  tube  (Hammamatsu-R1104),  which  we  filter  against  the  pump 
laser  and  other  luminosity  with  two  narrow-band  interference  filters  (CVI-F10-390-4-1).  The  PMT 
signals  are  passed  to  a  digitizing  oscilloscope  (HP  545 10A,  250  MHz  analog  bandwidth,  8-bit  ver¬ 
tical  resolution)  and  are  read  to  computer  with  custom  data  acquisition  software.  In  a  typical  ring- 
down  spectrum,  16  or  32  decay  curves  are  averaged  at  each  wavelength,  and  the  resulting  wave¬ 
form  is  fitted  with  an  exponential  to  yield  the  ring-down  time  x.  For  the  DC  measurements,  the  por¬ 
tion  of  the  ring-down  signal  used  in  the  fit  is  that  in  between  90%  and  10%  of  the  peak  (initial)  sig¬ 
nal  amplitude.  The  detuned  ring-down  time  to  is  determined  with  the  laser  tuned  off  the  absorption 
features.  Spectral  scans  use  a  step-size  of  0.001  nm.  When  performing  spatial  scans,  we  use  step- 
sizes  of  0.2  mm. 

The  CRDS  set-up  used  for  the  pulsed  measurements  differs  only  in  terms  of  data  fitting  and 
timing.  Because  sample  absorption  loss  is  no  longer  constant  during  the  measurement,  the  time  de¬ 
pendent  equation  (36)  is  used  to  fit  the  data.  Rather  than  computing  derivatives  we  fit  a  series  of 
line  segments  to  the  logarithm  of  the  ring-down  signal.  The  fitting  windows  are  1-p.s  in  length. 
This  time  interval  represents  a  good  compromise  in  making  the  window  short  compared  to  the 
timescale  of  the  process  studied  yet  affording  an  acceptable  signal-to-noise  level.  Because  our 
time  windows  are  long  compared  to  the  pulse  length  (~  10  ns),  we  do  not  resolve  the  build-up  of 
ionization  that  occurs  during  the  time  the  pulse  is  on;  yet  we  are  able  to  resolve  the  subsequent  re¬ 
combination.  We  synchronize  the  laser  relative  to  the  firing  of  the  discharge  with  an  external  tim¬ 
ing  circuit.  In  order  to  obtain  concentration  information  at  different  times  relative  to  the  firing  of 
the  high  voltage  pulse,  we  vary  the  delay  time  between  the  firing  of  the  laser  and  the  firing  of  the 
high  voltage  pulser.  Delay  times  (Tpuise-TLaser)  of -10, -8, -6, -4, -2, -1,0, 1,2  ps  are  used. 
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Figure  71  Schematic  diagram  of  CRDS  set-up.  The  ring-down  cavity  has  a  length  of  0.75  m,  and  uses  0.5-m  radius 
of  curvature  mirrors.  An  OPO  is  used  as  the  light  source,  and  a  photomultiplier  tube  (PMT)  detects  the  light  exit¬ 
ing  the  cavity. 


Implementing  CRDS  in  the  atmospheric  plasma  requires  special  care  in  the  choice  of  cavity 
geometry.  For  a  linear  cavity  formed  with  mirrors  of  equal  radius-of-curvature,  the  cavity  geometry 
is  determined  by  the  dimensionless  g-parameter,  defined  as  unity  minus  the  cavity  length  divided  by 
the  mirror  radius  of  curvature  [90],  Initial  attempts  to  form  a  ring-down  cavity  with  g=  0.875 
(length  75  cm,  6-m  ROC  mirrors),  resulted  in  distorted  and  irreproducible  profiles,  owing  to  beam 
steering  from  index-of-refraction  gradients  (similar  to  a  mirage).  Recent  work  by  Spuler  et  al.  [91] 
simulates  the  effect  of  cavity  geometry  on  beam  propagation  in  CRDS  experiments.  Their  results 
indicate  that  a  g-parameter  of  about  -0.5  represents  a  good  compromise  between  beam  waist  and 
beam  walk  in  environments  where  beam  steering  may  be  present.  Accordingly,  we  form  a  cavity  of 
length  75  cm,  with  50-cm  ROC  mirrors  (Research  Electro  Optics).  No  beam  steering  is  detected 
with  this  geometry.  Qualitatively,  this  geometry  tends  to  recenter  deviated  beams,  whereas  the  more 
planar  geometry  is  not  as  effective. 


7.3.1. c.  Electrical  Measurements. 

We  also  determine  the  electron  number  density  in  the  discharge  by  an  electrical  conductivity 
approach.  For  the  DC  discharge,  we  measure  the  time  independent  discharge  current  and  electric 
field,  and  use  Ohm’s  law  to  compute  the  product  of  the  average  electron  number  density  and  col¬ 
umn  area.  The  electric  field  is  found  as  the  slope  of  the  discharge  voltage  versus  electrode  separa¬ 
tion.  We  write  Ohm’s  law  as  j=i/area=(nee2/meEveh)E  where  veh  is  the  average  collision  frequency 
between  electrons  and  heavy  particles.  Because  of  the  low  ionization  fraction  (<~10'5),  veb  is  domi- 
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dominated  by  collisions  with  neutrals,  so  we  can  write  veh=nngeQen,  where  nn=P/kTg  is  the  number 
density  of  neutrals,  ge=(8kTe/7ime)0  5  is  the  thermal  electron  velocity,  and  Qen  is  the  average  momen¬ 
tum  transfer  cross-section  for  electron-nitrogen  collisions.  We  determine  Te  by  means  of  a  colli- 
sional  radiative  model  [5],  with  input  parameters  (vibronic  ground  state  population  and  rotational 
temperature)  obtained  from  the  CRDS  measurements.  Qen  is  determined  as  a  function  of  Te  from 
tabulated  values  [15]. 

We  follow  an  analogous  approach  to  determine  the  time  varying  electron  number  density  in 
the  pulsed  discharge.  In  this  case  we  measure  the  time-dependent  current  and  electric  field,  and  use 
these  to  determine  temporally  resolved  electron  number  densities. 


7.3.2.  N2+  Measurements  -  Results  and  Discussion 

7.3. 2. a.  N?+  Ring-Down  Spectra. 

N2+  ring-down  spectra  are  recorded  as  a  function  of  discharge  current  and  position.  We  dis¬ 
cuss  the  spectra  in  terms  of  measurement  accuracy  and  detection  sensitivity. 

Figure  72  shows  measured  and  simulated  absorption  spectra  in  the  vicinity  of  the  (0,0) 
bandhead  of  the  first  negative  system.  Rotationally  resolved  lines  from  the  P  and  R  branches  are 
visible.  The  lines  are  identified  using  tabulated  line  locations  [92,  93],  and  are  labeled  with  the  an¬ 
gular  momentum  quantum  number  N"  of  the  lower  state.  The  displayed  spectrum  is  recorded  along 
the  discharge  centerline,  at  a  current  of  187  mA,  and  averages  16  shots  at  each  wavelength.  The 
experimental  spectrum  is  plotted  in  terms  of  single-pass  cavity  loss  and  illustrates  the  high  sensitiv¬ 
ity  attainable  with  the  CRDS  technique.  The  cavity  loss  is  the  sum  of  mirror  reflective  loss  and 
sample  absorptive  loss.  (The  Rayleigh  scattering  losses  are  computed  to  be  negligible.)  Near  the 
bandhead,  the  signal  is  ~  280  ppm/pass,  the  baseline  reflective  loss  is  ~  200  ppm/pass,  and  the  base¬ 
line  noise  is  ~  5  ppm/pass,  so  that  the  signal-to-noise  (S/N)  ratio  is  ~  56.  The  S/N  ratio  suggests  an 
absorbance-per-pass  sensitivity  of  about  1  ppm,  which  corresponds  to  a  detection  limit  of  about 
7xl010  cm'3  for  N2+  ions  at  our  experimental  conditions.  The  baseline  reflective  loss  of  ~  200 
ppm/pass  corresponds  to  a  mirror  reflectivity  of  ~  0.9998  which  is  in  accord  with  the  manufac¬ 
turer’s  specifications. 
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Figure  72  Measured  and  simulated  N2+  absorption  spectra  near  the  (0,0)  bandhead  of  the  first  negative  system. 

Lines  from  the  P  and  R  branches  are  identified. 

Figure  73  shows  an  expanded  view  of  the  P(28)  and  R(l)  lines  after  baseline  subtraction. 
Fitted  Voigt  peaks  (constrained  to  have  the  same  shape  and  width)  are  shown  with  solid  lines,  and 
their  sum  is  shown  with  a  dotted  line.  For  the  P(28)  lines,  the  doublet  structure  due  to  the  unpaired 
electron  is  apparent.  The  fit  yields  a  doublet  spacing  of  0.005  nm,  in  good  agreement  with  the  lit¬ 
erature  [92,  93].  The  R(l)  lines  are  close  to  the  detection  limit  and  correspond  to  a  N2+  X  state 
population  of  about  1010  cm'3.  Their  splitting  is  much  less  than  the  linewidths  and  is  not  resolved. 
The  fitted  FWHM  of  each  peak  is  0.0042  nm,  or  7.9  GHz,  which  is  consistent  with  an  expected 
thermally  broadened  linewidth  of  ~  7  GHz  and  a  measured  laser  linewidth  of  ~  4  GHz.  Our  calcu¬ 
lations  indicate  that  for  these  linewidths  and  absorption  parameters,  a  measurement  of  the  area  of  an 
absorbing  feature  will  be  negligibly  affected  by  finite  laser  bandwidth  [88].  Note  that  unlike  in 
conventional  absorption,  such  an  assumption  is  not  always  correct  in  CRDS.  The  assumption  was 
verified  in  the  present  work  by  fitting  progressively  earlier  parts  of  the  ring-down  signal  and  con¬ 
firming  that  the  apparent  absorbance  remains  constant. 
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Figure  73  Expanded  view  of  P(28),  and  R(l)  lines  from  Fig.  3.  Background  absorption  has  been  subtracted. 
Voigt  peaks  fitted  to  the  doublet  are  shown  with  solid  lines,  while  their  sum  is  shown  with  a  dotted  line. 


7.3.2.b.  Spatial  Profiles  of  Ion  Concentration  and  Electron  Number  Density. 

We  obtain  spatial  profiles  of  the  N2+  concentration  by  displacing  the  discharge  perpendicu¬ 
larly  to  the  optical  axis.  CRDS  is  a  path-integrated  technique  and  the  discharge  has  axial  symme¬ 
try.  We  verify  the  symmetry  of  the  discharge  by  performing  measurements  with  the  plasma  rotated 
by  90  degrees,  and  find  that  the  cases  have  <2%  deviation.  We  use  an  Abel  inversion  to  recover  the 
radial  N2+  concentration  profile.  The  concentration  measurements  are  based  on  the  (frequency  inte¬ 
grated)  area  of  the  lines  P(9)-P(17)  in  the  (0,0)  bandhead  vicinity.  We  use  tabulated  line  strengths 
from  References  [92,  93],  Figure  74  shows  concentration  profiles  determined  for  different  values  of 
current  (i=52,  97,  142,  and  187mA).  We  find  peak  (centerline)  N2+  concentrations  of  7.8xlOn, 
1.5xl012,  2.4xl012,  and  3.6xl012  cm'3  for  i=52,  97,  142,  and  187  mA  respectively.  The  shape  of 
the  concentration  profile  remains  approximately  uniform  at  the  different  conditions,  though  we  ob¬ 
serve  that  the  radial  half-maximum  values  increase  slightly  with  current.  We  find  radial  half- 
maximums  of  0.80,  0.82,  0.93,  and  1.05  mm  for  i=52,  97,  142,  and  187  mA  respectively. 
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Figure  74  Radial  concentration  profiles  of  N2+  measured  by  CRDS  in  an  atmospheric  pressure  glow  discharge. 
Experimental  data  points  are  joined  with  line  segments  for  visual  clarity. 


The  error  bars  on  N2+  concentration  represent  one  standard  deviation  (lcr).  They  primarily 
arise  from  the  uncertainties  in  relating  the  measured  population  of  several  rotational  levels  in  the 
ground  vibronic  state,  to  the  overall  population  of  N2+.  Because  the  discharge  is  out  of  equilibrium, 
this  relationship  depends  on  the  how  the  rotational,  vibrational,  and  electronic  energy  levels  are 
populated.  The  rotational  levels  are  equilibrated  at  the  gas  temperature  owing  to  fast  collisional 
relaxation.  The  rotational  temperatures  used  in  the  analysis  are  obtained  from  Boltzmann  plots,  and 
are  Tr=3100,  3600,  4150,  and  4700  K  for  currents  of  i=52,  97,  142,  and  187  mA  respectively.  The 
vibrational  and  electronic  energy  levels  are  out  of  equilibrium,  and  a  collisional-radiative  (C-R) 
model  [5]  is  used  to  determine  the  fraction  of  the  population  in  the  ground  vibronic  state,  and  pre¬ 
dicts  0.37±0.02,  0.3510.02,  0.3310.02,  and  0.3110.02  for  i=52,  97,  142,  and  187  mA  respectively. 
Combining  the  rotational  temperature  uncertainties  with  those  from  the  C-R  model  and  those  from 
the  Abel  inversion  (~4%),  results  in  an  overall  experimental  uncertainty  in  concentration  of  -10%. 

The  spatial  resolution  of  our  measurements  is  determined  by  the  spatial  step-size  (0.2  mm). 
To  justify  this  claim,  we  need  to  verify  that  the  dimension  of  our  laser  beam  waist  does  not  influ¬ 
ence  the  measured  spatial  profiles.  The  simulations  by  Spuler  et  al.  [91]  indicate  that  our  expected 
beam  waist  is  approximately  160-320  microns,  depending  on  the  level  of  mode  matching  achieved. 
Deconvoluting  the  broader  case  has  an  effect  of  only  about  1%  (0.02  mm)  on  the  measured  profiles, 
which  is  negligible  compared  to  the  spatial  step-size.  Therefore  the  resulting  spatial  resolution  is 
about  0.2  mm. 
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We  incorporate  the  electrical  measurements  by  comparing  the  electron  number  density  in¬ 
ferred  from  the  CRDS  ion  measurements,  to  the  electron  number  density  from  the  electrical  con¬ 
ductivity  approach.  To  infer  electron  number  densities  from  the  CRDS,  we  need  to  know  the  frac¬ 
tion  of  positive  ions  that  are  N2+.  At  our  conditions,  the  C-R  model  predicts  that  96,  93,  89  and 
85%  of  ions  are  N2+,  and  the  remainder  is  N+,  for  i=52,  97,  142,  and  187  mA  respectively.  By 
charge  balance,  the  sum  of  the  N2+  and  N+  concentrations  equals  the  electron  number  density.  We 
convert  the  N2+  concentration  profiles  (found  by  CRDS)  to  electron  number  density  profiles  using 
these  percentages.  In  order  to  determine  electron  number  densities  from  the  conductivity  measure¬ 
ments  (which  yield  the  product  of  average  electron  number  density  with  area),  we  assume  that  the 
shape  of  the  electron  number  density  profile  is  the  same  as  that  for  the  ions.  The  electron  number 
densities  (at  the  radial  half-maximum)  found  in  this  way  from  electrical  measurements  are  com¬ 
pared  with  those  inferred  from  the  CRDS  ion  measurements  in  Table  4.  The  values  are  plotted  in 
Figure  75.  The  uncertainty  in  the  electrical  measurement  (10%)  is  primarily  from  uncertainties  in 
the  momentum  transfer  cross-section  (5%),  the  discharge  area  (4%),  and  the  average  gas  tempera¬ 
ture  (8%).  Column  4  of  Table  4  shows  that  the  electron  number  densities  found  from  optical  and 
electrical  measurements  overlap  within  their  error  bars.  This  excellent  agreement  gives  us  confi¬ 
dence  in  our  results  for  the  electron  number  density. 

Table  4  Comparison  of  electron  number  densities  (at  the  radial  half-maximum)  inferred  by  CRDS  to  those 
found  by  electrical  measurement,  for  the  DC  discharges.  The  last  column  is  the  ratio  of  the  electron  number 
density  inferred  by  CRDS  to  that  found  from  electrical  measurement. 


i(mA) 

iVcRDsfcm'3) 

ne.Eiec(cm'3) 

CRDS/Elec 

52 

4. 1+0.4  x  10" 

3. 8+0.4  x  10“ 

1.08±0.16 

97 

8. 2+0.8  x  1011 

7.8+0.8  x  1011 

1.05±0.16 

142 

1. 4+0.1  x  1012 

1 .4+0.1  x  1012 

0.96+0.14 

187 

2. 1+0.2  x  1012 

2. 0+0.2  x  1012 

1.06±0.16 

7.3.2.C.  Temporal  Profiles  of  Ion  Concentration  and  Electron  Number  Density. 

Figure  76  shows  ring-down  traces  obtained  with  and  without  firing  the  high  voltage  pulse, 
and  with  the  laser  tuned  to  the  N2+  bandhead.  In  the  absence  of  the  high-voltage  pulse  (dashed  line) 
the  absorption  losses  are  constant  in  time,  and  the  signal  decays  as  a  single-exponential.  In  the 
trace  with  the  pulse  (solid  line),  the  light  decays  more  steeply  after  the  pulse,  reflecting  an  increased 
concentration  of  N2+.  The  spike  in  the  latter  trace  coincides  with  the  firing  of  the  pulse,  and  is 
caused  by  RF  interference  generated  by  the  pulser.  To  verify  that  we  are  observing  changes  in  the 
N2+  concentration,  we  examine  the  analogous  traces  but  with  the  laser  detuned  from  the  absorption 
band  -  Figure  77.  These  traces  confirm  that  the  only  effect  of  the  high  voltage  pulse  on  the  ring- 
down  system  is  to  generate  the  interference  spike.  We  analyze  these  traces  to  determine  over  what 
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region  the  interference  spike  affects  the  data.  We  vary  the  delay  of  the  high-voltage  pulse  relative  to 
the  laser  shot  so  that  we  can  obtain  ion  concentrations  at  different  times. 
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Figure  75  Electron  number  densities  (at  the  radial  half-maximum)  as  a  function  of  discharge  current.  Number 
densities  are  derived  from  CRDS  ion  measurements  (squares),  and  from  electrical  measurement  (circles). 


We  quantify  the  time-varying  N2+  concentration  using  equation  (36)  with  a  1-ps  window. 
This  time  interval  represents  a  good  compromise  in  making  the  window  short  compared  to  the  time- 
scale  of  the  process  studied  yet  affording  an  acceptable  signal-to-noise  level.  The  empty-cavity 
losses  (mirror  reflectivity)  are  found  from  the  ring-down  signals  with  the  laser  detuned,  and  are 
subtracted  in  the  analysis.  Using  tabulated  linestrengths  and  the  discharge  dimensions,  we  find  the 
absolute  N2+  centerline  concentrations  as  a  function  of  time.  Figure  78  presents  the  time  varying 
concentrations  (symbols).  The  error  bars  reflect  uncertainties  in  the  population  fractions,  as  well  as 
uncertainty  associated  with  a  possible  change  in  shape  of  the  concentration  profile.  The  latter  un¬ 
certainty  is  estimated  by  chemical  kinetic  considerations.  One  microsecond  after  the  pulse,  the  N2+ 
concentration  is  ~1 .5  x  1013  cm"3,  and  then  N2+  recombines  to  the  DC  level  in  about  10  ps.  The  DC 
level  is  found  by  analyzing  the  pulsed  data  at  sufficiently  long  time  delays  after  the  pulse,  and  is 
consistent  with  that  found  in  the  DC  plasma  without  the  pulser. 
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Figure  76  Experimental  ring-down  traces  with  the  laser  tuned  to  the  N2+  absorption  bandhead  (inset)  with  the 
high-voltage  pulse  (solid  line)  and  without  the  high-voltage  pulse  (dashed  line). 


Figure  77  Experimental  ring-down  traces  with  the  laser  tuned  away  from  the  N2+  absorption  (inset).  We  slightly 
scale  (<5%)  the  amplitude  of  the  traces  for  visual  clarity.  The  detuned  trace  (dashed  line)  is  offset  by  0.2  mV  to 
make  it  more  visible. 


For  the  pulsed  discharge,  we  also  determine  the  electron  concentration  by  measuring  the 
electrical  conductivity.  The  temporally  resolved  electron  concentrations  are  shown  with  a  swath  in 
Figure  78.  The  uncertainty  in  the  DC  electron  concentration  reflects  uncertainties  in  the  profile 
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shape,  the  momentum  transfer  cross-section,  and  the  gas  temperature.  The  collisional-radiative 
model  predicts  that  N2+  is  the  dominant  ion  produced  by  the  pulse.  Thus,  the  agreement  between  the 
time-dependent  electron  and  N2+  concentrations  during  plasma  recombination  verifies  the  tempo¬ 
rally  resolved  CRDS  measurement.  The  measured  recombination  time  is  consistent  with  reported 
[10]  dissociative  recombination  rate  coefficients  for  N2+  (approximately  5  x  10'8  cm3/s). 


Figure  78  CRDS  measurements  of  N2+  concentrations  (circles)  and  conductivity  measurements  of  electron 
densities  (swath)  versus  time  following  the  firing  of  a  high-voltage  pulse  in  an  atmospheric  pressure  nitrogen  DC 
plasma.  The  DC  level  of  N2+  concentration  found  by  CRDS  is  shown  with  a  hatched  bar.  The  inset  shows  the 
ring-down  signals  (plotted  on  a  semi-log  scale)  with  the  HV  pulse  (solid),  and  without  the  HV  pulse  (dotted). 


7.3.2.d.  Discharge  Nonequilibrium 

To  have  a  measure  of  the  degree  of  non-equilibrium  in  the  DC  discharges,  we  examine  the 
ratio  of  the  measured  electron  number  density  (at  the  radial  half-maximum)  to  the  LTE  electron 
number  density  at  the  corresponding  gas  temperature.  These  ratios  are  given  in  column  3  of  Table 
5  for  the  four  conditions  studied  in  the  DC  discharge.  The  measured  ion  and  electron  concentra¬ 
tions  in  the  discharge  are  significantly  higher  than  those  corresponding  to  LTE  conditions  at  the 
same  gas  temperature.  The  results  quantify  the  degree  of  ionizational  nonequilibrium  in  the  dis¬ 
charges.  At  higher  values  of  discharge  current  the  LTE  concentration  of  charged  species  rises 
steeply,  so  that  the  ratio  of  measured  concentration  to  LTE  concentration  reduces.  Related  work  in 
our  laboratory  has  shown  that  by  flowing  the  gas  more  quickly,  comparable  electron  densities  may 
be  achieved  with  lower  gas  temperatures.  Clearly,  additional  non-equilibrium  is  generated  in  the 
pulsed  discharge.  The  high  voltage  pulse  has  a  negligible  effect  on  the  gas  temperature  (and  hence 
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corresponding  LTE  number  density)  yet  the  measured  electron  number  density  in  the  discharge  in¬ 
creases  by  a  factor  of  at  least  4  immediately  following  the  high  voltage  pulse. 


Table  5  Ratio  of  the  measured  DC  electron  number  density  to  the  concentration  corresponding  to  an  LTE 
plasma  at  the  same  gas  temperature. 


i  (mA) 

Tg(K) 

ne-CRDS  /  ne-LTE 

52 

3100 

2.8xl04 

97 

3600 

980 

142 

4200 

48 

187 

4700 

5.6 

7.4.  NO+ Measurements 


7.4.1.  Experimental  Setup 

7.4.1  .a.  RF  Air  Plasma 

The  experimental  set-up  is  shown  schematically  in  Figure  1.  Atmospheric  pressure  air  plas¬ 
mas  are  generated  with  a  50  kW  radio-frequency  inductively-coupled  plasma  torch  operating  at  a 
frequency  of  4  MHz.  The  torch  is  operated  with  a  voltage  of  8.9  kV  and  a  current  of  4.6  A.  The 
torch  has  been  extensively  characterized  at  similar  conditions,  and  the  plasma  is  known  to  be  near 
LTE  with  a  temperature  of  about  7000  K  [1], 


7.4,  l.b.  CRD S  Measurements 

Unlike  the  N2+  ion,  the  NO+  ion  does  not  have  optically  accessible  electronic  transitions.  To 
perform  CRDS  measurements,  the  ion  must  be  probed  by  accessing  its  infrared  vibrational  transi¬ 
tions.  The  strongest  vibrational  transitions  are  the  fundamental  bands,  and  for  these  transitions  one 
finds  that  the  absorbance  per  NO+  ion  is  about  20000  times  less  than  that  of  the  electronic  transi¬ 
tions  of  the  N2+  ion.  Figure  79  shows  the  modeled  absorbance,  as  a  function  of  wavelength,  for  the 
air  plasma  at  the  conditions  used.  The  simulation  is  performed  with  Specair  and  assumes  a  path- 
length  of  5  cm,  and  LTE  conditions  at  a  temperature  of  7000  K  (Tg=Tr=Tv=Teiectronic=7000  K).  The 
simulation  includes  the  IR  absorption  features  of  NO,  OH,  and  NO+.  The  absorption  by  NO  and 
OH  is  relatively  weak,  while  the  various  fundamental  bands  ofNO+have  stronger  predicted  absorb¬ 
ances. 


It  is  evident  that  the  NO+  absorption  begins  at  a  wavelength  of  about  3950  nm,  and  is  a 
maximum  at  about  4100  nm.  Accessing  these  infrared  wavelengths  is  challenging  in  terms  of  avail¬ 
able  laser  sources.  The  current  measurements  have  been  performed  using  a  Continuum-Mirage 
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OPO  system.  The  Mirage  laser  is  designed  to  operate  at  a  maximum  wavelength  of  4000  nm,  how¬ 
ever  we  optimized  the  alignment  in  a  manner  that  enabled  operation  in  the  vicinity  of  4100  nm,  in 
order  to  be  nearer  to  the  peak  NO+  absorption.  Ring-down  cavity  alignment  at  these  wavelengths  is 
challenging,  since  the  beam  (and  its  back-reflections)  are  not  readily  observable.  The  ring-down 
cavity  was  aligned  using  a  combination  of  LCD  (liquid  crystal  display)  paper  to  locate  the  beam,  and 
a  Helium-Neon  laser  to  act  as  a  reference.  With  the  plasma  off,  ring-down  times  of  about  1.2  us 
were  obtained,  corresponding  to  mirror  reflectivities  of  about  0.998  (approximately  an  order  of 
magnitude  worse  than  the  mirrors  used  for  the  N2+  experiments). 
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Figure  79  Modeled  absorbance  of  the  air  plasma  at  LTE  temperature  of  7000  K,  and  path  length  of  5  cm.  Ab¬ 
sorption  by  NO,  OH,  and  NO+are  included.  Rotationally-resolved  lines  of  the  vibrational  transitions  are  shown. 

Our  initial  attempts  to  perform  CRDS  measurements  in  the  plasma  torch  used  the  same  cav¬ 
ity-geometry  as  was  used  in  the  N2+  experiments  -  a  g-parameter  of  0.5.  With  the  plasma  off,  this 
geometry  yielded  excellent  stability  in  the  ring-down  times:  1%  standard  deviation  in  ring-down 
time  for  single  shot  ring-down  signals.  However,  with  the  plasma  on,  the  beam  steering  reduced  the 
stability  significantly.  In  the  RF  plasma,  as  compared  to  the  smaller  nitrogen  plasma,  the  cavity- 
geometry  considerations  are  different.  In  the  smaller  nitrogen  plasma  we  wanted  to  simultaneously 
minimize  the  cavity  beam-waist,  and  the  beam-walk,  leading  to  a  g-parameter  of  —0.5  (see  discus¬ 
sion  above).  On  the  other  hand  in  the  RF  plasma,  the  plasma  dimension  (about  5  cm)  is  signifi¬ 
cantly  larger  than  the  beam  dimension  (about  1  mm).  Therefore,  the  exact  beam  dimension  is  not 
critical,  and  the  cavity-geometry  may  be  selected  solely  to  minimize  beam-walk.  The  numerical 
modeling  of  Spuler  et  al.  [91]  indicate  that  minimizing  the  beam-walk  may  be  accomplished  with  a 
g-parameter  of  about  0.25,  which  we  implemented  by  using  a  cavity  of  length  75  cm,  and  mirrors  of 
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radius-of-curvature  of  1  m.  This  geometry  did  indeed  reduce  the  beam-walk  and  enabled  improved 
stability  (about  2%  standard  deviation  in  empty  cavity  ring-down  times). 

As  will  be  discussed,  the  identification  of  spectral  lines  in  the  analysis  of  the  air  plasma 
spectra  is  challenging.  In  order  to  assist  in  identifying  NO+  spectral  features,  we  also  collected 
CRDS  spectra  with  the  plasma  running  with  argon  and  nitrogen  (as  opposed  to  ari),  conditions  that 
are  not  expected  to  have  any  significant  NO+  concentration. 


7.4.2.  Results  and  Discussion 

Figure  80  shows  a  measured  absorbance  spectrum  along  the  centerline  of  the  air  plasma. 
The  experimental  data  was  obtained  by  averaging  16  laser  shots  at  each  spectral  position.  The  plot¬ 
ted  CRDS  data  has  been  converted  to  absorbance,  and  fitted  with  a  peak-fitting  program.  (Fitted 
peaks  are  shown  in  black,  while  raw-data  is  shown  with  blue  symbols.)  Also  shown  is  the  modeled 
NO  absorbance  assuming  the  expected  plasma  conditions  of  path  length  5  cm,  and  LTE  at  7000  K. 
The  modeled  contributions  from  OH  and  NO  absorption  are  negligible  on  this  scale.  Comparing 
the  CRD  spectrum  in  the  air  plasma  to  the  CRD  spectrum  in  the  argon/nitrogen  plasma  provides 
information  as  to  line  identities.  The  largest  spectral  feature  (at  -4127.7  nm)  is  present  in  both 
spectra,  and  therefore  is  presumed  to  be  not  from  NO+.  Comparing  the  other  observed  spectral  fea¬ 
tures  with  the  model  does  not  yield  good  agreement.  To  the  best  of  your  knowledge,  the  spectro¬ 
scopic  constants  used  in  our  modeling  are  the  most  recent  and  accurate  ones  available  [94].  The 
exact  locations  of  the  rotationally  resolved  lines  are  largely  determined  by  the  B  parameter,  and  the 
quoted  uncertainty  on  the  B  parameters  in  this  work  is  +/-  0.005  cm'1  (or  about  0.25%).  Based  on 
the  quoted  uncertainty  we  performed  an  uncertainty  analysis,  and  found  that  with  this  level  of  pre¬ 
cision  it  is  not  possible  to  accurately  predict  the  locations  of  the  rotational  lines.  Therefore,  any 
match  between  the  experimental  data  and  model  would  be  fortuitous.  Our  experimental  features  are 
repeatable  (to  within  experimental  uncertainty)  and  have  approximately  the  correct  integrated  area, 
so  that  we  do  believe  they  are  due  to  NO+. 
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Figure  80  Experimental  and  modeled  absorbance  spectrum  from  the  air  plasma  in  the  vicinity  of  41 00  nm.  Raw 
data  (blue  symbols)  as  well  as  fitted  peaks  (top  black  line)  are  shown,  as  well  as  the  modeled  NO+  lines  (plotted 
negative  for  visual  clarity).  The  precision  of  the  spectroscopic  constants  used  in  the  model  is  insufficient  to  accu¬ 
rately  predict  the  location  of  the  rotational  lines. 


7.5.  Conclusions 

Spatial  and  temporal  profiles  of  N2+  concentration  have  been  measured  in  DC  and  pulsed 
atmospheric  pressure  nitrogen  glow  discharges  by  cavity  ring-down  spectroscopy.  Special  care  in 
the  selection  of  cavity  geometry  is  needed  in  the  atmospheric  pressure  plasma  environment.  Sub¬ 
millimeter  spatial  resolution,  microsecond  temporal  resolution,  and  sub-ppm  concentration  sensitiv¬ 
ity  have  been  achieved.  The  signal-to-noise  ratio  suggests  a  DC  detection  limit  of  about  7xl010  cm'3 
for  N2+  ions  at  our  experimental  conditions.  Using  a  collisional-radiative  model  we  infer  electron 
number  densities  from  the  measured  ion  profiles.  The  values  of  electron  number  density  found  in 
this  way  are  consistent  with  those  found  from  spatially  integrated  electrical  measurements.  The 
spectroscopic  technique  is  clearly  favorable,  because  it  offers  spatial  resolution  and  does  not  require 
knowledge  of  other  discharge  parameters.  The  measurements  verify  that  the  atmospheric  pressure 
discharges  studied  here  are  two-temperature  plasmas  with  electron  number  densities  elevated  with 
respect  to  equilibrium,  which  is  characteristic  of  glow-discharge  plasmas.  The  measurements  also 
demonstrate  the  feasibility  of  repetitively  pulsing  the  discharge  as  a  means  to  elevate  the  time- 
averaged  ionization  fraction. 

Measurements  of  the  NO+  ion  in  air  plasmas  have  also  been  demonstrated.  The  accessible 
spectral  features  of  NO+  are  vibrational  transitions,  are  considerably  weaker  than  the  UV  electronic 
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transitions  used  to  probe  N2+.  Nevertheless,  CRDS  data  from  air  plasmas  was  obtained,  and  spec¬ 
tral  features  attributed  to  NO+  were  observed.  More  precise  knowledge  of  the  spectroscopic  con¬ 
stants  of  NO+  is  needed  to  match  the  experimental  data  with  modeled  spectra. 
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8.  Summary  and  Conclusions 

Two-temperature  (Te>Tg)  kinetic  models  accounting  for  ionizational,  chemical,  vibrational 
and  electronic  nonequilibrium,  and  incorporating  a  collisional-radiative  model  with  over  11,000 
transitions  have  been  developed  to  understand  the  mechanisms  of  ionizational  nonequilibrium  in 
atmospheric  pressure  air  and  nitrogen  electrical  discharges.  These  models  predict  that  (even)  at  at¬ 
mospheric  pressure  energetic  electrons  driven  by  the  discharge  can  establish  and  maintain  electron- 
density  nonequilibrium  of  over  six  orders  of  magnitude.  An  unexpected  result  is  an  “S-shaped”  de¬ 
pendence  of  ne  on  Te  at  steady-state  for  a  given  gas  temperature.  This  behavior  results  from  a  tran¬ 
sition  between  predominately  molecular  ions  to  atomic  ions  at  a  critical  value  of  Te  and  values  of  ne 
above  about  1014  cm'3.  Above  this  critical  value  of  Te,  the  electron  density  increases  dramatically 
so  that  three-body  recombination  can  maintain  a  steady  state.  Departures  from  a  Maxwellian 
distribution  of  the  free-electrons  were  found  to  have  a  negligible  effect  on  the  predicted  steady-state 
characteristics,  at  least  for  the  case  of  a  nitrogen  plasma. 

The  feasibility  of  such  nonequilibrium  discharges  was  demonstrated  in  atmospheric- 
pressure  nitrogen  and  air  at  both  room  temperature  and  around  2000K  with  electrode  spacing  of  cm 
scale.  Stable,  diffuse  DC  discharges  have  been  achieved  at  atmospheric  pressure  for  a  range  of  gas 
flow  and  temperature  conditions  including  those  which  produce  ne  of  1012  to  1013  cm'3  without  sig¬ 
nificant  gas  heating.  Good  agreement  between  theoretical  and  measured  discharge  characteristics 
has  been  obtained  for  both  air  and  nitrogen  discharges  over  a  wide  range  of  conditions  including 
electron  densities  greater  than  1012  in  air  and  1013  in  nitrogen. 

To  reduce  the  power  required  to  maintain  such  nonequilibrium,  the  finite  electron  recombi¬ 
nation  time  (~10ps)  has  been  exploited  by  means  of  pulsed  discharges  of  10  ns  duration.  Both  sin¬ 
gle-shot  and  repetitively  pulsed  diffuse  discharges  at  100  kHz  have  been  demonstrated,  with  power 
reductions  of  over  two  orders  of  magnitude  for  average  electron  densities  greater  than  1012  cm'3. 
Power  reductions  of  approximately  three  orders  of  magnitude  are  possible  with  a  1  ns  repetitive 
pulsing  technique. 

A  computational  fluid  dynamics  code  for  the  simulation  of  flowing  nonequilibrium  air 
plasmas  including  the  presence  of  a  DC  discharge  was  developed  and  compared  to  the  DC  experi¬ 
ments  conducted  at  Stanford  University.  The  code  uses  a  detailed  two-temperature  chemical  kinetic 
mechanism,  along  with  appropriate  internal  energy  relaxation  mechanisms.  The  discharge  region 
was  modeled  by  generalizing  the  channel  model  of  Steenbeck,  and  a  new  semi-implicit  time  inte- 
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gration  method  was  developed  to  reduce  the  computational  cost.  The  results  of  a  simplified  ap¬ 
proximate  computational  model  are  in  general  agreement  with  the  experimental  data.  Initial  simu¬ 
lations  of  the  pulsed  discharges  have  been  carried  out. 

The  CRDS  technique  has  been  used  as  a  means  to  probe  the  population  of  the  N2+  ion. 
CRDS  allows  ground  state  species  measurements,  and  thus  offers  the  most  direct  possible  meas¬ 
urement  of  electron  number  density.  We  have  performed  spatially  resolved  measurements  in  a  DC 
atmospheric  pressure  nitrogen  discharge.  We  find  a  peak  concentration  of  2.6x1 012/cm3  and  a  col¬ 
umn  diameter  of  approximately  2.2  mm,  consistent  with  electrical  measurements.  We  have  also 
extended  the  CRDS  technique  to  perform  temporally  resolved  measurements  appropriate  for  diag¬ 
nostics  in  pulsed  plasmas.  Measurements  of  N2+  recombination  in  a  pulsed  nitrogen  discharge  have 
been  performed  with  a  temporal  resolution  of  1  microsecond.  The  measurements  show  significant 
(>~10)  enhancement  of  the  ionization  following  the  pulse,  and  indicate  that  recombination  occurs 
over  about  5  microseconds.  Measurements  of  the  NO+  ion  in  air  plasmas  have  also  been  demon¬ 
strated.  The  accessible  spectral  features  ofNO+are  vibrational  transitions,  are  considerably  weaker 
than  the  UV  electronic  transitions  used  to  probe  N2+.  Nevertheless,  CRDS  data  from  air  plasmas 
was  obtained,  and  spectral  features  attributed  to  NO+  were  observed.  This  technique  shows  great 
promise  for  the  measurement  of  NO+  concentrations  once  more  accurate  spectroscopic  constants  of 
NO+  become  available. 
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ber  2,  2000. 

•  Yalin,  A.P.,  Lommatzsch  U.,  Zare,  R.N.,  Laux,  C.O.,  Kruger,  C.H.,  “Cavity  Ring-Down  in  At¬ 
mospheric  Pressure  Plasmas,”  Western  Spectroscopy  Conference,  Jan  31.-Feb.2,  2001. 
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•  Laux,  C.O.,  Kruger,  C.H.,  and  Zare,  R.N.,  “Diagnostics  of  Atmospheric  Pressure  Air  Plasmas,” 
Arbeitsgemeinschaft  Plasmaphysik  (APP)  Spring  Meeting,  Bad  Honnef,  Germany,  February  18- 
21,  2001  (invited).. 

•  Yalin,  A.P.,  Lommatzsch  U.,  Zare,  R.N.,  Laux,  C.O.,  Kruger,  C.H.,  “Cavity  Ring-Down  Spec¬ 
troscopy  of  N2+  in  Pulsed  and  DC  Atmospheric  Pressure  Discharges,”  Proceedings  of  the  28th 
IEEE  International  Conference  on  Plasma  Science,  p.  9,  Las  Vegas,  NV,  June  17-22, 2001. 

•  Nagulapally,  M.,  Candler,  G.V.,  Laux,  C.O.,  Kruger,  C.H.,  “Numerical  Modeling  of  Repeti¬ 
tively  Pulsed  Discharges  in  Flowing  Atmospheric  Pressure  Plasmas,”  Proceedings  of  the  28th 
IEEE  International  Conference  on  Plasma  Science,  p.  351,  Las  Vegas,  NV,  June  17-22,  2001. 

•  Yu,  L.,  Packan,  D.M.,  Laux,  C.O.,  Kruger,  C.H.,  “Direct-Current  Glow  Discharges  in  Atmos¬ 
pheric  Pressure  Air  and  Nitrogen  Plasmas,”  Proceedings  of  the  28th  IEEE  International  Confer¬ 
ence  on  Plasma  Science,  p.  350,  Las  Vegas,  NV,  June  17-22,  2001. 

•  Kruger,  C.H.,  Laux,  C.O.,  Packan,  D.M.,  Yu,  L.,  Yalin,  A.P.,  Zare,  R.N.,  Nagulapally,  M.,  Can¬ 
dler,  G.V.,  Kelley,  J.D.,  “Nonequilibrium  Discharges  in  Atmospheric  Pressure  Air,”  Proceed¬ 
ings  of  the  28th  IEEE  International  Conference  on  Plasma  Science,  p.  348,  Las  Vegas,  NV,  June 
17-22,2001. 

•  Packan,  D.,  Yu,  L.,  Laux,  C.O.,  Kruger,  C.H.,  “Repetitively  Pulsed  DC  Glow  Discharge  in  At¬ 
mospheric  Pressure  Air:  Modeling  and  Experiments  with  a  12  kV,  10  ns,  100  kHz  Pulse  Gen¬ 
erator,”  Proceedings  of  the  28th  IEEE  International  Conference  on  Plasma  Science,  p.  259,  Las 
Vegas,  NV,  June  1 7-22,  2001. 

•  Kruger,  C.H.,  Laux,  C.O.,  Yu,  L.,  Packan,  D.M.,  Pierrot,  L.,  “Nonequilibrium  Discharges  in  Air 
and  Nitrogen  Plasmas  at  Atmospheric  Pressure,”  15th  International  Symposium  on  Plasma 
Chemistry,  Orleans,  France,  July  9-13, 2001  (invited). 

•  Laux,  C.O.,  “DC  and  Pulsed  Discharges  in  Atmospheric  Pressure  Air  and  Nitrogen  Plasmas,” 
43rd  Annual  Meeting  of  the  Division  of  Plasma  Physics,  Long  Beach,  CA,  Oct.  29-  Nov.  2, 

2001.  Bulletin  of  the  American  Physical  Society,  Vol.  46,  No.  8,  pp.  99-100,  October  2001  (in¬ 
vited). 

•  Laux,  C.O.,  “Air  Plasma  Ramparts  Research,”  ONERA/DGA/  Industry  Workshop,  ONERA, 
Palaiseau,  France,  March  27,  2002  (invited). 

•  Kelley,  J.D.,  Packan,  D.M.,  Yu,  L.,  Laux,  C.O.,  and  Kruger,  C.H.,  “Repetitively  pulsed  glow 
discharge  in  atmospheric  pressure  air:  modeling  and  experiments,”  Fourth  Workshop  on  Mag¬ 
neto-  and  Plasma  Aerodynamics  for  Aerospace  Applications,  Moscow,  Russia,  April  9-11, 

2002. 

•  Laux,  C.O.,  “Nonequilibrium  Air  Plasmas  for  Scramjet  Ignition,  Bio-Decontamination,  and 
other  Engineering  Applications,”  Thermosciences  Seminar  Series,  Stanford  University,  April 
17,2002. 

•  Laux,  C.O.,  “Nonequilibrium  Air  Plasmas  for  Scramjet  Ignition,  Aerodynamic  Flow  Control, 
Bio-Decontamination,  and  other  Engineering  Applications,”  Joint  Electrical  Engineer¬ 
ing/Mechanical  Engineering  Seminar,  Old  Dominion  University,  April  26, 2002. 

•  Chauveau,  S.M.,  Laux,  C.O.,  Kelley,  J.D.,  and  Kruger,  C.H.,  “Vibrationally  Specific  Colli- 
sional -Radiative  Model  for  Nonequilibrium  Air  Plasmas,”  AIAA  2002-2229,  33rd  AIAA  Plas- 
madynamics  and  Lasers  Conference,  Maui,  HW,  May  20-23,  2002. 

•  Yalin,  A.P.,  Laux,  C.O.,  Kruger,  C.H.,  and  Zare,  R.N,  “Spatially  and  Temporally  Resolved 
Concentration  Measurements  of  the  N2+  Ion  in  Nitrogen  Discharges  by  Cavity  Ring-Down 
Spectroscopy,”  AIAA  2002-2245,  33rd  AIAA  Plasmadynamics  and  Lasers  Conference,  Maui, 
HW,  May  20-23,  2002. 
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•  George,  C.,  Candler,  G.V.,  Laux,  C.O.,  and  Kruger,  C.H.,  “Computational  Analysis  of  Diffuse 
Discharges  in  Atmospheric  Pressure  Air,”  AIAA  2002-2223,  33rd  AIAA  Plasmadynamics  and 
Lasers  Conference,  Maui,  HW,  May  20-23,  2002. 

•  Laux,  C.O.,  “Radiation  and  Nonequilibrium  Collisional-Radiative  Models,”  von  Karman  Insti¬ 
tute  Special  Course  on  Physico-Chemical  Modeling  of  High  Enthalpy  and  Plasma  Flows, 
Rhode-Saint-Genese,  Belgium,  June  4-7, 2002  (invited). 

•  Machala,  Z.,  Laux,  C.O.,  and  Kruger,  C.H.,  “Nonequilibrium  Microwave  and  DC-Discharge 
Plasmas  in  Atmospheric  Pressure  Air,”  poster  presented  at  the  Gordon  Research  Conference  on 
Plasma  Processing  Science,  Tilton,  NH,  July  21-26,  2002. 

•  Yalin,  A.P.,  Machala,  Z.,  Laux,  C.O.,  Kruger,  C.H.,  and  Zare,  R.N.,  “Optical  Emission  and  Cav¬ 
ity  Ring-Down  Spectroscopy  Measurements  in  Atmospheric  Pressure  Nitrogen  Glow  Dis¬ 
charge,”  poster  presented  at  the  Gordon  Research  Conference  on  Plasma  Processing  Science, 
Tilton,  NH,  July  2 1  -26,  2002. 

•  Yalin,  A.P.,  Laux,  C.O.,  Kruger,  C.H.,  and  Zare,  R.N.,  “Temporally  and  Spatially  Resolved 
Cavity  Ring-Down  Spectroscopy  Measurements  in  Atmospheric  Pressure  Plasmas,”  OSA  An¬ 
nual  Meeting  and  Exhibit  2002  and  Laser  Science  XVIII,  Orlando,  FL,  Sept.  29-  Oct.3,  2002. 

•  Laux,  C.O.,  Machala,  Z.,  Duten,  X.,  Packan,  D.M.,  Yu,  L.,  Kruger,  C.H.,  “Scaled-up  Nonequi¬ 
librium  Air  Plasmas,”  to  be  presented  at  the  5th  Weakly  Ionized  Gases  AIAA  Workshop,  Reno, 
NV,  January  6-9  2003. 

•  Chauveau,  S.M.,  Kelley,  J.D.,  Laux,  C.O.,  Kruger,  C.H.,  “Vibrationally  Specific  Collisional- 
Radiative  Model  for  Nonequilibrium  Air  Plasmas,”  to  be  presented  at  the  41st  AIAA  Aerospace 
Sciences  Meeting  and  Exhibit,  Reno,  NV,  January  6-9,  2003. 


11.2.  Consultative  and  Advisory  Functions  to  Other  Laboratories  and  Agencies 

•  Prof.  Candler  currently  serves  on  the  DoE/Sandia  ASCI  Advisory  Panel.  He  is  a  member  of  the 
NATO  Research  Technology  Organization,  Working  group  10,  Technologies  for  Hypersonic 
Flight.  He  also  serves  on  the  Thermophysics  committee  of  the  American  Institute  of  Aeronau¬ 
tics  and  Astronautics. 

•  Dr.  Laux  serves  on  the  Plasmadynamics  and  Lasers  Committee  of  the  American  Institute  of 
Aeronautics  and  Astronautics. 


12.  New  Discoveries,  Inventions,  or  Patent  Disclosures 

None 


13.  Honors/Awards 
Prof.  C.H.  Kruger’s  group: 

•  1997  Invited  Plenary  Lecture:  Kruger,  C.H.,  Owano,  T.G.,  Laux,  C.O.,  “Nonequilibrium  Plas¬ 

mas  and  Diagnostics,”  23rd  International  Conference  on  Phenomena  in  Ionized  Gases,  Toulouse, 
France,  July  17-22,  1997. 
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•  1998  Invited  Lecture:  Laux,  C.O.,  Gessman,  R.J.,  Packan,  D.M.,  Yu,  L.,  Kruger,  C.H.,  and  R.N. 
Zare,  “Experimental  Investigations  of  Ionizational  Nonequilibrium  in  Atmospheric  Pressure  Air 
Plasmas,”  25th  IEEE  International  Conference  on  Plasma  Science  (ICOPS),  Raleigh,  NC,  June 
1-4,  1998. 

•  1998  Best  Paper  Award,  AIAA  Plasmadynamics  and  Lasers:  Gessman,  R.J.,  Laux,  C.O.,  and 
Kruger,  C.H,  “Experimental  Study  of  Kinetic  Mechanisms  of  Recombining  Atmospheric  Pres¬ 
sure  Air  Plasmas,”  award  presented  at  the  29th  AIAA  Plasmadynamics  and  Lasers  Conference, 
Albuquerque,  NM,  June  1998. 

•  1999  Invited  Lecture:  Laux,  C.O.,  Yu,  L.,  Packan,  D.M.,  Gessman,  R.J.,  Pierrot,  L.,  Kruger, 
C.H.,  and  Zare,  R.N.,  “Ionization  Mechanisms  in  Two-Temperature  Air  Plasmas,”  AIAA  99- 
3476,  30th  AIAA  Plasmadynamics  and  Lasers  Conference,  Norfolk,  VA,  June  28-July  1,  1999. 

•  2000  Invited  Lecture:  Laux,  C.O.,  "Glow  Discharges  in  Atmospheric  Pressure  Air,"  Gordon 
Research  Conference  on  Plasma  Processing  Science,  Tilton,  NH,  August  13-18,  2000. 

•  2001  Invited  Plenary  Lecture:  Laux,  C.O.,  Kruger,  C.H.,  and  Zare,  R.N.,  “Diagnostics  of  At¬ 
mospheric  Pressure  Air  Plasmas,”  Arbeitsgemeinschafit  Plasma  Physik  (APP)  Spring  Meeting, 
Bad  Honnef,  Germany,  February  18-21,  2001. 

•  2001  Invited  Plenary  Lecture:  Kruger,  C.H.,  Laux,  C.O.,  Yu,  L.,  Packan,  D.M.,  Pierrot,  L., 
“Nonequilibrium  Discharges  in  Air  and  Nitrogen  Plasmas  at  Atmospheric  Pressure,”  15th  Inter¬ 
national  Symposium  on  Plasma  Chemistry,  Orleans,  France,  July  9-13,  2001. 

•  2001  Invited  Lecture:  Laux,  C.O.,  “DC  and  Pulsed  Discharges  in  Atmospheric  Pressure  Air  and 
Nitrogen  Plasmas,”  43rd  Annual  Meeting  of  the  Division  of  Plasma  Physics,  Long  Beach,  CA, 
Oct.  29-  Nov.  2,  2001.  Bulletin  of  the  American  Physical  Society,  Vol.  46,  No.  8,  pp.  99-100, 
October  2001. 

•  2002  Invited  Special  Course:  Laux,  C.O.,  “Radiation  and  Nonequilibrium  Collisional-Radiative 
Models,”  von  Karman  Institute  Special  Course  on  Physico-Chemical  Modeling  of  High  En¬ 
thalpy  and  Plasma  Flows,  Rhode-Saint-Genese,  Belgium,  June  4-7,  2002. 


Prof.  R.N.  Zare: 

•  1997  Eastern  Analytical  Symposium  Award  for  Outstanding  Achievements  in  the  Field  of  Ana¬ 
lytical  Chemistry. 

•  1 997  California  Scientist  of  the  Year. 

•  1 997  NASA  Exceptional  Scientific  Achievement  Award. 

•  1997  Fischer  Award  in  Analytical  Chemistry  given  by  the  American  Chemical  Society. 

•  Allan  V.  Cox  Medal  for  Faculty  Excellence  Fostering  Undergraduate  Research,  Stanford  Uni¬ 
versity,  1998. 

•  ACS  Award  in  Analytical  Chemistry,  1998. 

•  Centennial  Medal,  Graduate  School  of  Arts  and  Sciences,  Harvard  University,  1998. 

•  E.B.  Wilson  Award  in  Spectroscopy,  American  Chemical  Society,  1999. 

•  Honorary  Doctor  of  Science,  Columbia  University,  1 999. 

•  The  Welch  Award  in  Chemistry,  1999. 

•  Doctor  of  Philosophy,  honoris  causa,  Uppsala  University,  2000. 

•  Arthur  L.  Schawlow  Prize  in  Laser  Science,  American  Physical  Society,  2000. 

•  Nobel  Laureate  Signature  Award  for  Graduate  Education,  American  Chemical  Society,  2000. 
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•  Charles  Lathrop  Parsons  Award,  American  Chemical  Society,  2001,  "to  recognize  outstanding 
public  service  by  a  member  of  the  American  Chemical  Society." 

•  Faraday  Medal  &  Lectureship,  Royal  Society  of  Chemistry,  2001. 

•  2000-2001  Madison  Marshall  Award,  American  Chemical  Society,  North  Alabama,  2001. 

•  Honorary  Doctorate,  State  University  of  West  Georgia,  2001. 

•  Honorary  Doctorate,  D  Univ,  The  University  of  York,  2001. 

•  2000-200 1  CaSSS  (California  Separation  Science  Society)  Scientific  Achievement  Award,  2002. 

•  Distinguished  Chemist  Award,  American  Chemical  Society  Sierra  Nevada  Section,  2002. 


Prof.  G.V.  Candler’s  group: 

•  1998  Invited  Plenary  Paper:  Candler,  G.V.,  “High  Enthalpy  Flow  Simulation  Challenges,”  29th 
AIAA  Plasmadynamics  and  Lasers  Conference,  Albuquerque,  NM,  June  1998. 

•  1 999  Invited  Paper:  Candler,  G.  V.  and  Kelley,  J.  D.,  “Supersonic  Flow  of  Discharge-Excited 
Air  Around  Blunt  Bodies,”  Workshop  on  Perspectives  of  MHD  and  Plasma  Technologies  in 
Aerospace  Applications,  Moscow,  Russia,  March  24-25, 1999. 

•  1999  Invited  paper:  Campbell,  C.H.,  and  G.V.  Candler,  “Detailed  Simulation  of  Nitrogen  Dis¬ 
sociation  in  Shock  Waves,”  33rd  AIAA  Thermophysics  Conference,  Norfolk,  VA,  June  1999. 

•  2001  Best  Student  Paper  Award,  AIAA  Thermophysics  Committee. 

•  2002  George  Taylor  Distinguished  Research  Award:  Candler,  G.V. 


Dr.  J.  D.  Kelley: 

•  September  1998:  named  a  Boeing  Technical  Fellow. 

•  1999  Invited  Paper:  Candler,  G.  V.  and  Kelley,  J.  D.,  "Supersonic  Flow  of  Discharge-Excited 
Air  Around  Blunt  Bodies,"  Workshop  on  Perspectives  of  MHD  and  Plasma  Technologies  in 
Aerospace  Applications,  Moscow,  Russia,  March  24-25, 1999. 
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